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PREFACE 

This  work  was  carried  out  as  a part  of  an  In-House  Laboratory  Independent  Research 
project  entitled,  “Investigation  of  the  Strength  and  Stability  Characteristics  of  Pressurized 
Rib  Structures."  This  work  was  initiated  in  response  to  the  findings  of  a systems  analysis 
which  identified  the  pressurized  rib  concept  as  being  the  most  promising  structural 
alternative  for  meeting  the  Army  requirement  for  lightweight  highly  mobile  tentage.  As 
is  indicated  by  the  references  cited  in  this  report,  two  previous  reports  describing  a similar 
investigation  of  pressure-stabilized  beams  have  been  published.  In  the  reference  citation 
the  organizations  "US  Army  Natick  Laboratories"  and  "US  Army  Natick  Development 
Center"  refer  to  the  organization  now  called  the  "US  Army  Natick  Research  and 
Development  Command".  The  author  wishes  to  express  his  thanks  to  Colette  Klarman, 
formerly  of  the  Data  Analysis  Office,  for  her  effort  in  writing  the  computer  program 
for  carrying  out  the  solution  of  the  equations  developed  in  this  report  and  for  plotting 
the  theoretical  and  experimental  results. 
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THE  STRUCTURAL  BEHAVIOR  OF  PRESSURE-STABILIZED  ARCHES 


INTRODUCTION 

The  use  of  pressure-stabilized  structures  for  Army  field  shelters  has  been  the  subject 
of  considerable  attention  over  the  past  decade  with  emphasis  on  the  single-wall  and 
double-wall  type  of  shelter,  as  indicated  by  references  1 and  2.  A systems  analysis  of 
Army  field  shelter  needs  for  the  1985  time  frame  revealed  that  the  requirements  for 
lightweight  tentage  of  low  package  bulk  and  with  minimum  setup  and  disassembly  times 
might  be  most  effectively  achieved  with  a tent  having  a frame  of  highly  pressurized 
(compared  with  present  air-supported  tents)  structural  elements  supporting  a lightweight 
fabric  barrier,  as  illustrated  in  Figure  1.  As  a result  of  the  potential  indicated  by  the 
systems  analysis  an  investigation  of  the  structural  behavior  of  pressure-stabilized  structural 
elements  was  begun  with  the  objective  being  to  develop  a design  capability  which  would 
permit  the  concept  to  be  evaluated  for  use  in  Army  tentage.  The  results  of  a theoretical 
and  experimental  study  of  the  behavior  of  pressure-stabilized  beams  under  load  are 
presented  in  references  3 and  4.  However,  as  is  illustrated  in  Figure  1,  pressure-stabilized 
arches  will  be  used  in  this  concept  for  tentage  and  thus  design  data  for  these  structural 
elements  are  required.  Such  design  data  are  presented  in  this  report  in  the  form  of  a 
theory  for  the  deformation  of  pressure-stabilized  arches  together  with  experimental 
confirmation  of  the  prediction  made  with  this  theory. 

A review  of  previous  work  on  the  behavior  of  pressure-stabilized  structural  elements 
is  given  in  references  1 and  2.  The  only  additional  work  known  to  the  author  of  interest 
relative  to  the  behavior  of  arches  is  the  theoretical  study  of  the  stability  of 
pressure-stabilized  rings  given  in  reference  5.  This  work  dealt  entirely  with  the  stability 
problem  for  the  complete  ring  and  these  results  were  not  used  in  the  present  work  because 
it  was  believed  that  a more  consistent  theory  would  be  obtained  by  adopting  the  ideas 
of  superposition  of  small  displacements  on  larger  ones  as  was  done  for  the 
pressure-stabilized  beam  problem  in  reference  3. 


Dietz,  A.  E.,  R.  B.  Proffitt,  R.  S.  Chabot,  E.  L.  Moak,  and  C.  J.  Monego;  Wind 
Tunnel  Tests  and  Analyses  for  Ground-Mounted,  Air-Supported  Structures  (revised); 
US  Army  Natick  Laboratories  Technical  Report  70-7-GP;  1969. 

Dietz,  A.  E.,  R.  B.  Proffitt,  R.  S.  Chabot,  and  E.  L.  Moak;  Design  Manual  for 
Ground-Mounted  Air-Supported  Structures  (Single-and  Double-Wall)  (revised);  US 
Army  Natick  Laboratories  Technical  Report  69-59— GP;  1969.  ,* 

Steeves,  Earl  C.;  A Linear  Analysis  of  the  Deformation  of  Pressure  Stabilized  Beams; 
US  Army  Natick  Laboratories  Technical  Report  75-47-AMEL;  1975.  r 

Steeves,  Earl  C.;  Behavior  of  Pressure  Stabilized  Beams  Under  Load;  US  Army  Natick 
Development  Center  Technical  Report  75-82-AMEL;  1975.  , 

Weeks,  G.  E.;  Buckling  of  a Pressurized  Toroidal  Ring  Under  Uniform  External 
Loading;  NASA  TN  D-4124;  1967. 
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Included  in  this  report  is  the  derivation  of  the  equations  describing  the  behavior 
of  pressure-stabilized  arches  under  loads  both  in  the  plane  of  the  arch  and  normal  to 
the  plane  of  the  arch.  The  solution  of  the  equations  for  the  in-plane  problem  are  presented 
along  with  a computer  program  to  carry  out  this  solution  numerically.  The  experimental 
study  is  also  descirbed,  including  the  loading  tests  on  arches  and  the  measurement  of 
the  elastic  properties  of  the  fabric  from  which  the  arches  were  made.  The  results  of 
this  experimental  study  are  compared  with  the  predictions  of  the  theory. 

THEORETICAL  ANALYSIS 

This  section  of  the  report  is  concerned  with  development  of  the  theory  to  predict 
the  behavior  of  pressure-stabilized  arches  under  load.  We  begin  with  a derivation  of  the 
governing  equations  for  deformation  in  the  plane  of  the  arch  and  normal  to  the  plane 
of  the  arch.  These  sets  of  equations  are  uncoupled  and  the  solution  to  the  in-plane  set 
is  obtained  here  and  a computer  program  to  carry  out  the  solution  numerically  is  described. 

Derivation  of  Governing  Equations 

Fundamental  Principles 

We  consider  a surface  in  the  form  of  a circular  torus  with  a circular  cross-section, 
as  shown  on  Figure  2.  Points  on  the  surface  are  located  by  specification  of  the  coordinates 
0,  and  02,  and  the  differential  arc  element  and  radii  of  curvature  on  the  surface  are 
given  by 


ds2  = [R  + acos(0, )]  2d02  + a2d02 

da) 

1/Ri  = 1/a 

(1b) 

1/R2  = COS(0,  )/ [3005(0, ) + R] 

(1c) 

In  equation  (1),  a is  the  cross-section  radius  and  R is  the  radius  of  the  torus  centerline. 
Assuming  that  the  membrane  state  of  stress  is  present,  the  principle  of  minimum  potential 
energy  is  written  as 

5 fd^HNne,,  + N22e22  + 2N12e12  - 2Pv3]dA  = 0 <2> 

Where  Njj,  ejj  p an(j  are,  respectively,  the  stress  resultants,  the  surface  strains,  the 
inflation  pressure,  and  the  normal  displacement.  The  surface  strains  can  be  expressed 
in  terms  of  the  displacements  and  their  derivatives  for  the  toroidal  coordinates  as: 

Cm  = (1/a)[v,  , + v3  + (1/2a)(-v3jl  + v,)2]  (2a) 
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e22  * (1/a2)[v2(2  - v,sin(0, ) + v3cos(0,)  + (1/2a2)(-v3  2 + v2cos(0,))2] 

*i  2 = (1/2a2)[(R/a  + cos(0I))v2  l + vl  2 + v2sin(0,)  + 

(1/a)(-v3  , + v1)(-v3  2 + v2cos(0 1 ))  ] 


(2b) 


(2c) 


In  equations  (2)  a2  - (R  + acos(0,)).  To  complete  the  basis  for  the  theory  we  specify 
an  orthotropic  law  relating  the  stress  resultants  and  the  surface  strains  as 


Nii  = C] i€i  i + Ci 2e22 
N22  = Ci2ei!  + C22e22 
Nj2  = 2C336i2 


* (3a) 
(3b) 
(3c) 


It  has  been  assumed  in  writing  equations  (3)  that  the  lines  of  orthotropic  symmetry  are 
parallel  to  the  coordinate  lines. 

Energy  Principle  for  Deformation  about  Pressurized  State 

The  derivation  now  proceeds  identically  with  that  for  pressure-stabilized  beams  in 
reference  3 by  assuming  that  the  displacements  vj  can  be  expressed  as  the  sum  of 
displacements  Wj  due  to  pressurization  and  Uj  due  to  applied  load  as 


Vj  = Wj  + Uj 


(4) 


Substitution  of  this  displacement  assumption  into  the  strain-displacement  relation  (2)  yields 

(5) 


tjj  = ejj  + ejj  + ey 


i = 1,2 
j = 1,2 

These  supplementary  strain  measures  are  defined  in  terms  of  the  displacements  Wj  and 
Uj  as 

e°i  = (l/altw,!  + w3  + (1/2a)(-w3jl  + w,)J] 

eii  = (1/a)[u,i  + u3  + (1/a)(— w3l  + w,)(-u3(,  + u2 )] 

, , ' ' (6) 
eri  - (1/2aJ)(-u3fI  + u,)J 

elz  =(1/or2)fw2<2  - w,  sin(0, ) +w3cos(0,)  + (1/2a2)(-w32  + w2cos(0,))2] 
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e22  ~ (1/<*2)[u2(2  - u,sin(0,)  + u3cos(0,)  + 

(1/«j)(-wj#a  + WjCosWJK-us  j + u2oot(0,Ul 

e22  ~ (1/2o^)[-u3j2  + u2cos(0,)]J 

e°2  =(1/2a2K(R/a  + cos(0,))w2(I  + Wi2  + w2sin{0,)  + 

(1/aH-Wj  , + wx)(-w32  + w2 cos(0 1 ))]  (6) 

e,'2  = (1/2ck2  ) [(R/a  + cos(0,))u2i  + u1>2  + u2sin(0|)  + 

(1/a)(-w3  , + w,  )(-u3  2 + u2cos(0,))  + 

(1/a)(-u3f,  + u,)(-w32  + w2  cos(0 , ))] 

e,"2  = (1/2aa2)[— u3(,  + u,][-u3(J  + UlCos(0,)] 

The  ejj  are  the  complete  nonlinear  strain  measures  for  the  pressurized  state  expressed 
in  terms  of  the  displacements  Wj.  The  ejj  are  the  linear  strain  measures  in  the  displacements 
u;  due  to  applied  load  and  also  contain  the  terms  coupling  the  pressurized  state  and  that 
due  to  applied  load.  The  e"j  are  the  nonlinear  terms  of  the  strains  due  to  applied  load 
and  they  can  be  interpreted  as  the  products  of  the  rotations  of  a surface  line  element 
about  the  i^1  and  j^1  axes.  Corresponding  to  these  supplementary  strain  measures  is  a 
set  of  stress  measures  which  represent  the  total  stress  components  as 


Ni|  = Njj  + Njj  + N 


i = 1.2 

i = 1.2 


Expressions  for  the  stress  measures  Njj,  Njj,  and  Njj  in  terms  of  the  strain  measures  ejj, 
ejj,  and  e,j  are  found  by  substitution  of  equation  (5)  into  the  stress  strain  law  (3)  which 


yields 


~ Ci , e,  | 

+ Ci  2e22 

- t 

~ C(  i Ci  | 

+ Ci  2e22 

= Ci iCi , 

+ Ci2e22 

= C| 2ef i 

+ C22e22 

= C, 2e\ i 

+ C22e22 

= C| 26i 1 

+ C22e22 
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(8) 


N i 2 = 2C33ei2 

Nia  * 2C33e,3 

Nij  = 2C336|2 

Substitution  of  equations  (4),  (5)  and  (7)  into  the  energy  principle  (2)  and  retaining  terms 
up  to  those  quadratic  in  the  displacement  components  Uj  yields 

6 Jl/2[(e°X,  + e2jl^2  + 2e°i  jN° 2 - 2Pw3)  + 2(e'nN?,  + + 

2^12^12  — Puj)  + 1 1 N j 3 + €2 2 N*2 2 + 2e  1 2 2 + 2e'Iihf>ii  + 

2^22^22  ^€12^12  — 2FjUi  — 2F2U2  — 2F3u3)]dA  ■ 0 (9) 


In  writing  this  relationship  the  work  due  to  the  applied  forces  F;  which  do  work  only 
through  the  displacement  components  uj  (i  = 1,3)  has  been  added,  and  the  following 
identities  which  can  be  proven  through  the  stress-strain  law  have  been  used 


f i° , n ; , 

+ 

ea  2 Na  2 

= 

cI.Nm 

+ 

C*2  2 N1 

cl  2 N 1 2 

= 

Ci  2 Ni  2 

euNn 

+ 

eajNaa 

* 

cr.wt. 

+ 

C2jN 

ei  2 N|  2 

= 

cr2N°2 

(10) 


To  yield  the  governing  equation  for  deformation  about  the  pressurized  state,  which  is 
described  by  the  displacement  components  uj,  the  variation  in  equation  (9)  must  be 
interpreted  so  as  to  be  carried  out  with  respect  to  the  components  uj.  Since  the  expression 
inclosed  in  the  first  set  of  parenthesis  in  equation  (9)  is  independent  of  the  components 
Uj,  its  variation  must  vanish.  Carrying  out  the  variation  of  the  expression  in  the  second 
set  of  parenthesis  and  performing  the  required  integration  by  parts,  one  finds  that  the 
coefficients  of  the  variational  displacements  5u;  are  the  equilibrium  equations  governing 
the  pressurized  state  which  must  vanish  if  the  stress  resultants  N , , , f/a  a , and  N , a represent 
an  equilibrium  state.  Taking  these  facts  into  account,  the  following  is  the  energy  principle 
for  deformation  about  the  pressurized  state: 

6 f 1/2(«', , N;  , + + 2e;aN'ia  + 2e ” 1 N°  , + 2e$aN$,  + 

A . 0 (11) 

^12^11  2F,U|  2Fauj  2F3u3)dA  * 0 

00  o 

Since  the  stresses  N, , , Nya,  and  N,  a due  to  pressurization  are  considered  known  quantities 
like  the  elastic  constants,  this  energy  expression  is  quadratic  in  the  displacements  u;  and 
will  yield  a set  of  linear  equations.  Thus  the  displacement  assumption  (4)  has  provided 
a means  for  linearizing  the  problem. 
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Simplified  Displacement  Approximation 

As  indicated,  we  have  been  able  to  linearize  the  problem  while  including  the  effect 
of  pressurization  but  we  still  have  a problem  in  two  independent  coordinates  which  will 
be  described  in  terms  of  a set  of  linear  partial  differential  equations.  The  present  objective 
is  to  simplify  the  mathematical  form  by  introducing  a displacement  approximation  which 
specifies  the  form  of  the  displacement  components  with  respect  to  the  meridional 
coordinate,  0j . This  displacement  approximation  describes  the  deformation  in  terms  of  rigid 
body  motions  of  the  cross-section;  that  is,  the  cross-section  retains  its  shape  during 
deformation.  The  approximation  is  the  one  used  in  beam  theory.  Here,  however,  transverse 
shear  deformation  is  included.  This  approximation  is  stated  mathematically  by  expressing 
the  displacement  components  u;  in  terms  of  the  three  displacements  of  the  cross-section 
U,  V,  W and  the  rotations  of  the  cross-section  about  three  mutually  perpendicular  axes, 
#x>  ^y>  as: 


Ui 

= a^x  + Wsin(0,)  + Vcos(0,) 

(12a) 

= U - tf>yacos(0,)  - 0zasin(0,) 

(12b) 

U3 

= Vsin(0,)  — Wcos(0,) 

(12c) 

These  cross-section  displacements  which  are  functions  of  02  only  are  shown  graphically 
in  Figure  2.  These  displacement  approximations  can  be  used  to  formulate  one-dimensional 
strain  measures,  stress  resultants,  and  a one-dimensional  energy  principle. 

The  one-dimensional  strain  measures  are  obtained  by  substituting  equations  (12)  into 
equations  (6)  yielding 

e.'.  * 0 

el,  = *’/2  (13) 

e22  = (1/a2)[U'  - W)  - 0yacos(0,)  - (pz  + px)asin(0,)] 
e2"2  = (1/2a2)[(U  + W')cos(0,)  - Vsin(0,)  -^vacosI(0,)  - pzasin(0,  )cos(0,  )J 2 
*;2  = (1/2a2  )[a(0x  -4>z)  + (W'  + U + Rpy)sin(0,)  + (V'  - Rpz)cos(0,)] 
e,’2  = (1/2a2)[px((U  + W'  )cos(0, ) - V'sin(0,)  - Pyacos^O.)  - 

^asinlO,  )cos(0, ))] 

At  this  point  we  introduce  two  additional  assumptions.  The  first  of  these  assumes  that 
the  deformation  is  such  that  a uniaxial  state  of  strain  results  which  implies  ctl  = 0. 
This  is  consistent  with  the  displacement  approximation  (12)  to  the  first  order  end  requires 
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that  we  neglect  the  second  order  term  <p\  and  set  e"t  • 0.  The  second  of  these  additional 
assumptions  is  that  the  shear  strain  can  be  adequately  represented  by  a linear  expression 
in  the  displacements  and  their  derivatives.  This  assumption  implies  that  2 - 0.  This 
assumption  can  be  further  justified  by  noting  that  in  the  energy  expression  (11)  e'2  is 
a coefficient  of  N,2  which  is  the  shear  stress  component  for  the  pressurized  state.  Since 
this  stress  component  vanishes  for  internal  pressure  loading,  the  nonlinear  term  being 
neglected  in  the  shear  strain  would  make  no  contribution  to  the  energy. 

The  one-dimensional  energy  principle  is  obtained  by  incorporating  the  above 
assumptions  in  the  energy  principle  (11)  and  rewriting  it  as  a quadratic  in  the  strain 
measures  as: 

6 f 1/2  [C22(ei2)2  + 4C33(e,'2)2  + 2N^2e2"2  (14) 

-2(F,U|  + F2u2  + F3u3)]dA  = 0 

With  the  substitution  of  equations  (12)  and  (13)  into  this  expression  and  integrating  with 
respect  to  6,  from  0 to  2n  the  one-dimensional  energy  principle  is  obtained  as 

6 f */ 2[^[2(U'  - W)2  + Z2(i)2(U'  - W + R^)2  + a2Z3(02  + 0X>J)  + 
n R 

-Hi-3[2a2(^  - 0Z)2  + Z2(— )2  (a(0x  - *z)  - £ <V'  - R0Z))2  + 

R2  R 

o 

Z3(W'  + U + Rtfy)2)  + -^[(U  + W')2  + Z(1)2(U  + W'  + R*y)2  + 

(V')2  + Z4(±)2(V'  - R0Z)2]  - 2(Uf2  + Wf,  + Vf3  + 

R 

a0xf4  - a^yf5  - a0zf6]]aRd02  =0  (15) 

The  details  of  this  development  are  presented  in  Appendix  A. 

As  derived  in  Appendix  B,  expressions  for  the  one-dimensional  stress  resultants  and 
moments  in  terms  of  the  cross-section  displacements  and  rotations  are 

t = irC22ll<2  + <|)2Z2)(U'  - W)  + (1)2Z2R^1 

R R (16a) 

qy  = JrC33iz2[(V'  - R0Z)  - R(i)2(*'x  - 0Z)]  - 

^ajl-  Z3  V'  + (i  )2  RZ40z]  (16b) 
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(16c) 


L 


q.  = *C33Z31(W'  + U + R*y]  + 
z R 1 

7rN22|[Z2(U  + W’)  + (I^ZR^yl 
H H 

mx  = *C33a(i)2[R(2  + (1)2Z 2)(0X  - 4>z)  - Z2(V'  - R <f>z)]  (16d) 

R R 

my  = — jrC22Z2a(|)2[U'  - W + R0y) 
mz  = — ffC22Z3aR(  — )2  (0Z  + <px) 

n 

The  energy  principle  (15)  can  be  expressed  in  terms  of  these  stress  resultants  and  moments 
by  carrying  out  the  variation  and  substituting  equations  (16)  to  give 

t^t6(U'  - W)  - 1 my<Py  — lmz6(*z  + <f>x)  + 

R R 

imx6(*;  - <t>z)  + — qy6(V'  - Rtfz)  + *N2  2 Iv'6tR0z) 

R R R R 

+ -qz5(U  + W'  R«y)  --N22-|(U  + W')6(R0y)  (17) 

R R R 

— 7ra(f25U  + f,6W  + f3«V  + af460x 
-afs8tj)y  - af 6 60z) ] Rd 02  = 0 

Observation  of  the  expression  leads  to  the  following  definition  of  the  one-dimensional 
strain  measures 


(16e) 

(16f) 


e = i(U'  - W) 

(18a) 

1 a' 

Ky  = -0y 

(18b) 

Kz  = + ^x> 

(18c) 

«x  = -^  <*x  - *z> 

(18d) 

7z  = i (U  + w'  R <py) 

R 

(18e) 

7y  =-^  (V'  - R0Z) 

(18f) 
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This  completes  the  derivation  of  the  one-dimensional  theory.  The  governing  equations ; 
in  terms  of  the  cross-section  displacements  and  rotations  can  be  obtained  directly  from 
equation  (15)  by  application  of  the  calculus  of  variations.  Here  the  alternate  procedure 
of  obtaining  the  equilibrium  equation  in  terms  of  the  stress  resultants  and  moments  will 
be  used.  The  governing  equations  in  terms  of  the  displacements  and  rotation  are  then 
obtained  by  use  of  the  stress-displacement  relations  (16). 

Governing  Differential  Equations 


The  equilibrium  equations  in  terms  of  the  stress  resultants  and  moments  are  obtained 
from  equation  (17)  by  integrating  by  parts  where  possible  to  yield 


i 


^ (-t'  + qz)  - naf2  ] SU  - (i  (t  + q')  + 7raf,]6W+ 


(19) 


. 1 , 7ra  o . 1 

^R-^y  + 9z  ~ p-N2J(U  + W')  + 7ra2fs  ] 50y  - + 7raf3]6V  + 

.1,1  , na  o 

l-(mz  - mx)  - qy  + - N22V'  + wa2f6]5</>z  - 

R R 

[■^(mz  + mx)  + 7ra2f4]5<£x|  Rd02  + 

a 

(t6U  - my6<*>y  - mz60z  + mx6<px  + qy8V  + qz5W]  |_aF0 


Since  the  variational  displacements  are  independent,  satisfaction  of  the  above  equality 
requires  that 


^ (-t'  + qz)  - rraf2  = 0 
n 

- — (t  + q^)  - waf,  = 0 
H 

my  + qz  — rr^-N22(U  + W' ) + Tra2 f s = 0 

-- 1 qy  - jrafj  = 0 

-R(m'z  - mx)  - qy  + n N22V'+ffa2f6  = 0 

— B ,mz  + mx>  - ’raJf4  = 0 
n 


(20a) 

(20b) 

(20c) 

(20d) 

(20e) 

(20f) 
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over  the  region  — a<02<a.  On  the  boundaries  02=— a and  02*a  one  of  each  of  the 
following  pairs  of  parameters  must  be  specified 


t or  U 
my  or  0 y 
qz  or  W 

(21) 

mz  or  0Z 
mx  or  0X 
qy  or  V 

To  complete  the  development,  we  substitute  the  stress-displacement  relations  (16)  into 
the  equilibrium  equations  to  give  the  following  governing  equations  for  the  displacements 
and  rotations 

-C 22-|j-  (2  + (-|)2Z2)U"  +-|5-  (C33Z3  + N2°2Z2)U 

— (— )3C22Z20y  + — (C33Z3  + (—  )2N22Z)0y  (22a) 

+ p7  [C22(2  + |-j^-)2Z2)  + C33Z3  + N22Z2]W'  — af2  = 0 

-1^-  (C33Z3  + N22Z2)W"  + C22±-  (2  + (1)2Z2)W 

tc2 2 (2  + (^t)Z2)+  C33z3  + N2JZ2]U'  (22b) 

C22Z2  + C33Z3  + (|)JN22Z]0y'  - af,  = 0 
R R R 

-(|)2aC22Z20;'  + (C33Z3  + (|)2N22Z)a0y 

-<!>3C22Z2U»  + -|<C33Z3  + (|)2N22Z)U  (22c) 

+ i[(I)2c22z2  + C3JZ3  + (1)2N22Z1W’  + a2fs  = 0 
R R R 
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kL 


-L_  (c33z2  + n|2z3)V"  + ( ->3c33z20"x 
R2  R 

+ -[C33Z2(1  - (-V)  + l|)2N22Z4]*'  - af3  = 0 
R R R 

— (j^2aC2 2Z302  + [C33(2(|)2  + (1  - {^)l)2Z2)  + (l)iN22Z4la0z 

-n|)2C22Z3  + C33(l)2(2  - Z2  + (l)2Z2)]atfx 
R HR 

-i[C33Z2(1  - (l)2)-  N2°2(1  - Z3)]V'  + a2f6  = 0 
R R 

~(i)2aC33(2  + (|)2Z2)0x"  + (|)2aC22Z3*x 

n R R 

+(|)2[C22Z3  + C33(2  - Z2  + (±)2Z2)]a «' 

R R 


(22d) 


(22e) 


(22f) 


+ {^)3C33Z2V"  - a2 f4=  0 
R 

This  gives  a set  of  six  linear  second  order  differential  equations  with  constant  coefficients 
in  the  six  cross-section  displacement  parameters.  The  first  three  of  these  equations 
describing  the  deformation  in  the  plane  of  the  arch  are  uncoupled  from  the  last  three 
which  describe  the  deformation  out  of  the  plane  of  the  arch.  This  completes  the 
development  of  the  theory  for  the  structural  behavior  of  pressure-stabilized  arches.  The 
theory  is  put  in  nondimensional  form  in  Appendix  C.  Although  the  motivation  for  this 
derivation  was  to  establish  equations  to  describe  or  predict  the  behavior  of 
pressure-stabilized  arches  under  applied  static  loads,  these  equations  also  represent  a theory 
for  the  elastic  stability  of  arches.  That  is,  if  the  applied  loads,  fj,  i=1,  6,  are  set  equal 
to  zero  and  homogenous  boundary  conditions  areQused,  we  have  a set  of  homogenous 
equations  for  the  determination  of  the  values  of  N22  giving  nontrivial  solution  functions 
for  the  displacement  parameters.  These  values  of  N22  are  the  buckling  stresses,  and  the 
solution  functions  are  buckling  mode  shapes. 

Discontinuity  Conditions  for  Green's  Function: 

One  additional  item  that  will  be  required  is  a derivation  of  the  discontinuity  conditions 
required  for  obtaining  a Green's  function  solution.  The  Green's  function  for  structural 
problems  can  be  interpreted  as  the  deformation  at  a general  point  02  due  to  a unit  load 
applied  at  02  The  discontinuity  conditions  involve  this  unit  load  and  internal  stress 
resultants  and  moments.  To  derive  these  conditions  we  consider  an  arch  having  an  angular 
span  of  2a  with  the  origin  of  the  coordinate  02  positioned  at  the  center  of  the  span 
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as  shown  in  Figure  3.  The  load  applied  at  02=#2  has  some  distribution  with  respect 
to  the  coordinate  0,  as  shown  in  Figure  3 and  requires  three  components,  G,  (0, ),  G2  (0, ), 
and  G3(0,),  for  complete  specification.  The  arch  is  divided  into  two  regions,  one  on 
either  side  of  the  applied  load  or  mathematically  -a<02<82  and  &2<d2<a.  To  obtain 
the  needed  conditions  the  energy  principle  (11)  is  rewritten  integrating  over  each  of  the 
regions  separately  and  including  the  potential  energy  of  the  forces  applied  at  02**$2 


.7, 


,2n 


S f f E(0,  02)(1  + a/Rcos(0,))Rad0,d02  + 
a * o 

a -2n 

8J  ] E(0i,M(1  + a/Rcos(0,))Rad0,d02  - 

7-, 


(23) 


f 2?r 

6 I (G,  (0, ) ui  (0,  52 ) + G2(0,)u2(0,02)  + G3  (0, ) u3(0l02)]ad01  = 0 

* n ' 9 • 


where  E(0|,02)  is  the  energy  density  function  and  the  displacement  components  in  the 
third  integral  are  evaluated  at  02=72.  The  simplified  displacement  approximations  are 
substituted  as  done  previously  and  the  integration  with  respect  to  0 , carried  out  to  give 


r C 

S I E(02)aRd02  + 5 I E(02)aRd02  - 

J*  Jtj2 

5 { ™lg2U$2)  + g,W$2)  + g3V(02)  + 


(24) 


g4a0x(02)  - gsa0y(^2)  - g6a0z(<5“2)^  = 0 


where 


9i 

= 1 
ir 

fill 

(G,  (0,  )sin(0, ) - G3(0,  )cos(0,  ))d0, 

(25a) 

= i 

it 

f nt 

9i 

J G2  (0, )d 0, 

(25b) 

, ± 

f 11T 

93 

]a  (G,  (0,  )cos (0,)  + G3(0|  )sin(0,  ))d0, 

(25c) 

= 1 

7T 

flit 

94 

Jo  g,(0,)<J0, 

(25d) 

= 1 
it 

ent 

9s 

I G2  (0,  )cos(0,  )d0, 

* O 

(25e) 

= j. 

it 

f 11T 

96 

] G2  (0,  )sin(0,  )d0, 

(25f) 
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Over  the  region  -a  < 02  <&2,the  energy  density  is  expressed  in  terms  of  the  displacement 
parameters  U|,  W,,  V,,  <t>xi , 0yi,  and  0Zl  and  over  the  region  £2  < 02  < a in  terms 

of  U2 , W2,  V2,  <t>XJi  and  1^  tj)e  variation  in  (24)  is  carried  out,  the  result 

expressed  in  terms  of  stress  resultants  and  moments,  and  the  integration  by  parts  performed, 
a set  of  six  governing  equations  identical  in  form  to  (22)  are  obtained  from  each  of  the 
two  regions  being  considered.  What  is  of  interest  here  are  the  boundary  terms  and  the 
applied  force  a\V2.  This  boundary  term  which  must  vanish  is 


t,6U,  - myi6</>y,  _mzi60Zl  + + dyi5Vi  + qZl6W, 


3, 


t26U2  - my 26«y2  - mZ26<*>22  + mX2^XJ  + qy28V2  + q 


,Z28W2 


-arr[g  28U,  (02)  + g,8W|(02)  + g36V,(02)  + g4aS 0X|(52) 


(26) 


— gs a50y , (02 ) - g6a80z|  (02)  )=  0 

# 

In  writing  (26)  subscript  1 refers  to  a variable  defined  on  the  region  -a  < 02  < 52 
and  subscript  2 to  a variable  defined  on  ?2  < 02  < a.  The  boundary  terms  at  02=±a 
give  the  boundary  conditions  at  the  ends  of  the  arch  which  are,  not  unexpectedly,  identical 
with  those  derived  previously.  The  remaining  boundary  terms  will  yield  the  conditions 
that  must  hold  at  the  applied  load,  but  it  is  first  necessary  to  state  the  conditions  of 
continuity  of  deformation  between  the  two  regions.  These  conditions  are  simply 

U,(02)  = U2(?2) 

VMfrj)  = W2(S2) 


V,(£2)  = V2(?2) 

0x,<^)  = #*(*,» 


#y,(^"2)  <t>y2  (5*2  ) 

0Z|  (^j)  * 0Z2(7 j) 
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Thus  the  variations  of  these  parameters  must  be  equal  and  the  boundary  terms  at  W2 
become 

[t,  - t2  - a7rg2]6U,(02)  + [-myi  + myj  + wa2 gs ] 60y,  (S2) 

fqzi  ~ ^Z2  ~ rrag,  ] 6W,  (02 ) + [~mzi  + ™Z2  + 7ra2 g6  J 50Zl  ) ^8) 

lmx,  " mx2  ' ^2g4]«0xl(^)  + [qyi  _ Qy2  _ jrag3]  6V,  (?2)  = 0 
Each  of  the  terms  must  vanish  independently  giving  the  discontinuity  conditions  as 

MM  ~ t2(02)  = irag2  (29a) 

myi(02)  - my2(02)  = 7raJgs  (29b) 

«*«(?« ) ~ qZ2(02)  = Trag,  (29c) 

mz,(52)  - %j(52)  = rra2g6  (29d) 

™xi(52)  = ffaJ94  (29e) 

qyi(52)  - qy2  (&2 ) = jrag3  (29f) 

These  are  the  required  conditions  and  they  are  expressed  in  nondimensional  form  in 
Appendix  C. 

Solution  of  the  Governing  Equations 

Having  derived  a set  of  governing  equations  in  the  previous  section  we  turn  now 
to  the  solution  of  these  equations.  In  what  follows  only  the  inplane  deformation  will 
be  considered  with  solutions  obtained  in  the  form  of  a Green's  Function  and  for  an  arch 
with  a uniform  normal  load.  The  Green's  Function  provides  a means  for  obtaining  solutions 
for  quite  general  loading  by  integration. 

Homogenous  Solution 

The  equations  to  be  solved  are  (22a),  (22b),  and  (22c)  which  are  given  in 
nondimensional  form  in  Appendix  C as  (C22),  (C23)  and  (C24).  These  are  rewritten  here 
in  more  compact  form  to  aid  in  writing  the  solutions. 


-H,lT  + H 2U  - H 3«"  + H40y  + (H,  + H2)|iV'  = 0 (30a) 

-HjW"  + H.JV  - (H3  + H4)^  - (H,  + Hj)i/  - 0 (30b) 

-pH30y  + pH40y  - H 3U"  + H 4U  + (H3  + H4)W'  = 0 (30c) 

where 

H,  = 2d/p  + dZj/p3  (31a) 

Hj  = cZ3/p  + nZ2/p  (31b) 

H3  = d Z2/p2  (31c) 

H4  = cZ3  + nZ/p1  (3  Id) 

Since  the  homogenous  solution  is  sought,  the  force  parameters  are  set  equal  to  zero.  The 
equations  have  constant  coefficients,  so  the  solution  is 

U = Aew02  (32a) 

«y  = Bew02  (32b) 

W = Cew0J  (32c) 


Substitution  of  this  solution  into  (30)  yields  the  following  set  of  homogenous  equations 
for  A,  B,  and  C: 

H2  - H,  cu2  H4  - H3to2  (H,  + H2)w 
H4  — H3oj2  p(H4  — H3oj2)  (H3  + H4)w 
— (H,  + H2)lo  — (H3  + H4)cj  H,  — H2co2 
Existence  of  a nontrivial  solution  requires  that  u be  a solution  of 

H,H4(H4  - pH2)  + cj2  [H2 H3 (pH,  - H3)  + 2H,H4(H4  - pH2)]  + 

w4IH,H4(H4  - pH2)  + 2H2  H3(pH|  - H3)J  + (34) 

w6[H2H3(pH,  - H3)l  = 0 

This  is  a sixth  order  equation  in  cj  or  a cubic  equation  in  u>2  which  can  be  written 
in  factored  form  as 

(w2  + 1)(oj2  + 1) [H2 H3 (pH,  - H3)o;2  + H,  H4(H4  - pH2)]  = 0 (35) 


25 


for  which  the  roots  are  (-i,  -i,  i,  i,  X,  -X)  where 


H,  H4  (pH?  - H4)  np2  (2  + Z,/p2)(cZ3  + nZ /p2) 

X ~ H2H3(pH,  - H3)  " 2Z2d(cZ3  + nZ2)  (36> 

Thus  X2  is  always  positive  and  the  homogenous  solution,  taking  account  of  the  repeated 
roots,  is 


U = A,cosh(X02)  + A2sinh(X02)  + A3sin<02)  + A4cos(02)  + 

A502sin(02)  + A602cos(02) 


(37a) 


= B,cosh(X02)  + B2sinh(X02)  + B3sin(02)  + B4cos(02)  + 

B502sin(02)  + B602cos(02) 


(37b) 


W = C,cosh(X02)  + C2sinh(X02)  + C3sin(02)  + C4cos(02)  + 

C502sin(02)  + C602cos(02) 


(37c) 


The  satisfaction  of  (27)  implies  a relationship  among  the  constants  Aj,  Bj,  and  C-,  such 
that  two  sets  can  be  expressed  in  terms  of  the  third.  Here  Aj  and  Bj  are  expressed 
in  terms  of  Cj  as 

A|  = 0i  C2 

A2  = 0i  C, 


A3  = P3C4  + psC5 
A4  = — 03  C3  + 05  C6 


(38) 


As  = 03  ^6 


Aft  ~ 03  ^5 


r 


i 


B I = 02  c 2 

Bj  = PjC, 

B3  = 04^4  + 06^5 

(39) 

B4  = —04^3  + 06  C6 

B5  * 04  C6 

B6  = —04  C5 

Expressions  for  the  determination  of  the  0-,  are  given  in  Appendix  D.  With  these  relations 
among  the  constants  the  solution  is 

U = C 1 0|  sinh(X0 2 ) + C20,cosh(X02)  - C303cos(02)  + C403sin(02)  + 


C5  (05 sin(02 ) - 0302cos(02))  + C6(05cos(02)  + 0302sin(02)) 


(40a) 


dy  = C,02sinh(X02)  + C202cosh(X02 ) - C3/J4cos(02)  + C404sin(02)  + 

C5(06sin(02)  - 0402cos(02 ))  + C6(06cos(02)  + 3402sin(02)) 


(40b) 


W =C,cosh(X02)  + C2sinh(X02)  + C3sin(02)  + C4cos(02) 


(40c) 


C502sin(02)  + C602cos(02) 

This  is  the  homogenous  solution  containing  six  unknown  constants  to  be  determined  by 
the  boundary  conditions  for  a given  set  of  applied  loads. 

Green's  Function 


Instead  of  attempting  to  find  particular  solutions  for  a variety  of  loading  situations, 
a Green's  function  is  obtained  which  will  provide  a solution  for  general  load  through 
integration.  As  stated  above,  the  Green's  Function  has  for  structural  problems  the 
interpretation  of  the  deflection  at  a general  point  02  due  to  a unit  load  applied  at  V2- 
As  in  the  derivation  of  the  discontinuity  conditions  we  consider  a beam  having  an  angular 
span  of  2a  with  the  origin  located  at  midspan,  as  shown  in  Figure  3.  Over  the  region 
-a  < 02  < a solution  Ut,  0yi,  and  W , is  defined  in  terms  of  integration  constants 
Cj,  i=  1,6  and  over  the  region  £2  < 02  < a a solution  U7.  and  W7  is  defined  in 
terms  of  integration  constants  Cj,  i = 7,12.  Thus  the  complete  solution  is  defined  by 
12  constants  and  12  conditions  must  be  specified  for  their  determination.  Six  of  these 
come  from  the  boundary  conditions,  three  on  each  end.  The  other  six  come  from  the 
three  conditions  of  continuity  of  displacements  and  the  three  discontinuity  conditions 
on  the  internal  stress  resultants  and  moments.  If  both  ends  are  simply  supported  these 
conditions  are 
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Ui  (-a)  = 0 


Myi  (—a)  = 0 
W,  (— ot)  = 0 
U2(a)  = 0 
Myj(a)  = 0 
W2  (a)  = 0 

Uy{&2)  - U2(V2)  = 0 


(41) 


tyyfc)  ~ 0y2<82)  = 0 
Wy(V2)  - MM?2)  = 0 
T,(ffa)  — T2(?2)  — 7rg2 

My,(^2)  - My,^)  = *9s 

Qz,(5’2)  - Q21(?2)  = -*g, 

and  if  both  ends  are  fixed,  the  conditions  Myi  (—  a)  = 0 and  My2  (a)  = 0 are  respectively 
replaced  by  0 y,  (— «)  = 0 and  <£yj  (a)  = 0.  Other  boundary  conditions  at  ±o  are  also  possible. 

These  conditions  give  a system  of  12  linear  equations  for  the  constants  Cj,  i = 1,12.  The 

expressions  for  the  coefficients  of  this  system  of  equations  are  given  in  Appendix  E. 

Because  only  inplane  deformation  of  the  arch  is  being  considered,  the  applied  forces 
G, , G2,  and  G3  must  be  symmetric  about  the  plane  of  the  arch.  Two  loading  situations 
are  of  particular  interest,  normal  load  and  vertical  load.  For  normal  loading,  the 
components  Gi  and  G2  vanish,  therefore  the  nondimensional  generalized  forces  §2  and 
Vs  vanish  and  the  load  is  described  in  terms  of  9, , and  for  a unit  load  g,  is  set  equal 
to  unity.  A vertical  load  is  defined  to  be  parallel  to  the  line  02  = 0.  For  a vertical 
load  distribution  G(0t)  located  at  #2,  the  force  components  are 

G,  = -G(0i  )sin(0!  )cos(32) 

G2  = — G(0,  )sin(?2)  (42) 

G3  = G(0)  )cos(0,  )cos((F2) 
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and  the  generalized  forces  are  found  by  carrying  out  the  integrals  in  the  definition  (25). 
For  example,  if  we  have  two  concentrated  forces  of  magnitude  G/2  located  at  0 , = ± ft, 
then 

G(0,)  = (G/2) [6(0,  + £,)  + 6(0,  - £,)]  (43) 

where  6 here  denotes  the  Dirac  delta  function.  Substituting  (43)  and  (42)  into  (25) 
and  carrying  out  the  integrals  using  the  integral  property  of  the  delta  function  yields 

7j,  = -gcos(82) 

g2  = -gsin(5’2) 

$3  = — gsin(?2)cos(£, ) 

for  a unit  load,  g = G/7rC22  is  set  equal  to  unity. 

Uniformly  Distributed  Normal  Load 

The  solution  for  a uniform  normal  load  could  be  obtained  by  integration  of  the 
Green's  function  derived  in  the  previous  section.  However,  since  a particular  solution 
of  the  nonhomogenous  equations  for  this  loading  condition  is  easily  found  it  was  decided 
to  give  the  solution  explicitly  as  the  sum  of  this  particular  solution  and  the  homogenous 
solution.  For  this  loading  condition  we  deal  with  the  forces  F, , F2,  and  F3fwithF,  and 
F2  vanishing.  F3  is  independent  of  02  and  has  a variation  in  0,  which  is  symmetrical 

about  0,  = 0 and  is  either  constant  or  zero  to  give  a distribution  similar  to  that  shown 

in  Figure  4.  Such  a distribution  will  cause  all  the  nondimensional  force  parameters  to 
vanish  except  F,  which  will  be  a constant.  Thus  we  need  a particular  solution  of  30a, 
30b,  and  30c  with  F,  appearing  as  a constant  nonhomogenous  term  in  (30b).  The 
required  particular  solution  is 

U = 0 

0y  = 0 (44) 

W = T,  /H, 

The  complete  solution  is  the  sum  of  this  particular  solution  and  the  symmetric  part  of 
the  homogenous  solution 

U = C,0,sinh(X02)  + C403sin(02)  + C5 (05sin(02)  - 0302cos(02))  (45a) 

0y  * C,02sinh(X02)  + C404sin(02)  + Cs (06sin(02 ) - 0402cos(02 ))  (45b) 

W = C,cosh(X02)  + C4cos(02)  + C502sin(02)  + pf,/d(2  + Z 2/p2)  (45c) 


30 


FIGURE  4.  ILLUSTRATION  OF  UNIFORM  NORMAL  LOAD  DISTRIBUTION 


Only  the  symmetric  part  of  the  solution  is  included  because  the  loading  is  symmetric 
about  02  = 0,  thus  the  deformation  will  be  symmetric.  Here  symmetric  deformation 
is  used  to  characterize  the  situation  where  the  displacement  W is  symmetric  and  U and 
<j>yi  are  antisymmetric,  thus  the  integration  constants  C2,  C3,  C6  are  9et  equal  to  zero. 
The  remaining  three  constants  C,,  C4,  Cs  are  determined  from  the  boundary  conditions 
which  for  simply  supported  end  conditions  are 

U(  a)  = 0 

My  (a)  = 0 (46) 

lY(a)  = 0 

and  for  fixed  end  conditions  My(a)  = 0 is  replaced  by  0y( a)  = 0.  In  both  of  these  cases 
symmetric  boundary  conditions  nave  been  assumed.  If  the  boundary  conditions  to  be 
applied  to  each  of  the  two  ends  are  different,  the  complete  homogenous  solution  must 
be  used  in  conjunction  with  the  particular  solution.  Equations  (46)  for  the  simply 
supported  arch  or  with  the  substitution  for  fixed  ends  gives  a system  of  three  equations 
for  Ci,  C4,  and  C5.  Expressions  for  the  coefficients  of  this  system  are  given  in 
Appendix  E.  A computer  program  to  carry  out  the  computations  associated  with  this 
solution  and  the  Green's  function  is  described  in  Appendix  F. 

EXPERIMENTAL  ANALYSIS 

To  provide  some  data  on  which  to  base  a judgment  of  the  validity  of  the  theory 
developed  above,  the  deformation  and  load-carrying  capability  of  a series  of 
pressure-stabilized  arches  were  determined  experimentally.  The  elastic  and  shear  moduli 
of  the  fabric  used  in  making  the  arches  were  also  measured  experimentally.  These  moduli 
are  needed  to  carry  out  the  computations  of  the  theoretical  predictions  for  comparison 
with  experimental  results. 

Arch  Loading  Tests 

The  objective  of  these  tests  was  to  measure  the  deformation  behavior  and  load-carrying 
capacity  of  pressure  stabilized  arches  of  several  sizes  and  inflation  pressure  levels.  The 
flexibility,  that  is,  the  deflection  per  unit  load,  is  taken  as  the  measure  of  the  deformation, 
and  the  wrinkling  load  is  used  to  measure  the  load-carrying  capability.  The  wrinkling 
load  is  defined  as  the  load  at  which  the  maximum  compressive  stress  due  to  applied  load 
is  equal  in  magnitude  to  the  tensile  stress  due  to  pressurization.  Thus  any  further  increase 
in  load  would  cause  wrinkling  of  the  fabric  because  it  cannot  support  compressive  stresses. 
The  arches  consist  of  a fabric  skin  which  gives  the  structural  strength  to  support  the 
pressure,  a bladder  to  prevent  leakage  of  the  inflation  gas  and  end  caps.  The  fabric  skin 
was  woven  with  beta  fiberglass  yarn  using  a three-dimensional  weaving  technique  that  results 
in  a fabric  having  the  natural  contour  of  a torus.  That  is,  when  the  arch  is  inflated 
it  has  the  shape  of  a torus  and  no  forces  are  required  to  maintain  it  in  that  shape.  This 
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weaving  process,  which  is  described  in  references  6 and  7,  produces  a fabric  in  which 
the  circumferential  or  fill  yarn  count  is  higher  at  the  inner  radius  of  the  arch  than  at 
the  outer  radius.  This  nonuniformity  of  yarn  count  and  the  resulting  natural  contour 
or  shape  causes  some  difficulty  in  the  determination  of  the  fabric  stiffness  properties  as 
will  be  discussed  subsequently.  The  fabric  used  was  a plain  weave  with  an  average  yarn 
count  in  both  warp  and  fill  of  512/m  using  133  Tex  fiberglass  yarns.  Arches  having 
cross-section  radii  of  0.038,  0.051,  0.076,  and  0.089  m,  all  having  an  inner  arch  radius 
of  0.914  m,  were  woven  and  tested.  The  arch  radius  to  cross-section  radius  ratios  for 
these  arches  are  then  respectively  25,  19,  13,  and  11.3. 

A polyethylene  film  bladder  was  used  to  provide  retention  of  the  inflation  gas.  These 
bladders,  made  from  4-mil-thick  film,  were  shaped  to  conform  to  the  arch  shape  by  cutting 
two  segments  of  a circular  annulus  and  joining  them  together  by  heat-sealing  at  the  inner 
and  outer  radii.  These  bladders  were  made  with  the  cross-section  about  10%  larger  than 
that  of  the  corresponding  fabric.  In  doing  this  the  strength  requirements  of  the  bladder 
were  not  a consideration.  However,  because  of  the  oversizing  of  the  bladders  it  was  possible 
for  folds  to  form  in  them,  and  when  pressurized,  the  folds  were  flattened  out  resulting 
in  cracks  that  caused  leakage.  This  is  believed  to  be  a baldder  material  problem  that 
could  be  cured  by  choosing  a material  with  greater  crack  resistance. 

The  ends  of  the  arches  were  closed  and  restrained  with  end  cap  assemblies  shown 
in  Figure  5.  A schematic  of  this  assembly  is  shown  in  Figure  6.  As  is  shown,  the  fiberglass 
fabric  and  polyethylene  bladder  go  through  the  inside  of  and  wrap  around  the  sealing 
ring.  This  arch-sealing  ring  combination  is  seated  in  the  end  cap  on  the  rubber  gasket. 
The  arch  and  sealing  ring  are  captured  by  the  retaining  ring  which  is  bolted  to  the  end 
cap  to  hold  the  assembly  together  and  apply  pressure  to  seal  the  arch  for  inflation.  Inflation 
is  accomplished  through  the  port  shown  in  Figure  6.  As  is  seen  in  Figure  5,  the  end 
cap  assembly  was  bolted  to  a steel  plate  so  that  it  could  not  rotate.  Thus  the  end  restraint 
simulated  fixed  or  clamped  ends.  Each  arch  had  two  end  caps  placed  so  that  the  arches 
had  an  angular  span  of  n radians. 


6.  Koppelman,  Edward  and  Arthur  R.  Campman;  Woven  Fabrics  and  Method  of  Weaving; 
US  Patent  No.  2998030,  Aug  29,  1961. 

7.  Koppelman,  Edward  and  Arthur  R.  Campman;  Method  of  and  Apparatus  for  Weaving 
Shaped  Fabrics  and  Articles  Woven  Thereby;  US  Patent  No.  313671;  May  12,  1964. 
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Figure  5.  End  Closure  and  Restraint  for  Fiberglass  Arches. 
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Figure  6.  Schematic  of  the  End  Cap  Assembly. 


35 


I 


A schematic  of  the  test  apparatus  is  shown  in  Figure  7 and  a photograph  of  the 
actual  test  apparatus  is  shown  in  Figure  8.  The  loading  for  all  tests  was  a concentrated 
load  at  midspan.  This  load  was  applied  using  a 2.5-cm-wide  strap  wrapped  around  the 
arch.  Tests  were  conducted  with  the  load  directed  toward  the  arch  center  of  curvature 
and  away  from  the  center  of  curvature.  Wrinkling  was,  as  expected,  observed  only  when 
the  load  was  directed  toward  the  center  of  curvature,  as  this  loading  condition  is  the 
one  which  puts  the  arch  in  a compressive  state  of  stress.  The  magnitude  of  the  applied 
load  was  measured  with  a strain  gage  force  transducer.  The  other  measurement  made 
was  of  the  arch  deflection  at  the  point  of  the  applied  load.  This  was  accomplished  with 
a linear  variable  differential  transformer  positioned  to  measure  the  normal  component  of 
the  deflection.  The  output  from  these  transducers  was  amplified  when  required  and 
recorded  on  a strip  chart  recorder,  a sample  of  which  is  shown  in  Figure  9,  for  both 
the  flexibility  and  the  wrinkling  load  tests.  Three  tests  were  carried  out  for  each  pressure 
level  with  good  repeatability. 

To  facilitate  comparison  of  these  experimental  results  with  the  theoretical  prediction, 
the  measured  force  and  displacement  were  converted  to  the  nondimensional  parameters 
used  in  the  theory.  The  measured  normal  component  of  displacement,  W,  is  put  in 
nondimensional  form  by  the  relation 


W = W/a 


(47) 


where  W is  the  nondimensional  displacement  and  a is  the  cross-section  radius.  Conversion 
of  the  measured  force  Fg  to  the  nondimensional  force  parameters Tfj,  i = 1,6  is  illustrated 
in  Figure  10.  As  shown  in  this  figure,  G3  is  the  only  nonvanishing  component  of  the 
line  load  applied  to  the  arch  and  this  is  in  the  negative  direction.  Thus  the  nondimensional 
force  parameters  of  concern  are  "5,  = g,/C22  and  ^ = g2/C22.  Using  equations  (25a) 
and  (25c)  we  have 
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(1/ttC 


22 


r 
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4 

+ n 


G3(0,  )cos(0]  )d0, 


■g3  =-(1/7re22)T  GjtfJsinlMdO, 


(48) 


(49) 


where  C2  2 is  a reference  value  of  the  elastic  modulus  and  the  negative  direction  of  the  load 
G3(0, ) has  been  accounted  for.  Since  G3  (0  2 ) is  nonvanishing  only  over  a portion  of  the 
periphery  of  the  arch,  the  limits  and  integrals  can  be  changed  to  give 
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FIGURE  7.  SCHEMATIC  OF  TEST  APPARATUS 


Figure  8.  Arch  Model  in  Test  Apparatus. 
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FLEXIBILITY  MEASUREMENT 
INWARD  LOAD 

\ DISPLACEMENT 
2.54  X I0-3  m/DIV. 


\ FORCE 

4.45  N/DIV. 


TIME 


diameter:  0.102m 

PRESSURE!  69  kPa 


Figure  9.  Typical  Recordings  of  the  Force  and  Displacement  From  Arch  Loading  Tests. 
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G3(0,  )sin{0,  Id©, 


(51) 


Because  of  the  symmetry  of  G3(0,)  about  0 , = 0,  §3  vanishes  and  we  need  only  be 
concerned  with  If, . To  obtain  a relationship  between  and  Fg  we  observe  that  G3(0, ) 
also  acts  on  the  webbing  used  to  apply  the  load  and  for  equilibrium  must  have 


piir-u 

Fa  =f  G3 (0,  )cos(0|  )ad0,  (52) 

* CJ 

Substitution  of  (52)  into  (50)  gives 

= F g/ira^  2 

the  expression  used  to  convert  the  measured  applied  load  to  the  nondimensional  force 
parameter.  If  the  subscript  w denotes  wrinkling  then  the  nondimensional  wrinkling  load 
is  given  in  terms  of  the  physical  wrinkling  load  as 

w - F aw^3^2  2 (53) 

the  magnitude  of  Fgw  is  determined  experimentally  by  finding  the  value  of  Fg  for  which 
the  loading  curve  becomes  nonlinear.  As  can  be  seen  from  Figure  9,  the  transition  from 
linear  behavior  to  nonlinear  behavior  is  smooth,  making  the  evaluation  of  Fgw  rather 
difficult.  This  difficulty  is  especially  pronounced  in  comparison  to  the  behavior  of 
pressure-stabilized  beams  which,  as  reported  in  reference  4,  exhibit  a constant  load  jump 
in  deflection  as  wrinkling  begins.  With  such  behavior  it  was  quite  easy  to  evaluate  the 
wrinkling  load  from  the  loading  curve. 

The  flexibility  is  the  transverse  displacement  due  to  a unit  load  which  in 
nondjrriensional  form  is: 

y = IV/g,  ' (54) 

Substituting  the  expressions  for  IV  and  U,  gives  the  following  expression  for  the 
nondimensional  flexibility  in  terms  of  the  measured  force  and  displacement: 

7 = 7tS2  2 W/Fa  (55) 
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Determination  of  Material  Stiffness  Properties 
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In  order  to  accomplish  the  objective  of  comparing  the  theory  with  experimental  results 
obtained  by  the  above  described  procedure,  the  elastic  and  shear  moduli  of  the  fabric 
used  in  the  arches  must  be  measured.  These  moduli  are  needed  for  the  theoretical 
computation  of  the  behavior  of  these  fiberglass  arches.  These  moduli  or  stiffnesses  were 
determined  from  the  results  of  tension  and  torsion  tests.  Although  these  tests  are  very 
similar  to  those  commonly  used,  they  were  performed  on  specimens  in  their  pressurized 
state  of  stress.  Testing  in  this  manner  yields  stiffnesses  relative  to  a state  of  stress  very 
close  to  that  present  in  the  arch  loading  tests,  that  is,  uniaxial  tension  or  shear  superimposed 
on  biaxial  tension  resulting  from  internal  pressure.  The  moduli  were  measured  using 
relatively  small  excursions  from  the  pressurized  state  of  stress,  and  the  behavior  was 
assumed  to  be  linear  over  these  excursions.  Thus,  for  a given  value  of  pressure,  the  elastic 
and  shear  moduli  are  constants.  However,  because  of  the  nonlinear  behavior  of  fabrics, 
these  moduli  were  found  to  be  dependent  on  the  pressure  level.  This  testing  procedure 
had  the  additional  advantage  of  yielding  the  moduli  of  the  fabric  skin-polyethylene  film 
composite,  and  although  the  bladder  is  thought  to  make  little  contribution  to  the  stiffness, 
any  contribution  is  included.  It  should  also  be  recalled  that  the  arches  were  woven  with 
a natural  curved  contour,  and  it  was  not  possible  to  devise  a simple  test  using  curved 
specimens  that  would  be  suitable  for  measuring  the  elastic  and  shear  moduli.  We  therefore 
obtained  some  straight  woven  fabric  tubes  having  a cross-section  radius  of  0.038  m for 
use  in  determination  of  the  moduli  in  tension  and  torsion  tests.  These  straight  tubes 
were  woven  with  the  same  fiberglass  yarns  and  plain  weave  design  used  in  the  contoured 
arches  so  the  stiffness  characteristic  of  the  straight  tubes  should  closely  represent  that 
of  the  arches. 

Elastic  Modulus 

As  indicated  above,  the  elastic  modulus  was  determined  from  the  results  of  a tension 
test,  as  is  shown  schematically  in  Figure  11.  After  pressurizing  the  specimen  the  test 
is  performed  in  the  usual  manner,  applying  an  axial  force  F and  measuring  its  magnitude 
along  with  the  corresponding  elongation,  e.  The  theory  requires  the  constant  C2  2 relating 
the  stress  resultant,  N22,  and  the  strain  e 22  as 

Nj2  - C22e22 

Thus,  C22  is  the  slope  of  the  stress  resultant-strain  curve  and  can  be  computed  from 
the  force-elongation  curve  generated  from  the  tension  test  by  the  expression 

C22  = lgAF/27raAe 

where  AF  and  Ae  are  respectively  increments  of  applied  force  and  specimen  elongation. 
In  the  tests  conducted,  the  specimen  gage  length,  lg,  and  the  radius,  a,  were  0.343  m and 
0.038  m,  respectively.  Typical  force-elongation  plots  from  these  tests  are  shown  in 
Figure  12. 
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Figure  11.  Schematic  Diagram  of  the  Tension  Test 


43 


FIBERGLASS  FABRIC  WITH  4 MIL 
POLYETHYLENE  BLADDER 

GAUGE  LENGTH  * 0.343m 
RADIUS  = 0.038m 
PRESSURE  = 138,000  Po 

Figure  12.  Typical  Force-Elongation  Recordings  From  Tension  Test. 
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Shear  Modulus 

The  shear  modulus  was  determined  from  the  results  of  a torsion  test,  as  shown 
schematically  in  Figure  13.  This  test  is  performed  by  inflating  the  specimen  to  a given 
pressure  and  applying  a torque,  T,  about  the  axis  of  the  cylinder  and  measuring  0,  the 
rotation  of  the  end  of  the  cylinder.  The  theory  requires  the  constant,  C33,  relating  the 
shear  stress  resultant,  N12,  to  the  shear  strain,  e12,  as 

^1 2 = C33C| 2 

The  parameter  C33  is  the  slope  of  the  stress  resultant-strain  curve  and  can  be  computed 
from  the  torque-rotation  curve  generated  from  a torsion  test  by  the  expression 

C33  = lgAT/27ra3A0 

where  AT  and  A0  are,  respectively,  increments  of  applied  torque  and  specimen  rotation. 
The  parameters,  l„  and,  a are  as  defined  above  and  for  the  tests  conducted  have  the  values 
0.267  m and  0.038  m,  respectively.  Typical  torque  and  rotation  plots  from  these  tests 
are  shown  in  Figure  14. 


RESULTS  AND  DISCUSSION 

In  this  section  we  will  present  the  experimental  results  obtained  with  the  procedures 
just  described  and  compare  the  arch  performance  characteristics  obtained  experimentally 
with  those  predicted  by  the  theory  developed  in  this  report. 

Material  Properties 

The  material  properties  required  by  the  theory  are  the  elastic  modulus  and  the  shear 
modulus.  As  is  indicated  above,  these  parameters  are  found  as  a function  of  pressure, 
as  shown  in  Figures  15  and  16.  In  these  figures  are  plotted  the  nondimensional  elastic 
and  shear  moduli  as  a function  of  the  nondimensional  pressure  parameter.  These  are 
the  nondimensional  parameters  used  in  the  theory  with  the  value  of  the  elastic  modulus, 
C22,  for  zero  pressure  being  used  as  the  reference  value.  It  is  denoted  as  C22.  The 
pressure  parameter  which  is  used  throughout  this  report  is  the  nondimensional  axial  stress 
resultant  due  to  internal  pressure  and  is  given  in  terms  of  the  pressure,  P,  and  the 
cross-section  radius,  a,  as: 


n = Pa/2C2 2 

where  n denotes  the  pressure  parameter.  In  both  Figures  15  and  16  the  experimental 
results  which  are  tabulated  in  Tables  1 and  2 are  denoted  by  the  symbols,  and  the  least 
squares  linear  fits  of  the  data  used  for  the  theoretical  prediction  of  the  arch  behavior 
are  shown  both  graphically  and  mathematically.  The  reference  value  of  the  elastic  modulus, 
C22,  is  also  given. 
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Figure  14.  Typical  Torque  and  Rotation  Recordings  From  Torsion  Test. 
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FIGURE  15.  ELASTIC  MODULUS  AS  A FUNCTION  OF  THE  PRESSURE  PARAMETER 
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TABLE  1 

ELASTIC  STIFFNESS  DATA  FOR  FIBERGLASS  FABRIC 


Pressure  Elastic  Stiffness 


Force 

Elongation 

Modulus 

P,  kPa 

n,  Nondim 

AF,  N 

Ae,  m 

Cj  2 = N/m 

34.5 

0.012 

49. 

1.0  x 10'3 

70400 

58. 

0.8 

104100 

57. 

1.0 

81900 

63. 

0.8 

113100 

59. 

1.0 

84800 

79. 

1.0 

113500 

61. 

1.0 

87600 

83. 

1.0 

119200 

115. 

1.0 

165200 

162. 

1.0 

232700 

69.0 

0.024 

81. 

1.0 

116400 

122. 

1.0 

175300 

94. 

0.9 

150000 

132. 

1.0 

189600 

82. 

0.7 

168300 

138. 

1.0 

198200 

80. 

0.7 

164100 

124. 

1.0 

178100 

81. 

0.7 

166200 

129. 

1.0 

185300 

81. 

0.7 

166200 

128. 

1.0 

183900 

173. 

1.0 

248500 

182. 

1.0 

261400 

Stress 


N22,N/m 

864 

1200 

780 

1280 

905 

1280 

864 

1280 

1910 

2830 

1520 

2270 

1600 

2230 

1560 

2190 

1520 

2100 

1520 

2100 

1520 

2100 

2980 

4200 
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TABLE  1 


ELASTIC  STIFFNESS  DATA  FOR  FIBERGLASS  FABRIC  (cont'd) 


Pressure 

Elastic  Stiffness 

Stress 

Force 

Elongation 

Modulus 

P,  kPa 

n,  Nondim 

AF,  N 

Ae,  m 

C2  2 # N/m 

n22,n  / 

103.5 

0.036 

64. 

0.5  x 10  3 

183900 

2300 

156. 

1.0 

224100 

3300 

76. 

0.5 

218400 

2220 

180. 

1.0 

258600 

3390 

80. 

0.5 

229800 

2220 

186. 

1.0 

267200 

3390 

78. 

0.5 

224100 

2260 

180. 

1.0 

258600 

3140 

214. 

1.0 

307400 

4390 

138.0 

0.048 

80. 

0.5 

229800 

2870 

200. 

1.0 

287300 

3710 

200. 

0.9 

319200 

4500 

98. 

0.5 

281600 

2870 

208. 

1.0 

298800 

3620 

243. 

1.0 

349100 

4630 

102. 

0.5 

293100 

2870 

218. 

1.0 

313200 

3710 

224. 

0.9 

357600 

4550 

172.5 

0.060 

120. 

0.5 

344800 

3530 

232. 

1.0 

333300 

4530 

218. 

1.0 

313200 

5370 

104. 

0.5 

298800 

3380 

244. 

1.0 

350500 

4450 

242. 

1.0 

347600 

5450 

88. 

0.4 

316000 

3530 

252. 

1.0 

362000 

4360 

248. 

1.0 

356300 

5450 
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TABLE  2 


SHEAR  STIFFNESS  DATA  FOR  FIBERGLASS  FABRIC 


Pressure 

Shear  Stiffness 

Torque 

Rotation 

Modulus 

P,  kPa 

n,  Nondim. 

AT,  N— m 

A 0,  rad. 

C3  3 , N/m 

34.5 

0.012 

2.32 

0.319 

5584. 

2.17 

0.305 

5510. 

2.32 

0.346 

5192. 

69.0 

0.024 

2.82 

0.346 

6312. 

3.15 

0.388 

6287. 

2.71 

0.332 

6312. 

103.5 

0.036 

3.11 

0.332 

7254. 

3.08 

0.332 

7184. 

3.11 

0.332 

7254. 

138.0 

0.048 

3.13 

0.319 

7599. 

3.22 

0.332 

7511. 

3.58 

0.360 

7701. 

172.5 

0.060 

3.94 

0.332 

9190. 

3.76 

0.319 

9128. 

3.66 

0.305 

9293. 

207.0 

0.072 

3.93 

0.305 

9979. 

3.93 

0.305 

9979. 

3.93 

0.305 

9979. 

In  reducing  the  data  for  the  elastic  modulus,  increments  of  axial  force  and  elongation 
were  taken  from  two  or  three  segments  of  the  plot  as  indicated  in  Figure  12.  Once 
the  modulus  was  computed  using  these  increments  a decision  had  to  be  made  regarding 
the  value  of  the  pressure  parameter  or,  equivalently,  axial  stress  associated  with  this 
modulus.  In  previous  work,  reference  4,  the  pressure  parameter,  or  axial  stress  due  to 
pressure  at  which  the  test  was  run,  was  used.  In  those  tests,  however,  only  a short  segment 
of  the  load-elongation  curve  beyond  the  pressurized  state  was  used.  In  the  present  work, 
as  indicated  above,  several  segments  beyond  the  pressurized  state  were  used,  and  it  was 
decided  to  use  a corrected  pressure  parameter  equal  to  the  sum  of  the  axial  stress  due 
to  pressure  and  the  axial  stress  due  to  axial  load  at  the  beginning  of  the  increment  used 
to  calculate  the  modulus.  This  corrected  value  of  the  pressure  parameter  is  given  in  the 
last  column  of  Table  1 which  is  headed  "Stress."  The  data  from  Table  1 is  plotted 
on  Figure  15  using  this  corrected  pressure  parameter. 

The  data  in  Figure  15  appears  to  have  some  curvature  indicating  that  a higher  order 
curve  might  provide  a better  fit.  However,  it  was  decided  that  in  the  present  context 
a linear  approximation  would  be  sufficiently  accurate. 

It  is  interesting  to  compare  the  values  of  elastic  stiffness  used  in  the  analysis  with 
the  limiting  value  obtained  by  considering  the  yarns  as  axial  springs  in  parallel.  The  axial 
stiffness  of  a fiberglass  yarn  is  given  as  28.2  N/Tex.  Thus  for  512  yarns  of  133  Tex 
in  parallel,  a stiffness  of  1.92  x 106  N/m  is  obtained.  The  maximum  value  of  stiffness 
used  in  the  analysis  is  0.62  x 106  N/m,  for  the  0.089-m  cross-section  radius  at  207  kPa 
pressure.  Thus  the  stiffnesses  were  all  less  than  one-third  of  the  limiting  stiffness.  The 
fiberglass  was  chosen  for  this  work  because  it  has  a linear  stress-strain  behavior,  and  it 
was  believed  that  this  would  simplify  the  analysis  problems.  As  can  be  seen  from  the 
data  shown  here,  the  resulting  fabric  is  far  from  linear.  The  stress-strain  behavior  exhibits 
a stiffening  effect  as  stress  is  increased,  which  is  typical  of  fabrics.  Over  the  stress  range 
shown  in  Figure  15  the  stiffness  increases  by  a factor  of  nearly  seven.  This  is  attributed 
to  the  yarn  and  fabric  structure  because  fiberglass  does  have  a linear  stress-strain  behavior 
in  comparison  with  other  materials  used  in  fabrics  such  as  nylon  and  dacron.  This  increase 
in  stiffness  of  the  fiberglass  fabric  is  a marked  contrast  to  very  small  increases  for  a dacron 
fabric  and  the  absence  of  increases  in  a nylon  fabric  reported  in  reference  4.  These  facts 
illustrate  some  of  the  difficulties  in  testing  fabrics  as  continuum  structures. 

The  results  from  the  torsion  test  are  presented  in  Figure  16  in  terms  of  a plot  of 
the  shear  modulus  as  a function  of  the  nondimensional  pressure  which  is,  as  noted  above, 
the  nondimensional  axial  stress.  These  results  are  rather  unremarkable,  showing  a shear 
stiffness  of  10  to  18  percent  of  the  reference  elastic  stiffness  or  modulus.  These 
percentages  are  somewhat  larger  than  those  of  the  nylon  and  dacron  fabrics  of  reference  4. 
Those  fabrics  had  shear  stiffness  of  3 to  6 percent  and  4 to  12  percent,  respectively, 
of  their  reference  elastic  stiffness.  However,  the  elastic  stiffness  of  the  nylon  and  dacron 
fabrics  were  found  to  be  approximately  three  times  the  elastic  stiffness  of  the  present 
fiberglass  fabric,  so  on  an  absolute  basis  the  three  fabrics  have  nearly  equal  shear  stiffness. 
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Behavior  of  Arche*  Under  Loed 


The  results  from  the  arch  loading  tests  along  with  companion  theoretical  predictions 
are  given  in  graphical  form  on  Figures  17  through  24  for  the  four  arches  described  earlier. 
The  experimental  data  in  numerical  form  on  which  the  plots  are  based  is  given  in  Tables 
3 through  6.  For  each  arch  a plot  for  both  the  flexibility  and  the  wrinkling  load  as 
a function  of  pressure  are  included  with  the  experimental  results  denoted  by  symbols 
and  the  theoretical  predictions  by  the  curves.  The  flexibility  is  defined  as  the  deformation 
per  unit  applied  force.  The  wrinkling  load  is  the  load  for  which  the  maximum  compressive 
stress  resulting  from  the  applied  load  is  equal  in  magnitude  to  the  tensile  stress  resulting 
from  pressurization.  Since  the  fabric  cannot  resist  compressive  stress,  any  further  increase 
in  the  applied  load  would  cause  wrinkling  of  the  fabric.  The  flexibility  and  the  wrinkling 
load  are  then,  respectively,  measures  of  the  deformation  behavior  and  the  load-carrying 
capability. 

Flexibility 

The  results  for  the  arch  flexibility  are  presented  in  Figures  17,  19,  21  and  23. 
Experimental  results  are  given  for  the  flexibility  with  the  applied  force  directed  inward, 
that  is,  toward  the  center  of  curvature  of  the  arch  and  outward,  away  from  the  center 
of  curvature.  One  would  not  expect  these  two  cases  to  be  significantly  different,  but 
for  the  arch  with  a radius  ratio  of  25.0,  the  flexibility  for  the  inward  force  is  significantly 
larger  than  that  for  the  outward  force.  A similar  situation  exists  for  the  arch  having 
a radius  ratio  of  11.3  except  the  differences  are  not  as  great  in  this  case.  For  the  other 
two  arches,  the  difference  is  of  the  order  of  magnitude  of  the  experimental  error.  In 
general,  the  difference  in  the  flexibility  resulting  from  the  change  in  loading  direction 
is  greatest  at  low  pressures  and  becomes  significantly  smaller  with  increasing  pressure.  It 
will  be  observed  that  the  flexibilities  associated  with  the  inward  force  are  greater  than 
those  associated  with  the  outward  force  except  for  the  arch  with  a radius  ratio  of  13.0 
and  for  this  arch  the  opposite  is  true,  but  the  differences  are  also  the  smallest  for  this 
case.  A possible  explanation  of  this  behavior  is  a differenital  stiffness  effect  caused  by 
the  stress  or  internal  force  due  to  the  applied  load.  In  the  arch,  in  contrast  to  the  straight 
beam,  lateral  deformation  is  accompanied  by  an  internal  axial  force  which  is  positive  for 
the  outward  force  and  negative  for  the  inward  force.  This  internal  force  is  analogous 
to  the  internal  force  resulting  from  the  pressure  and  if  it  is  large  enough  can  have  the 
same  influence.  The  pressure  internal  force  is  positive  and  decreases  the  flexibility  of 
the  arch.  Therefore,  if  the  internal  force  due  to  applied  load  were  to  have  this  same 
effect,  it  would  be  expected  that  under  the  outward  force  the  arch  would  be  less  flexible 
than  under  the  inward  force,  and  this  is  what  is  observed.  This  phenomenon  is  not  included 
in  the  theory  since  only  the  internal  axial  force  due  to  pressure  has  an  effect  on  the 
stiffness  or  flexibility  of  the  arch.  It  should  also  be  pointed  out  that  this  phenomenon 
is  of  no  aid  in  explaining  the  disagreement  between  theory  and  experiment. 
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TABLE  3 

TEST  DATA  FOR  ARCH  WITH  p = 11.3 
Flexibility 

Pressure  Wrinkling  Load 

Force  Outward  Force  Inward 


P,  kPa 

n,  Nondim. 

w/Fa,  m/N 

y,  Nondim. 

w/Fa,  m/N 

y,  Nondim. 

Faw,  N 

gw,  Nondim. 

34.5 

0.028 

6.79  x 10  s 

11.7 

9.59  x 1Q'S 

16.5 

365. 

2.38  x 10’2 

6.79 

11.7 

9.59 

16.5 

365. 

2.38 

6.62 

11.4 

9.59 

16.5 

365. 

2.38 

69.0 

0.056 

5.31 

9.15 

6.74 

11.6 

587. 

3.83 

CT> 

5.08 

8.75 

6.79 

11.7 

587. 

3.83 

CO 

4.91 

8.46 

6.74 

11.6 

587. 

3.83 

103.5 

0.084 

4.28 

7.38 

5.31 

9.15 

756. 

4.93 

4.22 

7.27 

5.08 

8.76 

756. 

4.93 

4.11 

7.08 

5.02 

8.65 

756. 

4.93 

138.0 

0.112 

3.71 

6.39 

4.51 

7.77 

1023. 

6.67 

3.65 

6.29 

4.34 

7.48 

1023. 

6.67 

3.60 

6.20 

4.28 

7.38 

1023. 

6.67 

172.5 

0.140 

3.25 

5.60 

,3.88 

6.69 

1179. 

7.69 

3.20 

5.52 

3.71 

6.39 

1179. 

7.69 

3.20 

5.52 

3.71 

6.39 

1179. 

7.69 

207.0 

0.168 

2.85 

4.91 

3.43 

5.91 

2.85 

4.91 

3.43 

5.91 

2.85 

4.91 

3.48 

5.91 

TABLE  4 


TEST  DATA  FOR  ARCH  WITH  p = 13.0 


Pressure 


Flexibility 


Force  Outward 


Force  Inward 


Wrinkling  Load 


P.  kPa 

n,  Nondim. 

w/Fa,  m/N 

7,  Nondim. 

w/Fa,  m/N 

7,  Nondim. 

Faw,  N 

gw,  Nondim. 

34.5 

0.024 

10.9  x 10  5 

18.8 

11.4  x 10  s 

19.6 

231. 

1.76  x 10  2 

10.9 

18.8 

10.9 

18.8 

231. 

1.76 

10.9 

18.8 

11.2 

19.3 

69.0 

0.048 

8.62 

14.8 

7.94 

13.7 

356. 

2.72 

8.39 

14.5 

7.82 

13.5 

356. 

2.72 

8.11 

14.0 

7.82 

13.5 

103.5 

0.072 

6.62 

11.4 

6.22 

10.7 

476. 

3.63 

6.91 

11.9 

6.05 

10.4 

467. 

3.56 

6.91 

11.9 

6.05 

10.4 

138.0 

0.096 

5.94 

10.2 

5.31 

9.15 

600. 

4.58 

5.82 

10.0 

5.08 

8.75 

600. 

4.58 

5.71 

9.84 

5.19 

8.94 

172.5 

0.120 

5.25 

9.05 

4.39 

7.57 

756. 

5.77 

5.08 

8.75 

4.28 

7.38 

756. 

5.77 

4.91 

8.46 

4.28 

7.38 

207.0 

0.144 

4.57 

7.88 

3.88 

6.69 

934. 

7.13 

4.45 

7.67 

3.83 

6.60 

912. 

6.96 

4.44 

7.65 

3.83 

6.60 

912. 

6.96 

r 


TABLE  5 

TEST  DATA  FOR  ARCH  WITH  p = 19 
Flexibility 

Pressure  Wrinkling  Load 

Force  Outward  Force  Inward 


P,  kPa 

n,  Nondim. 

w/Fa,  m/N 

y,  Nondim. 

w/Fa,  m/N 

7,  m/N 

Faw,  N 

gw,  Nondim. 

34.5 

0.016 

26.9  x 10  s 

46.4 

38.6  x 10  s 

66.5 

62.3 

0.71  x 10' 2 

27.2 

46.9 

38.8 

66.9 

62.3 

0.71 

26.2 

45.1 

38.3 

66.0 

62.3 

0.71 

69.0 

0.032 

20.7 

35.7 

22.9 

39.5 

116. 

1.32 

20.1 

34.6 

22.0 

37.9 

116. 

1.32 

20.5 

35.3 

22.0 

37.9 

120. 

1.36 

103.5 

0.048 

16.1 

27.7 

18.1 

31.2 

173. 

1.97 

15.9 

27.4 

17.5 

30.2 

173. 

1.97 

15.7 

27.1 

17.3 

29.8 

173. 

1.97 

138.0 

0.064 

13.8 

23.8 

14.8 

25.5 

218. 

2.48 

13.6 

23.4 

14.6 

25.2 

213. 

2.42 

13.3 

22.9 

14.4 

24.8 

213. 

2.42 

172.5 

0.080 

12.2 

21.0 

12.4 

21.4 

253. 

2.88 

11.9 

20.5 

12.2 

21.0 

258. 

2.93 

11.8 

20.3 

12.1 

20.8 

253. 

2.88 

207.0 

0.096 

11.1 

19.1 

10.4 

17.9 

325. 

3.70 

10.8 

18.6 

10.2 

17.6 

325. 

3.70 

10.7 

18.4 

9.99 

17.2 

325. 

3.70 

TABLE  6 


TEST  DATA  FOR  ARCH  WITH  p = 25 


Pressure 


Flexibility 


Force  Outward 


Force  Inward 


P,  kPa 

n,  Nondim. 

w/Fa,  m/N 

7,  Nondim. 

w/Fa,  m/N 

7,  Nondim 

34.5 

0.012 

40.8  x 10  s 

70.3 

81.6  x 10  s 

140.6 

39.1 

67.4 

81.6 

140.6 

39.1 

67.4 

81.6 

140.6 

69.0 

0.024 

32.8 

56.5 

45.7 

78.8 

30.7 

52.9 

46.4 

80.0 

29.4 

50.7 

46.4 

80.0 

103.5 

0.036 

24.2 

41.7 

37.3 

64.3 

22.8 

39.3 

35.5 

61.2 

22.4 

38.6 

35.9 

61.9 

138.0 

0.048 

19.3 

33.3 

29.0 

50.0 

18.7 

32.2 

27.9 

48.1 

18.7 

32.2 

27.7 

47.7 

172.5 

0.060 

16.3 

28.1 

24.8 

42.7 

15.6 

26.9 

23.6 

40.7 

15.9 

27.4 

23.8 

41.0 

207.0 

0.072 

15.6 

26.9 

21.0 

36.2 

14.8 

25.5 

20.4 

35.2 

14.8 

25.5 

20.3 

35.0 

1 


Wrinkling  Load 
Faw,  N gw,  Nondim. 


Both  the  theory  and  experiment  show  the  same  general  behavior  of  the  flexibility 
with  pressure.  Much  like  that  found  for  pressure-stabilized  beams,  the  flexibility  of  the 
arch  varies  inversely  with  the  pressure.  The  theory  is  in  better  agreement  with  the 
experimental  results  for  the  inward  force,  and  the  quality  of  the  agreement  ranges  from 
good  over  the  entire  pressure  range  for  the  arches  with  radius  ratios  of  25.0  and  19.0 
to  poor  for  the  other  two  arches  at  the  lower  pressures.  In  general,  agreement  between 
theory  and  experiment  is  better  at  higher  pressures  and  for  larger  radius  ratios.  The  theory 
predicts  the  arches  to  be  more  flexible  than  they  are  found  to  be  experimentally.  This 
is  somewhat  unexpected  because  we  are  dealing  with  fixed  end  restraints,  and  it  is  always 
an  experimental  problem  to  attain  full  fixity  of  the  ends.  Since  fixity  in  the  theory 
can  by  attained  perfectly,  it  is  expected  that  the  theoretical  prediction  will  be  less  flexible 
than  the  experimental  results.  This  anomaly,  along  with  other  differences  between  theory 
and  experiment,  may  result  from  the  treatment  of  the  nonlinear  stress-strain  behavior  of 
the  fiberglass  fabric,  that  is,  linearizing  it  about  the  pressurized  state.  An  additional  factor 
that  could  contribute  to  the  lack  of  agreement  between  theory  and  experiment  is  the 
use  of  straight  woven  tubes  to  measure  the  fabric  moduli  or  stiffnesses  used  in  the 
theoretical  representation  of  the  contoured  arch  behavior. 

Wrinkling  Load 

The  wrinkling  load  results  are  presented  in  Figures  18,  20,  22  and  24  which  show 
the  wrinkling  as  determined  both  experimentally  and  theoretically,  increasing 
monotonically  with  pressure.  In  the  lower  pressure  range  the  increase  in  the  winkling 
load  is  seen  to  be  more  rapid  than  a linear  function,  but  at  higher  pressures  the  behavior 
is  essentially  linear.  It  should  be  noted  that  these  experimental  results  are  for  the  applied 
force  directed  inward,  since  this  loading  puts  the  arch  in  a compressive  state  of  stress, 
thus  making  wrinkling  possible.  Experimental  wrinkling  load  data  were  not  collected  for 
the  arch  having  a radius  ratio  of  25.0  because  of  experimental  difficulties  in  getting  the 
arch  inflated  without  distortion.  The  fiberglass  fabric  was  very  easily  distorted  when  not 
inflated  and  any  such  distortion  had  to  be  corrected  as  the  arch  was  inflated.  This  difficulty 
was  more  pronounced  with  the  arches  having  the  larger  radius  ratios.  For  the  arch  with 
radius  ratio  equal  to  25.0  we  were  able  to  inflate  the  arch  satisfactorily  for  the  flexibility 
tests  but  could  not  inflate  it  without  distortion  for  the  wrinkling  load  determination. 
The  theory  is  in  good  agreement  with  the  experimental  results  in  the  low  and  middle 
pressure  ranges  but  beyond  this  middle  range  the  agreement  becomes  poorer  as  the  pressure 
increases.  It  is  speculated  that,  as  with  the  flexibility,  the  disagreement  is  partly  due 
to  the  linearization  of  the  stress-strain  behavior  of  the  fabric  and  the  use  of  straight 
fiberglass  fabric  tube  to  determine  the  stiffnesses  used  to  represent  the  three-dimensionally 
woven  fabric  arches.  As  will  be  shown  below,  increases  in  stiffness  decrease  the  wrinkling 
load;  therefore,  the  disagreement  between  theory  and  experiment  for  the  wrinkling  load 
is  consistent  with  that  for  the  flexibility.  This  work  was  begun  under  the  assumption 
that  the  mode  of  failure  would  be  wrinkling  and  not  geometric  instability,  a possible 
failure  mode  for  arch  structures.  The  conduct  of  the  experiments  confirmed  this 
assumption.  The  wrinkle  formation  could  not  be  observed  because  it  was  covered  by 
the  strap  used  in  the  load  application.  However,  the  failure  was  a gradual  one  with  no 
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out-of-plane  deformation.  This  is  in  contrast  to  the  character  of  instability  failures  which 
occur  rapidly  and  are  accompanied  by  out-of-plane  deformation.  In  addition,  the  agreement 
between  the  theoretical  prediction  and  the  experimental  determination  of  the  wrinkling 
load  substantiates  the  belief  that  the  experimental  failures  occurred  by  the  wrinkling  mode 
and  not  as  the  result  of  instability. 

It  will  be  recalled  that  the  elastic  modulus  of  the  fiberglass  fabric  used  in  these  arches 
varied  widely  with  the  pressure  parameter.  This  variation  has  an  impact  on  the  load-carrying 
capability  of  the  arches  as  measured  by  the  wrinkling  load.  The  nondimensional  wrinkling 
load  is  computed  from  the  theory  by  the  relation 

9w  = n/cmE<"> 

where  gw  is  the  nondimensional  wrinkling  load,  n is  the  pressure  parameter,  em  is  the 
absolute  value  of  the  maximum  compressive  strain  and  E(n)  is  the  nondimensional  elastic 
modulus,  a function  of  the  pressure  parameter.  If  E(n)  is  a constant,  that  is,  not  a function 
of  pressure,  the  wrinkling  load  increases  with  pressure  by  two  factors;  first,  the  increase 
in  pressure  itself  and  secondly  the  decrease  in  em  as  a result  of  the  increased  stiffness 
caused  by  the  pressure  increase.  When  the  elastic  modulus  is  a function  of  the  pressure 
there  are  two  competing  effects.  The  strain  em  tends  to  decrease  because  of  the  stiffness 
increase  due  to  both  increases  in  pressure  and  the  corresponding  increase  in  E(n).  This 
decrease  in  strain  magnitude  is,  however,  diminished  by  the  increased  magnitude  of  E(n). 
This  situation  is  shown  graphically  in  Figure  25,  where  we  have  a plot  of  the  wrinkling 
load  as  a function  of  pressure.  The  different  curves  in  this  plot  represent  the  behavior 
of  the  wrinkling  load  for  different  rates  of  increase  of  the  elastic  modulus.  The  modulus 
is  a linear  function  of  the  pressure  parameter,  and  the  slope  of  this  linear  function  is 
the  parameter  that  changes  in  Figure  25.  Thus  the  curve  o = 0 represents  the  case  where 
E(n)  is  a constant  and  that  for  a = 62.0  is  the  behavior  of  the  arch  described  in  this 
report.  The  remaining  curve  is  for  intermediate  rate  of  increase.  It  is  clear  from  these 
results  that  the  increase  in  the  stiffness  more  than  overcomes  the  decrease  in  strain  resulting 
in  a net  loss  in  load-carrying  capability.  Thus  to  obtain  the  most  benefit  in  load-carrying 
capability  from  increases  in  pressure,  it  is  desirable  to  have  the  smallest  increase  in  elastic 
stiffness  due  to  pressure  of  the  fabric  used  to  make  the  arch. 

These  results  illustrate  the  deformation  and  load-carrying  capability  of 
pressure-stabilized  arches,  and  the  comparisons  of  the  experimental  and  theoretical  results 
give  confidence  that  the  theory  can  predict  the  bahavior.  The  theory  can  then  serve 
as  a useful  design  tool. 


CONCLUDING  REMARKS 

The  development  of  a linear  theory  for  the  behavior  of  pressure-stabilized  arches 
under  static  load  has  been  presented  and  solutions  to  the  resulting  governing  equations 
obtained.  Solutions  were  obtained  in  the  form  of  a Green's  Function  from  which  solutions 
to  general  loading  can  be  obtained.  In  addition,  a particular  solution  for  the  uniform 
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FIGURE  25.  EFFECT  OF  MATERIAL  STIFFNESS  INCREASE 

WITH  PRESSURE  ON  THE  WRINKLING  LOAD 


normal  load  was  obtained.  Because  of  the  complexity  of  the  solutions,  they  are  not 
written  down  explicitly  but  a computer  program  is  given  to  carry  out  the  solution  and 
compute  the  results  in  terms  of  the  deformation  and  stress  behavior.  An  experimental 
program  to  measure  the  deformation  and  load-carrying  capability  of  a series  of  arches 
is  described  along  with  measurements  of  elastic  and  shear  stiffnesses  of  the  fabric  used 
for  the  arches.  Results  of  this  experimental  program  are  presented  and  compared  with 
predictions  from  the  theory.  These  results  are  given  in  terms  of  the  flexibility  and  wrinkling 
load  which  describe  the  deformation  and  load-carrying  capability  of  the  arches.  The 
flexibility  and  wrinkling  load  are  presented  for  inflation  pressures  ranging  from  34.5  to 
207.0  kPa  for  four  arches  having  radius  ratios  of  11.3,  13.0,  19.0,  and  25.0.  The 
agreement  between  the  experimental  and  theoretical  results  establishes  the  theory  as  being 
adequate  for  the  prediction  of  the  deformation  and  load-carrying  capability  of 
pressure-stabilized  arches.  Thus  the  study  can  provide  a rational  basis  for  the  structural 
design  of  shelters  using  pressure-stabilized  arches. 
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APPENOIX  A 


DEVELOPMENT  OF  THE  ONE-DIMENSIONAL  ENERGY  PRINCIPLE 

The  two-dimensional  energy  principle  in  terms  of  the  cross-section  displacements  and 
rotation  is 


ft  <4 

of  O l 


C22/a|  [(U'  - W)  - 0yacos(0,)  - 


+ d>x)asin(0 , )] 2 + C33/aHa(0x  ~ ^ + 

(W'  + U + R0y)sin(0, ) + (V'  - R0z)cos(0,)]2  + 

Njj/aiKU  + W')cos(0,)  - V'sin(0j)  - 0yacos2(0,)  - 

(A  1 ) 

<f>zasin(0,  )cos(0, )] 2 - 2F|[a0x  + Wsin(0,)  + Vcos(0,)]  - 
2Fj[U  - 0yacos|0,)  - 0zasin(0,))  - 
2F3  (Vsin(0 , ) - Wcos(0,  )]^  a2ad0|  d02 
Expanding  this  expression  we  obtain 


sff  1/2  I — (<u'  - W)2  + (a0y)2cos2  (0, ) + 

of  o \R2(1  + a/Rcos(0, )) 


(a0z  + a0x)Jsin2  (0, ) - 2(U'  - W)0yacos (0,) 

-2(U'  - W)(0Z  + 0x)asin(0,)  + 2 0y(0z  + 0x)a2sin(0,  )cos(0, )]  + 


[(a <t>'  - a0z)2  + (W'  + U + R0v)2sin(0, ) + 

R2(1  + a/Rcos(0,)l  V 

(V'  - R0z)2cos2(0,)  + 2(W'  + U + R0y)(0x  - 0z)asin(0t ) + 

2(0X  - 0Z)|V'  - R0z)acos(0, ) + 2(W*  + U + R0y)(V'  - R0z)cos(0,  )sin(0, )] 


(A2) 

+ 
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N' 


1 

O 

21 [(U  + W')2cos2(0,)  + (V')*sin2 (Ox)  + (0va)2 cos4 (0 , ) + 

R2  (1  + a/ pcos(0 1 )) 

(a0z)2sin2  (0,  )cos2  (0 ( ) - 2(U  + W')V'cos(0,  )sin(0, ) - 
2(U  + W')0yacos3(0, ) - 2(U  + W')0zasin(0,  )cos2  (0, ) + 

2V'0yasin(0 1 )cos2  (0, ) + 2V'</>  zasin2  (0  ( )cos(0| ) + 

20y0za2sin(0,  )cos3(0, ))  - 2(1  + a/pcos(0,  ))[UF2  + 

W(F,sin(0,)  - F3cos(0,))  + V(F,cos(0,)  + F3sin(0 , ))  + 

#xaF,  - <£yaF2cos(0, ) - </>zaF2sin(0,  )]^  aRd0,d02  = 0 

The  0,  dependence  is  all  in  explicit  form  in  equation  (A2),  and  we  can  thus  evaluate 
the  integrals  with  respect  to  0,  and  this  is  done  as  follows,  starting  with 


/ 


2tt 

o 


(a/p)4COS4(0,) 

"1  + (d/p)cos(0, ) 


d0. 


-"V" z 


(A3) 


which  is  a definition  of  Z for  which  a formula  will  be  given  later.  It  should  be  noted 
that  Z has  a nonvanishing  limit  as  a/p  becomes  very  small.  To  aid  in  the  evaluation 
of  the  remaining  integrals  note  that  (A3)  can  be  written  with  the  aid  of  the  binomial 
expansion  as: 


-2tt  (a/R)4cos4(01 ) 

J Q 1 + (a/p)cos(0,)  1 


[ ( — )4  COS4  (0 1 ) 

R 


( — )5cos5  (0j ) + (l)6cos6(0,)  - . . .]d0, 
R R 


(A4) 


Taking  next  the  integral  having  the  third  power  of  the  cosine  and  using  the  binomial 
expansion  we  obtain: 


rv  (a/p)3cos3(0, ) d0, 
O 1 + (a/p)cos(0,) 


( — )3cos3(0,  )d0, 
R 


( 1 j6cos6(0, ) - . . .] d0 , = -7r(  — )4Z 
K R 


[(  l)4cos4  (0,)  - ( — )5 cos*  (0 j ) + 

R R 


(A5) 
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Thus  all  the  integrals  can  be  evaluated  in  terms  of  the  parameter  Z,  and  this  parameter 
can  be  evaluated  by  the  direct  integration  of  (A8)  as: 


f. 


1 + (a/R)cos(0,) 


d0,  = 2tt[1  - (*  )*r,/a 

R 


(A14) 


Equating  this  result  with  that  of  (A8)  gives 


Z = ( — )4 [2[1  - (I)2]-,/2 

a R 


2 - (_)2] 

R 


(A  15) 


Using  these  integral  formulas  in  (A2)  and  rearranging,  we  obtain  the  one-dimensional  energy 
principle 


s/"2  (5 


[2(U'  - W)2  + Z2(1)2(U'  - W + R0')3  + a3Z 3(<t>'  + <t>x)2]  + 
R v * 


'3  3 


fF  [2a 


!(0^  - 02)3  + Z2(_)3[a(^  - tf2)  - _(V'  - R0Z)] 3 + 
R a 

N°22 


Z3(W'  + U - R0J2]  + — - [(U  + W')  + 
y R2 


Z<_)2(U  + W'  + R0J2  + (\/')2  + 


Z4(-  )2(V'  - R0J2]  - 2 [ Uf 2 + Wf,  + 
R 


(A16) 


Vf3  + a<M4  - a0vf5 


- a0zf6]  | 


aRd02  = 0 


The  generalized  force  components  associated  with  the  cross-section  displacement  and 
rotations  are  defined  by 


1 #2* 

f,  =-[  (F, sin(0, ) - F3cos(0,))(1  + (a/R)cos(0,  ))d0, 

nJ  o 


(A  1 7) 
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U = - / Fj(1  + (a/RjcoslMdfl, 

n J o 


(F,  co$(0, ) + F3sin(0,  >)(  1 + (a/R)cos(0,  ))d0 
F,  (1  + (a/R)cos(0,  ))d0, 

F2cos(0,)(1  + (a/  R)cos(0 1 ))d0 1 
F2  sin(0 1 )( 1 + (a/R)cos(0 1 ))d0 , 


1 


This  completes  the  detail  development  of  the  one-dimensional  energy  principle. 
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APPENDIX  B 


DERIVATION  OF  ONE-DIMENSIONAL  STRESS  RESULTANTS  AND  MOMENTS 

The  one-dimensional  stress  resultants  and  moments  are  defined  as  integrals  of  N' 
and  N,  j over  the  cross-section  and  are  expressed  in  terms  of  the  cross-section  displacement 
and  rotations  by  use  of  the  stress/strain  law  and  the  strain  displacement  relations  which 
are  repeated  here  for  convenience. 

n;,  - c„.;,  (B1| 

Ni  2 = 2C33ei  2 /Qoi 


e'22  = . — MU'  — W ) - 0'acos(0,)  - 

R(1  + (a/R)cos(0,) 

+ </>x)asin(0, )] 


c/rii«5m  Ial**  ' * + lw' + u * * 


(V'  — R0z)cos(0,)] 


f 21  = 


2R2(1  + (a/ R)cos(0 1 ))‘ 


MU  + W')cos(0, ) - V'sin(0, ) - 


0yacos2  (0, ) - 0zasin(0,  )cos(0, )] 2 
The  axial  stress  resultant  is  defined  as: 


* n 


t =/  c22  - [ 


(U'  - W) 


0yacos(0, 


R 1 + (a/R)cos(0,)  1 + (a/R)cos(0,) 


< <t>'z  + 0x>asin(0|) 

1 + (a/R)cos(0,) 


80 


Using  (A5)  through  (A  13)  for  the  evaluation  of  the  integrals  we  have 


t = ttC22  1 [(2  + d)2Z2)(U'  - W)  + (±  )2Z2  R0' ] 

R R R y 

The  moments  about  the  x,  y,  and  z axes,  respectively,  are  defined  as: 


.2* 

mx 

« 

fo 

N',  2a2d0, 

f2ir 

my  = 

7 « 

(o 

N'2  2a2cos(0, 

)d0 

r2* 

mz  = 

1 

L 

N'22a2sin(0| 

)d0 

Substituting  the  stress/strain  law  and  the  strain  displacement  relations  we  have: 


nv 


a(4>'x  - <pz) 

1 + (a/R)cos(0,) 


(V'  — R</>z)cos(0, ) 
(1  + (a/R)cos(0,) 


m, 


(U'  - W)cos(0,) 

[ 

1 + (a/R)cos(0,) 


m. 


a 2 (U'  - W)sin(0,) 

C2  2 — t 

R 1 + la/R)cos(0,) 


IW'  + U + R<*> v)sin(0, ) 

] d0, 

1 + (a/R)cos(0 , ) 

0yacos2(0, ) 

1 + (a/R  )cos(0 , ) 

(<t>'z  + 0x)acos(0,  )sin(0, ) 

]d0, 

1 + Ia/R)cos(0|) 

0yacos(0,  )sin(0, ) 

1 + (a/R)cos(0, ) 


(0z  + 0x)asinJ  (0, ) 

] d0, 

1 + (a/R )cos(0 1 ) 
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(B6) 


I 


1 


Carrying  out  the  integral  with  respect  to  0,  using  the  formulas  (A5)  through  (A  13)  we 
have 


mx  = *c33  -2  tR(2  + ( — )2Z2)[<p'x  + 0Z)  - Z(V'  - R0  )] 
R ‘ R 

a3 

m = -7rC22Z2  _ [U'  - W + R0'] 
y Rz  y 


m = -ttC22Z3  - [R(0*  + 0 )] 
R2 


(B7) 

(B8) 

(B9) 


The  remaining  two  stress  resultants  are  the  transverse  shear  resultant  defined  by: 
•2tt 

[N',  2sin(0, ) + N°22  J2e'i2  cos(0,)]ad0, 


f 

f. 


[N',  2cos(0 1 ) - N22  sin(0,)ad0, 


Substituting  the  stress/strain  law  and  the  strain  displacement  relations  we  have: 


f2  a »x 

a 'J  0l  " TT 


- 0z)asin(0,)  (V'  - R0z)cos(0,  )sin(0, ) 


+ .. 


(a/p)cos  (0, 


1 + (a/R)cos(0,) 


(W'  + U + R0jsin20  0 a (u  + W')cos2(0.) 

V ] + N^2  _ [. 

1 + (a/R)cos(0 , ) 

V'sin(0,  )cos(0 1 ) 


R 1 + (a/R)cos(0 , ) 
d7asin(0,  )cos2(0. 


0yacos3(0, ) 

1 + (a/R)cos(0,) 


1 + (e/R)cos(0, ) 1 + (a/R)cos(0, ) 


1 


d0, 
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h. 


f.HZ 


- <pz) acos(0,)  (W'  + U + R0v)cos(0,  )sin(0, ) 


(a/R)cos(0, ) 


1 + (a/  R)cos(0 1 ) 


(V'  - R0_)cos 2 «?, ) 


- N^2  _ 


(U  + W')cos(0,  )sin(0, ) V'sin2(0,) 


1 + (°/R)cos(0 1 ) R 1 + (a/R)cos(0, ) 1 + (a/R)cos(0, ) 


0yacos2(0j  )sin(0, ! 


0zasin2  (0,  )cos(0, ) 


1 + (a/R)cos(0,)  1 + (a/R)cos(0 , ) 


j d0, 


Again  we  integrate  with  respect  to  0,  using  (A5)  through  (A  13)  to  give 
qz  = 7rC33Z3  1 (W'  + U + R0j  + 


(810) 


*Nf2  _ [Z2(U  + W')  + (—  )2ZR0v] 
R R y 


qy  = *c33z2 1 «v'  - R0Z>  - (1)2r(^  _ 0z)]  _ 


’rNa:  - (~Z3V’  + (l)2RZ40z) 
R R z 


(B11) 


These  stress  resultants  and  moments  are  shown  graphically  in  Figure  2.  In  addition  to 
being  the  definitions  of  the  stress  resultants  and  moments,  the  above  expressions  are  also 
the  force  displacement  relations  for  the  one-dimensional  theory. 
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APPENDIX  C 


THEORY  IN  NONDIMENSIONAL  FORM 


When  obtaining  solutions  and  results  from  theory  it  is  useful  to  have  the  equations, 
including  the  basic  definitions,  in  nondimensional  form.  To  do  this  for  the  present  theory, 
we  begin  by  defining  nondimensional  displacements. 


u = 

U/a 

(Cl) 

w = 

W/a 

(C2) 

V = 

Mia 

(C3) 

The  nondimensional  strain  measures  are 

obtained  from  equations  (18)  as 

E 

= ((/ 

- W)/p 

(C4) 

Ky 

= R«y 

= *y 

(C5) 

= R*c2 

= 0k  + K 

(C6) 

Kx 

X 

id 

cr 

n 

n 

■o- 

X ' 

1 

(C7) 

rz 

= ( u ^ 

• IV'  + P0y)/P 

(C8) 

rv 

= (V1 

- P0z)/P 

(C9) 

where  p = R/a.  The  nondimensional  stress  resultant  and  moments  are  defined  as  follows, 
and  the  nondimensional  stress-displacement  relations  follow  from  equations  (16): 


T = t/C22a  = (rrd/p) ((2  + Z 2/p2)(U'  ~ W)  + Z 2 0y/p]  (CIO) 

Qv  = q V/C22a  = (wcZ2  lp)[[V  - p<t>z)  - (<t>'x  - 4>z)/p]  - 

v y (Cl  1) 

lirn/p)l-Z3V  + Z 4<t>z/p] 

Q,  = q7/C22a  = {ircZ3/p){W'  + U + p</>v]  + 

(Cl  2) 

(irn  lp)[Z2(U  + W')  + Z0y/p] 

Mx  = mx/C22a2  = (ttc/p)  ((2  + Z2/p2)(<p'x  - <t>z)  - Z2(V  - p<t>z)/p]  (C13) 
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My  = my/C22a2  = -(?rd Z2/p2)[U'  - W + p0y]  (C14) 

Mz  = mz/C22a2  = -|?rdZ3/p2)l0z  + 0X)  (Cl  5) 

The  equilibrium  equations  in  terms  of  the  nondimensional  stress  resultants  and  moments 
are  obtained  by  dividing  equations  (20)  by  C22a  to  yield 

-T'  + Q2  - 7rf2  = 0 (Cl 6) 

-(T  + Q'z)  — ttT,  = 0 (Cl  7) 

My/p  + Qz  - (nn/p)(U  + W')  + 7rTs /p  = 0 (Cl 8) 

-Qy  - jtF j=0  (C19) 

(Mz  - Mx)/p  - Qy  + (irn/p)V'  + nJ6/p  = 0 (C20) 

-(Mz  + Mx)/p  - nTjp  = 0 (C21) 

The  nondimensional  forces  are  given  as  7j  = Rf,/C22  and  n = N22/C22  is  the 


nondimensional  stress  due  to  the  internal  pressure.  The  parameters  d=  C22/C22  and 
c = C33/C22  are  the  nondimensional  stiffnesses  and  C22  is  a reference  value  of  C22. 
The  governing  equations  in  terms  of  the  nondimensional  displacement  parameters  are 
obtained  from  equations  (22)  by  dividing  by  C22a  to  give 

-[2d/p  + dZ2/p3]  U"  + fcZ3 /p  + nZ2/p]  U - [dZ2/pJ]0y  + 

(C22) 

[cZ3  + nZ/p2]0y  + [2d/p  + dZ2/p3  + c Z3/p  + nZ2  /p]W'  - f2  = 0 
“IcZ 3/p  + nZ2  /p]W"  + [d(2  + Z2/p2)/p]W  - [dZ2/p2  + cZ3  + nZ/p2]0y  - 

(C23) 

[d(2  + Z2/p2)/p  + cZ3/p  + nZ2/p]U'  - T,  = 0 

— [dZ2/p]  0y  + [cZ3p  + nZ/p]  0y  - [dZ2/p2  ] U"  + (cZ3  + nZ /p2]U  + 

(C24) 

[dZ2/p‘  + cZ3  + nZ/p2]lV'  + Ts  = 0 
— [cZ2/p  + nZ3/p]  V"  + [cZ2  /p2  ] 0X  + 

(C25) 

[cZ2  (1  - 1/p2)  + nZ4/p2 ] <t>'z  - T3  = 0 

-IdZ 3/p2]0z  + [c(2/p2  + Z2(1  - 1/p2)2J+  nZ4/p2]<fiz  - 

[dZ3/p2  + c|2  - Z2  (1  - 1/p2))/p2]0x  - (C26) 

[cZ2  (1  - 1/p2)/p  - n(1  - Z3 )/p)  V"  + Tjp  = 0 
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(C27) 


— fc(2  + Z 2/p2)/p2]0"  + [dZ3/p2]0x  + [cZ2  /p3]  V"  + 

[dZ3/p2  + c|2  - Z2(1  - 1/p2))/p2]^  - Tjp  = 0 

The  boundary  conditions  are  similarly  expressed  in  terms  of  the  nondimensional 
displacement  parameter  or  the  nondimensional  stress  resultants  and  moments.  The  Green's 
function  discontinuity  conditions  are  put  in  nondimensional  form  by  dividing  the  stress 
resultant  equations  by  aC22  and  the  moment  equations  by  a2C22  to  give 


T,  — T2  = 7t§2 

(C28) 

My,  - Myj  = 7rgs 

(C29) 

Oz,  - QZJ  = 7rg, 

(C30) 

MZ1  ~ MZ2  = *96 

(C31) 

MX,  _ ^X2  = ^54 

(C32) 

Qyi  ~ Oy2  = 

(C33) 

where  the  nondimensional  forces  are  g-,  = g,/C22(  i = 1,6. 
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APPENDIX  D 

RELATIONSHIP  AMONG  THE  CONSTANTS  OF  INTEGRATION 

The  homogenous  solution  as  expressed  by  equations  (37)  contains  18  unknown 
constants.  But,  because  the  characteristic  numbers  cause  the  determinate  of  equation  (33) 
to  vanish,  there  is  for  each  of  these  characteristic  numbers  a relationship  among  these 
constants  which  is  used  here  to  express  the  Aj  and  Bj  in  terms  of  the  C;.  To  obtain 
these  relationships  the  solution  (37)  is  substituted  into  the  governing  differential 
equations  (30)  and  causing  the  coefficients  of  the  functions  cosh(X02 ),  sinh(X02 ),  cos(02 ), 
sin(02),  0 2 cos(0 2 ) and  02sin(02)  to  vanish  independently  yielding  18  homogenous 
equations  in  18  unknowns.  This  set  of  equations,  however,  has  a rank  of  only  12,  thus  12 
of  the  unknowns  can  be  expressed  in  terms  of  the  remaining  6.  The  18  equations  in 
matrix  form  take  the  form  shown  on  the  following  page 

where 

J,  = H2  - X2H, 

J2  = H4  - X2  H3 

J3  = H,  + H2  (2D) 

J4  = H3  + H4 

J5  = H,  - X2H2 

Examination  of  these  equations  reveals  that  the  relation  among  the  constants  A , 

B[  and  C|  is  the  same  as  that  among  A2,  B2,  and  C2;  thus  we  obtain  from  two  of 
the  three  equations 


and 


I 
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where 


Similarly,  the  relation  among  A3,  B3,  A6,  B6,  C4,  and  C5  is  the  same  as  that  among 
A4,  B4,  -As,  -B5,  — C3 , and  C6;  thus  we  obtain  from  equations  (7),  (8),  (10),  and  (11) 


where 
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APPENDIX  E 

SYSTEMS  OF  EQUATIONS  FOR  DETERMINATION  OF  THE 
CONSTANTS  OF  INTEGRATION 

equations  C°effidents  of  ^ems  of 

equations.  The  coefficients  for  the  Grepn'c  f * ^ *°r  the  solution  °f  tfle  arch 

for  the  uniform  “ad  “ton.  b*  ""*■  »v  -hose 

are  required  for  ’the'fl1' Snrtrata.  ' Si»°of,ath^e°IoSB  ' T * T’SMom 

solution  over  -a  < 0,  < 8 anH  ,h»  • cons  . ts-  C‘  trough  C6,  define  the 

5>  <».  <«  “““  - 

IS]  C = b 

r“n“  °'  ime9ra“<"’  a"d  ,h'  °<  ««*!»«.  S are 

p sented  below  The  first  three  equations  impose  the  boundary  conditions  aT  n - n 

and  equations  (4)  through  (6)  impose  the  boundary  conditions at  fl  - J r V ■ ' 

F>  ttr  5 ~ srsf?r  Eft: 

> . ' ' ' are  evaluated  with  ease  by  evaluation  of  U IV  ar>H  a,  ,*  a 

or  the  fixed  boundary  condition  and  Ult  Wt  and  7u!  - W ‘l  ot  \ In  = “* 

or  the  simply  supported  boundary  conditions.  For  s--  i = 4 6-  i - ItoVL  °2  = ““ 

U J(  Wj,  and  0 replace  U W ann  r*  ° " 1 '6' 1 ~ 7-12>  the  Parameters 

u 2 vV*  piace  t/"  wt>  and  ^y,  and  are  evaluated  at  «,  = „ ,n  ... 

we  have  s;;  = 0 for  i = ir  i - 710  ■ „ „ . a-  ,n  addition, 

171  through  ',9>  impose  continuity  A displacement  a'.  I,  %'  L sT-lT 

fi  - 7.9;  i- 7.12,)  :are  found  by  evaluating'!/,,  h^'and  *'  i - » Vh  ^ 

of  the  last  three  equations  are  given  in  detail  belou,  v ’ ' The  '««'cients 


Si  0 1 = L,cosh(X02) 

5102  = L,sinh(X02) 

5103  = LjsinOj ) 

s10>4  = Ljcos(^j) 

si05  = L3cos(#2)  + L202sin(02) 
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1 

Sio(6  = -LjsinlS'j)  + L2(r2cos(?2) 

L,  = (j3j A - 1)H,  + H302X 

L2  = |/33  - 1)H,  + H304 

L3  = (P5  ~ P3  )H j + (06  - 04)H3 

Si  i i = L4cosh(A02) 

s,  i 2 = L4sinh(A<72) 

Si  i 3 = Ls sin(02 ) 

Sl  l 4 = L5cos(B2) 

S|  1 5 = L6cos(52)  + Ls02sin(02) 

S|  i 6 = — L6sin(02 ) + Ls02cos($2  ) 

/ 

L4  = [«3j  A - 1)  + p/32  A]  H3 
Ls  = W3  - D + P/34]  H3 
L6  = Ms  ~ 03  > + PiPe  ~ P4 ) I H 3 
s,  2 1 = L7sinh(X$2) 

1 

S,  2 2 = L7cosh(A02) 

Si2(3  = L8COS(52) 

Si  2 4 = -LgSin(^j) 

St  2 5 = L9sin(52)  + L802cos«F2) 

Si2,6  = L9cos(02)  - L852sin(^2) 

L7  = (X  + Pt  )H2  + H4/32 

L8  = (1  - 03)H2  - H404 

L,  = |1  + /3s)H2  + H406 
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The  last  six  colums  of  these  three  equations  are  given  by  the  relation 


si,  j + 6 = -si,j 

i = 10,12 
j = 1.6 

The  system  of  equations  is  completed  by  specification  of  the  elements  of  the  force  matrix, 
b,  as  follows 

bj  =0  i - 1,9 

b| o = 92 
b\  i = —9s 
bi  2 = -S. 

where  g, , g2,  g5  are  the  nondimensional  applied  forces. 

For  the  uniform  load,  the  symmetry  of  the  structure  and  the  load  about  02  = 0 
makes  it  possible  to  reduce  the  number  of  constants  of  integration  to  three,  which  are 
determined  from  the  boundary  conditions.  The  equation  for  these  three  constants  can 
be  written  in  matrix  form  as: 


sl  l 

Sl  2 

S1  3 

V 

0 

S2  1 

S2  2 

s2  3 

) 

< 

c4 

> = < 

0 

S3  1 

S3  2 

s3  3 

C5 

L J 

pf,/d(2  + Z2/p2) 

J 

For  fixed  boundary  conditions  the  coefficients  are 

s,  , *>  /3,  sinh(Xor) 

Si  2 = /33 sin(or) 

s,  3 = 0ssin(a)  - 03cw:os(a) 

s2 1 = 02sinh(Xa) 

s22  = ft,  sin  (a) 

S23  = 06sin(a)  - 04acos(a) 
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s3,  = cosh(Xo) 


s3  2 = cos(a) 
s33  = asin(a) 

and  for  simply  supported  boundary  conditions,  the  second  equation  is  changed  as  follows 
to  set  the  moment  equal  to  zero. 

s2.  = ((0.X  - 1)H,  + H302X]cosh(Xa) 

s22  = 1(03  - 1)H,  + H304]  cos(a) 

s23  = 1(05  - 03)H,  + (06  - 04)H3]  cos(a)  + 

1(03  - 1)H,  + H304]asin(a) 
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COMPUTER  PROGRAM  FOR  SOLUTION  OF  GOVERNING  EQUATIONS 

Our  purpose  here  is  to  describe  the  computer  program  which  carries  out  the 
computations  involved  to  obtain  specific  results  from  the  solution  described  in  the  body 
of  the  report  and  Appendix  E. 

General  Organization  of  Program 

The  flow  of  calculations  in  the  program,  a copy  of  which  is  given  at  the  end  of 
the  appendix,  is  illustrated  by  the  flow  chart  in  Figure  F-1.  This  flow  chart  is  not 
intended  to  depict  the  details  of  the  calculation  but  only  the  general  organization  or 
sequence  of  the  calculations.  The  program  includes  a main  program  which  calls 
2 subroutines  peculiar  to  this  program,  all  of  which  carry  out  computations  in  double 
precision.  The  program  also  calls  a system  subroutine  for  the  solution  of  system  of  linear 
equation.  Because  the  program  is  in  double  precision  this  system  subroutine  must  also 
be  double  precision.  The  program  was  written  for  the  UNI  VAC  1106  computer  with 
the  EXEC-8  operating  system  and  the  system  subroutine  called  is  associated  with  this 
computer  system.  It  is  believed  that  the  program  written  in  FORTRAN  could  be  easily 
adapted  to  other  systems. 

The  program  contains  options  for  the  type  of  loading,  the  boundary  conditions  and 
plotting  the  output.  The  loading  option  includes  a uniform  normal  loaded,  a concentrated 
load  either  normal  to  the  arch  or  vertical  and  at  an  arbitrary  position,  and  a general 
distributed  load.  The  concentrated  load  solution  is  a Green's  function  which  is  used  to 
compute  the  solution  for  the  general  distributed  load  by  quadrature.  The  functional  form 
of  the  distributed  load  must  be  specified  in  SUBROUTINE  CALPAR.  The  quadrature 
for  the  general  load  is  carried  out  using  the  Gaussian  numerical  technique.  The  coordinates 
and  weighting  factors  for  this  technique  must  be  on  logical  unit  10  and  written  according 
to  FORMAT  (2E18.8).  Two  boundary  condition  options  are  provided.  One  has  both 
ends  fixed  and  the  other  both  ends  simply  supported.  Other  options  could  be  included 
by  modification  of  the  program.  The  plotting  option  when  exercised  provides  for  the 
plotting  of  the  computed  flexibility  and  wrinkling  load  as  a function  of  the  pressure 
parameter.  In  addition,  experimental  results  will  be  plotted  along  with  the  corresponding 
calculated  results.  This  experimental  data  must  be  located  on  logical  unit  16  and  written 
according  to  FORMAT  (5E15.6).  The  CALCOMP  plotter  is  used,  and  it  operates  a plot 
tape  which  must  be  assigned  as  logical  unit  20.  The  program  also  uses  a scratch  file 
which  must  be  assigned  as  logical  unit  15. 

Input  and  Execution 

The  input,  both  control  parameters  and  data,  required  to  operate  the  program  is 
described  in  this  section. 
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INPUT 


CYCLE  ON  PRESSURE  INCREMENT 


CYCLE  ON  RADIUS  RATIO  INCREMENT 


COMPUTE  CHARACTERISTICS  VALUES 
AND  ASSOCIATED  CONSTANTS 


GENERATE  AND  SOLVE  EQUATIONS  FOR  CONSTANTS 
OF  INTEGRATION  IN  ACCORDANCE  WITH  80UNDARY 
CONDITION  CODE 


COMPUTE  FOR  OUTPUT  THE  STRAINS  STRESSES  AND 
DISPLACEMENTS 


SUBROUTINE  CALCA 


SUBROUTINE  DGJR 


SUBROUTINE  CALPAR 


IF  INTEG-O 

TyeT" 


CYCLE  ON  INTEGRATION  INCREMENT 
FOR  GENERAL  LOAD 


READ  e AND  ASSOCIATED  GAUSSIAN 
WEIGHTS  FOR  INTEGRATION 


GENERATE  AND  SOLVE  EQUATIONS  FOR  CONSTANTS 
OF  INTEGRATION  IN  ACCORDANCE  WITH  BOUNDARY 
CONDITIONS  AND  CONCENTRATED  LOAD  POSITION 


<T1mnteg-o'"'>-^ 


YES 

CARRY  OUT  STEP  IN  INTEGRATION  OF 

THE  GREEN  S FUNCTION  FOR  THE 
CURRENT  VALUE  OF  9 

<^^INTEGRATION  CYCIe'''^>7 



Tyes 

COMPUTE  FOR  OUTPUT  THE  STRAINS  STRESSES  ANO 
DISPLACEMENTS 

1 SUBROUTINE  CALPAR  | 

IF  RADIUS  RATIO 
CYCLE  COMPLETE 


IF  PRESSURE 
CYCLE  COMPLETE 


FIGURE  F I FLOW  CHART  OF  ARCH  ANALYSIS  PROGRAM 
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CARD  #1 


Format  (14) 


Parameter 

Column 

Description 

IP  LOT 

1 

plot  control  IPLOT  = 1 output  is  not  plotted 

IP  LOT  = 0 output  is  plotted 

ILOAD 

2 

load  direction  control  ILOAD  = 1 vertical  load 

ILOAD  = 0 normal  load 

this  parameter  does  not  apply  to  the  uniform 
load 

INTEG 

3 

Integration  control  INTEG  = 1 general  load 

INTEG  = 0 uniform  and 

concentrated  load 

IBOUN 

CARD  #2 

Format  (7F10.0) 

4 

boundary  condition  control  IBOUN  = 0 fixed 

IBOUN  = 1 simply 
supported 

Parameter 

Column 

Description 

RLIM(1) 

1-10 

initial  value  of  the  radius  ratio 

RLIM(2) 

11-20 

final  value  of  the  radius  ratio 

RLIM(3) 

21-30 

radius  ratio  increment 

XNLIM(I) 

31-40 

initial  value  of  the  pressure  parameter 

XNLIM(2) 

41-50 

final  value  of  the  pressure  parameter 

XNLIMI3) 

51-60 

pressure  parameter  increment 
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CARD  #3 
Format  (3F10.0) 


Parameter 

Column 

Description 

ALPHA 

1-10 

half  span  angle  of  the  arch 

THETA 

11-20 

location  of  concentrated  load  on  the  arch 

THEHAT 

21-30 

location  of  load  on  the  cross-section,  see 
equation  (43) 

CARD  #4 

Format  (2F10.0,  15) 

Parameter 

Column 

Description 

F 

1-10 

intensity  of  uniform  load;  F = 0 causes  the 
uniform  load  calculation  to  be  omitted 

SMG 

11-20 

concentrated  load  magnitude;  SMG  - 0 causes  the 
concentrated  load  computation  to  be  omitted. 

For  the  general  distributed  load  SMG  must  be 
put  equal  to  unity 

NMAX 

21-25 

number  of  output  intervals;  if  NMAX  - 18  and  the 
arch  spans  180°,  program  will  print  output 
every  10°. 

CARD  #5 

Format  (7 FI 0.7) 

Parameter 

Column 

Description 

AC 

1-10 

the  program  allows  the  nondimensional  shear 

BC 

11-20 

modulus,  c,  and  the  nondimensional  elastic 

AE 

21-30 

modulus,  d,  to  be  linear  functions  of  the 

BE 

31-40 

pressure  parameter,  XN.  These  fj%ir  parameters 

define  these  linear  functions  as 
c - AC'XN  + BC 


d = AE'XN  + BE 
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CARD  #6 


Description 

an  identification  of  the  job  to  be  printed  out 
with  the  output 

An  example  of  the  coded  input  is  given  in  Table  F-1.  This  input  was  used  to 
generate  the  sample  output  to  be  discussed  subsequently.  The  execution  of  the  program 
on  the  UNI  VAC  1106  computer  with  the  EXEC-8  operating  system  can  be  accomplished 
with  the  runstream  shown  on  page  107.  This  includes  reduction  of  experimental  data 
for  plotting,  computation  of  Gaussian  integration  constants,  assignment  of  the  required 
files,  compilation  and  assembly  of  the  arch  analysis  program,  execution  of  the  program 
and  input  data.  The  computer  output  resulting  from  the  execution  of  this  runstream 
is  present  on  pages  108  to  144.  The  runstream  begins  with  the  assignment  of  tape 
no.  1962  which  contains  the  raw  experimental  data  presented  in  this  report  and  a program 
entitled  NONDIM  which  reads  the  raw  data  from  file  15,  reduces  this  data  to 
nondimensional  form,  and  writes  the  nondimensional  data  on  file  16.  If  it  is  desired 
to  plot  different  flexibility  and  wrinkling  load  data  against  pressure,  the  raw  data  must 
be  entered  into  file  15  before  execution  of  the  data  reduction  program.  If  no  plotting 
is  desired,  statements  1-16  of  the  runstream  may  be  omitted.  Upon  execution  of  the 
data  reduction  program,  the  reference  value  of  the  elastic  modulus  is  called  for  as  input. 
This  input  is  the  number  313.0  appearing  as  statement  12  in  the  runstream.  This  is 
the  English  unit  system  equivalent  of  the  number  given  in  the  body  of  the  report  and 
is  given  in  these  units  because  the  raw  data  is  in  English  units.  This  number  must  also 
be  changes  for  a different  set  of  data.  If  plots  are  to  be  made,  this  experimental  data 
must  be  present.  Following  the  data  reduction  we  have  the  generation  of  the  Gaussian 
integration  constants.  This  is  accomplished  with  a short  FORTRAN  program  which  calls 
the  subroutine  LGAUSS  to  generate  the  constants,  and  upon  return,  writes  them  on  file  10. 
As  presented,  90  integration  increments  are  used.  A change  in  the  number  of  increments 
is  accomplished  by  changing  the  parameter  N on  line  20  of  the  runstream.  The  subroutine 
LGAUSS  computes  the  integration  coordinates  over  the  range  -1  to  1.  These  coordinates 
are  transformed  to  the  range  —a  to  a at  the  time  they  are  used  in  the  arch  analysis. 
If  either  a concentrated  load  or  a uniform  load  is  being  treated  which  does  not  require 
integration,  then  statements  17  through  32  may  be  omitted  from  the  runstream.  After 
generation  of  the  integration  constants  tape  no.  1964  is  called  and  the  arch  analysis 
programs  are  read  into  a mass  storage  fil$v<  These  programs  are  compiled  and  collected 
using  the  MAP  processor.  The  absolute  executable  program  is  put  in  B.  If  a general 
distributed  load  is  being  treated,  its  functional  form  must  be  specified  in  subroutine 
CALPAR  which  is  the  symbolic  element  TPF$.ARCHRESULTS.  This  must  be  done  before 
compilation  of  this  subroutine.  If  plots  are  to  be  made,  a plot  tape  assigned  as  file  20 
must  be  assigned  before  execution  of  the  program. 


Format  (IX,  13A6) 

Parameter  Column 

QIDEN  2-79 


104 


Output 


A sample  of  the  output  generated  by  the  program  is  given  in  pages  141  to  144. 
This  output  corresponds  to  the  input  given  in  Table  FI.  As  is  seen,  the  heading  contains 
a title  followed  by  specification  of  the  type  of  loading  and  the  boundary  conditions. 
The  user  provided-problem  description,  the  linear  expressions  for  the  shear  and  elastic 
moduli  and  the  problem  parameters  then  appear.  The  main  body  of  the  output  consists 
of  the  normal  and  longitudinal  deflection,  the  rotation,  the  strains  at  the  outer  and  inner 
arch  radii,  the  axial  force,  and  the  moment  for  a series  of  equal  spaced  angular  positions. 
As  can  be  seen,  these  positions  begin  with  the  smallest  and  increase  to  the  largest  position. 
All  of  the  output  are  in  nondimensional  form. 
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ECS* JUNK ( 1 ) . RSAREP 

t PASG.T  TP..6C9.1962  . TAPE  WITH  DATA  & REDUCTION  PROGRAM 

2 ©MOVE  TP. 2 

3 PCOPIN  TP.. 

4 PFREE  TP. 

5 PASG.T  15..F14  . SCRATCH  FILE 

6 PASG.T  16..F14  . FILE  WITH  EXPERIMENTAL  DATA  FOR  PLOTTING 

7 PD AT  A , I L 15. 

8 PADD.D  TPFS.RAWADA 

9 PEND 

10  PFOR.S  TPF$ . NONDIM, N 

11  PXQT 

12  313.0 

13  PDATA.L  16. 

4 PENO 

15  PERS 

16  PERS  15. 

17  PASG.T  10..F14  . FILE  WITH  GAUSSIAN  INTEGRATION  CONSTANTS 

18  PFOR.IS 

19  DIMENSION  X( 100) ,W(100) 

20  N=90 

21  CALL  LGAUSS(X,W.N,$1 ) 

22  1 CONTINUE 

23  WRITE! 10,2)  N 

24  2 FORMAT (14) 

25  WRITE(10,3)  (X(I),W(I),I*1 ,N) 

26  3 FORMAT ( 2E1 8 . 8) 

27  STOP 

28  END 

29  PXOT 

.0  PDATA.L  10. 

31  PENO 

32  PERS 

33  PASG.T  TP..6C9.1964  . TAPE  WITH  ARCH  ANALYSIS  PROGRAMS 

34  PMOVE  TP. 7 

35  PCOPIN  TP.. 

36  PFREE  TP. 

37  PFOR.S  TPFS.MAINARCH.M 

38  PFOR.S  TPF$. ARCHRESULTS, R 

39  PFOR.S  TPF$ . BOUNDARY ARCH , B 

40  PMAP.IN  A.B 

41  IN  TPFS.M 

42  LIB  CALCOMP*BASIC. 

43  END 

44  PASG.T  20. .6C9.0928W  . PLOT  TAPE  WITH  APPROPRIATE  NO. 

45  PXOT  B 

46  1001 

47  20.0  30.0  10.0  0.05667  0.06667  0.01 

48  90. C 0.0  0.0 

49  0.0  -0.0000663  10 

50  0.0  0.05  0.0  1.0 

51  FINITE  ELEMENT  CHECK  CASE 


PADD.P 


J. RSAREP 


©ASG.T 


TP. ,6C9. 1962 


TAPE  WITH  DATA  & REDUCTION  PROGRAM 


©MOVE  TP. 2 

FURPUR  27R2  RL72R1  11/03/77  15:26:14 

©COPIN  TP., 

27  SYM  2 REL 

©FREE  TP. 

©ASG.T  15..F14  . SCRATCH  FILE 

©ASG.T  16..F14  . FILE  WITH  EXPERIMENTAL  DATA  FOR  PLOTTING 


PDATA.IL  15. 

DATA  T7  RL70-5  11/03-15:29:22 


1 . 

4 

2. 

38.0 

2.0 

18 

3. 

C.O 

0.0472 

0.0676 

14.0 

4. 

5.0 

0.0476 

0.068 

14.0 

5. 

5.0 

0.0459 

0.C671 

14.0 

6. 

10.0 

0.0362 

0.0402 

26.0 

7. 

10.0 

0.0352 

0.03e6 

26.0 

B. 

10.0 

0.036 

0.0385 

27.0 

9. 

15.0 

0.0282 

0.0317 

39.0 

10. 

15.0 

0.0279 

0.0306 

39-0 

1 1 . 

15.0 

0.0275 

0.0303 

39.0 

12. 

20.0 

0.0241 

0.026 

49.0 

13. 

20.0 

0.0238 

0.0256 

48.0 

14. 

20.0 

0.0233 

0.0253 

48.0 

15. 

25.0 

0.0213 

0.0217 

57.0 

16. 

25.0 

0.0208 

0.0214 

58.0 

17. 

25.0 

0.0206 

0.0212 

57.0 

18. 

30.0 

0.0194 

0.0182 

73.0 

19. 

30.0 

0.0189 

0.0179 

73.0 

20. 

30.0 

0.0187 

0.0175 

73.0 

21 . 

37.5 

1 .5 

18 

22. 

5.0 

0.0714 

0.143 

23. 

5.0 

0.0634 

0.143 

24. 

5.0 

0.0684 

0.143 

25. 

10.0 

0.0574 

0.03 

26. 

10.0 

0.0538 

0.0813 

27. 

10.0 

0.0515 

0.0813 

28. 

15.0 

0.0424 

0.0654 

29. 

15.0 

0.04 

0.0621 

30. 

15.0 

0.0393 

0.0629 

31  . 

20.0 

0.0338 

0.0508 

32. 

20.0 

0.0327 

0.0488 

33. 

20.0 

0.0327 

0.0485 

34. 

2S.0 

0.0286 

0.0435 

35. 

25.0 

0.0273 

0.0413 

36. 

25.0 

0.0279 

0.0417 

37. 

30.0 

0.0273 

0.0368 

38. 

30.0 

0.0260 

0.0357 

39. 

30.0 

0.0260 

0.0355 

40. 

39.0 

3.0 

18 

41 . 

5.0 

0.0192 

0.02 

52.0 

42. 

5.0 

0.0192 

0.0192 

52.0 

43. 

5.0 

0.0192 

0.0196 

44. 

10.0 

0.0151 

0.0139 

80.0 

45. 

10.0 

0.0147 

0.0137 

80.0 

46. 

10.0 

0.0142 

0.0137 

47. 

15.0 

0.0116 

0.0109 

107.0 

48. 

15.0 

0.0121 

0.0106 

105.0 

4'-. 

15.0 

0.0121 

0.0106 

50. 

20.0 

0.0104 

0.0093 

135.0 

51 . 

20.0 

0.0102 

0.0089 

135.0 

52. 

20.0 

0.0100 

0.0091 

53. 

25.0 

0.0092 

0.0077 

170.0 

54. 

25.0 

0.0089 

0.0075 

170.0 

55. 

25.0 

0.0086 

0.0075 

30.0 

0.0080 

0.0068 

21C.0 

30.0 

0.0078 

0.0067 

205.0 

30.0 

0.0076 

0.0067 

205.0 

39.5 

3.5 

18 

5.0 

0.0119 

0.0168 

82.0 

5.0 

0.0119 

0.0168 

82.0 

5.0 

0.0116 

0.0168 

82.0 

10.0 

0.0093 

0.0118 

132.0 

10.0 

0.0089 

0.0119 

132.0 

10.0 

0.0086 

0.0118 

132.0 

15.0 

0.0075 

0.0093 

170.0 

15.0 

0.0074 

0.0089 

170.0 

15.0 

0.0072 

0.0088 

170.0 

20.0 

0.0065 

0.0079 

230.0 

20.0 

0.0064 

0.0076 

230.0 

20.0 

0.0063 

0.0075 

230.0 

25.0 

0.0057 

0.0068 

265.0 

25.0 

0.0056 

0.0065 

265.0 

25.0 

0.0056 

0.0065 

265.0 

30.0 

0.005 

0.006 

30.0 

0.0050 

0.0060 

30.0 

0.0050 

0.0061 

■ I III  !■— li— 


MFCR.S  TPFS.NONDIM.N 

FOR  00E3-1 1/03/77-15:31 :29  (1.) 


MAIN  PROGRAM 

STORAGE  USED:  C0DE(1)  000156;  OATA(O)  000135;  BLANK  COMMON(2)  000000 


EXTERNAL  REFERENCES  (BLOCK,  NAME) 


303 

NINTRS 

0004 

NRDU$ 

0005 

N 1 02$ 

0006 

NWOUS 

0007 

NWEF$ 

0010 

NSTOPS 

STORAGE  ASSIGNMENT  (BLOCK,  TYPE,  RELATIVE  LOCATION,  NAME) 


0000 

00<"  023 

IF 

0000 

000025 

1 OOF 

0000 

000102 

1 0 1 F 

0000 

00011 1 

1 02F 

0000 

000100 

1 03F 

0001 

000025 

1 13G 

0001 

000054 

1 32G 

0000 

000076 

2F 

0001 

000145 

20  1 L 

0000 

000021 

3F 

0000 

000022 

4F 

0001 

000140 

999L 

0000 

R 

000000 

C22B 

0000 

R 

000012 

flexi 

0000 

R 

00001 1 

FLEXO 

0000 

I 

000007 

ID 

0000 

I 

000002 

IS 

0000 

I 

oooeos 

JDAT  a 

0000 

1 

000001 

dSETS 

0000 

R 

000010 

P 

0000 

R 

000014 

PI 

0000 

R 

000C16 

QFLEXI 

0000 

R 

000017 

QFLEXO 

0000 

R 

000015 

QNX 

0000 

R 

000020 

ONLOAD 

0000 

R 

000003 

R 

0000 

R 

000006 

RHO 

0000 

R 

000004 

SR 

0000 

R 

000013 

wlqaq 

00101 

1* 

RE AD ( 5 , 3 ) C22B 

000000 

00105 

2* 

3 FORMAT ( F8 . 2 ) 

000007 

00106 

3* 

READ( 15.4)  dSETS 

000007 

00111 

4* 

4 FORMAT (14) 

000020 

00112 

5* 

DO  201  I S= 1 , dSETS 

000020 

00115 

6* 

RE AD( 15.1)  R.SR.JDATA 

000025 

00122 

7* 

1 FORMAT  ( 2F8 .2,14) 

000037 

00123 

B* 

RHO*R/SR 

000037 

00124 

9* 

WR I TE ( 6,100)  RHO ,C226 

000042 

00130 

10* 

100  FORMAT ( 1H1 ,//  1 0X , 1 NONDIMENSIONAL  DATA  FOR  PRESSURE  STABL1ZE0  ARCH 

000054 

00130 

11* 

1*  /// , 5X ' RADIUS  RATIO* ' Fl 0 . 4 , /2X , ' REFEEENCE  RAT  10* ' FI 0 . 4 , // 

000054 

00130 

12* 

21 5X, * NONDIMENSIONAL1 /2X , 1 PRESSURE • . 6X * FLEXI BI LI TY 1 , 6X , ' FLEXIBILITY 

000054 

00130 

13* 

3’ ,4X, ' WRINKLING  LOAD ' /2l X , • OUT ' , 18X, 1 IN' / ) 

000054 

00131 

14* 

DO  200  I D=  1 , JDATA 

000054 

00134 

15* 

READ( 15,2, END* 999 ) P. FLEXO , FLEXI .WLOAD 

000054 

00142 

16* 

PI*3. 1415 

000066 

00143 

17* 

2 FORMAT ( 4F8 . 2 ) 

000070 

00144 

18* 

<j;.X*P*SR/(2.0*C22B) 

000070 

00145 

19* 

QFLEXI*  F LEX  I *PI »C22B 

000074 

00146 

20* 

QFLEX0=FLEX0«PI*C22B 

000077 

00147 

21* 

QWL0AD*NL0A0/(PI*SR*C22B) 

000103 

00150 

22* 

NR  1 TE ( 6 . 1 0 1 ) QNX , QF LEXO , QF LEX I , QWLOAD 

0001 1 1 

00156 

23* 

WRITE! 16.103)  RHO, QNX, QFLEXO, QFLEXI, QNLOAO 

000122 

00165 

24* 

103  FORMAT ( 5E1 5 . 6 ) 

000136 

00166 

25* 

101  FORMAT( IX, El  3.6, 3X. El  3.6, 3X, El  3.6, 3X. El  3. 6) 

000136 

00167 

26* 

200  CONTINUE 

000136 

00171 

27* 

GO  TO  201 

000136 

00172 

28* 

999  WRITEI6. 102) 

000140 

00174 

29* 

102  FORMAT! IX, 1 END  OF  FILE  ENCOUNTERED  BEFORE  ALL  EXPECTED  DATA  WAS  RE 

000146 

00174 

30* 

1 AD 1 ) 

000146 

00175 

31  • 

201  CONTINUE 

000146 

00177 

32  • 

END  FILE  16 

000146 

00«00 

33* 

STOP 

000151 

00201 

34* 

END 

000155 

• 

END  OF  i 

COMPILATION:  NO 

DIAGNOSTICS. 

•XQT 

MAP20R1 

RL71-3 

11/03/77  15:32:57 

ADDRESS 

LIMITS 

001000  013735 

5598  IBANK  WORDS  DECIMAL 

040000  044202 

2179  DBANK  WORDS  DECIMAL 

STARTING  ADDRESS  012142 

1 


SEGMENT 

SMAINS 

00 

000  013735 

040000  044202 

NRWNDS/F0R-E3 

S(1) 

0010C0 

001063 

S<2) 

040000 

04001 1 

NFTCMS/F0R-E2 

S(  1 ) 

001064 

001346 

S ( 2 ) 

040012 

040025 

1 1 NBFOOS 

S(2) 

040026 

042253 

M NCNVTS/F0R68 

S(1) 

001347 

001570 

S(2) 

042254 

042350 

NFTV$/F0R~E2 

SO) 

001571 

001613 

NBDCVS/FOR-EO 

S(1) 

001614 

001744 

S(2) 

042351 

042426 

NCL0SS/F0R-E3 

SO) 

001745 

002202 

S(2) 

042427 

042454 

NSWTCS/F0R69 

SO) 

002203 

002227 

NWBLKS/F0R68 

SO) 

002230 

002341 

NBSBLS/F0R68 

SO) 

002342 

002402 

NUPDAS/F0R68 

SO) 

002403 

002436 

NRBLKS/F0R-E2 

SO) 

002437 

002461 

NOT  I N$/ F0R-E3 

S(1) 

002462 

002756 

S(2) 

042455 

042460 

NININ/FOR— E3 

SO) 

002757 

003205 

$ ( 2 ) 

042461 

042466 

N1NPTJ/F0R-E3 

SO) 

003206 

004575 

S(2) 

042467 

042522 

NFCHKS/F0R-E3 

SO) 

004576 

005567 

S(2) 

042523 

042673 

S(3) 

005570 

005570 

S ( 4) 

042674 

042745 

FORCOMS/FORFTN 

S ( 2 ) 

042746 

042753 

FORVCOMS/FOR-TE3 

S ( 2 ) 

042754 

042763 

ERUS/SYS72-8 

NERCPMS/F0R-TE3 

SO) 

005571 

005650 

$(  2) 

042764 

042777 

NI0ERS/F0R-E3 

SO) 

005651 

006070 

S ( 2 ) 

043000 

043147 

NFMTS/F0R-E3 

SO) 

006071 

006753 

S<2) 

043150 

043224 

N0UTS/F0R-E3 

SO) 

006754 

010470 

5(2) 

043225 

043266 

NTABS/FOR— TE3 

$(2) 

043267 

043325 

NSTO' S/F0R-TE3 

SO) 

010471 

010533 

$(2) 

043326 

043335 

NWEFS/F0R-E2 

SO) 

010534 

010741 

S(  2) 

043336 

C433S5 

N0BUFS/F0R68 

so ) 

010742 

011002 

NIBUFS/F0R-E2 

SO) 

011003 

011042 

$(2) 

043356 

C43356 

NINTRS/F0R-E3 

SO) 

011043 

011115 

S(  2) 

043357 

043374 

SQRTS/F0RS9 

SO) 

011116 

011156 

S(2) 

043375 

043406 

NERRS/F0R-E3 

SO) 

011157 

011520 

S(2) 

043407 

043566 

NIERS/F0R-E3  5(1)  011521  011676 
NOSYMS/ F0R-E3  S(1)  011677  012141 
BLANKSCOMMON ( COMMGNBLOCK) 

N $(1)  012142  012317 
FPMIN  SI1)  012320  013545 
WRITEMATRIX  5(1)  013546  013735 


SYS$*RLI8$.  LEVEL  72-8 
END  MAP 


S(2)  043567  043706 

5(2)  043707  043710 

5(0)  043711  044045 

5(2)  BLANKSCOMMON 

5(0)  044046  044127 

5(2)  BLANKSCOMMON 

5(0)  044130  044202 

5(2)  BLANKSCOMMON 


NGNDIMENSIONAL  DATA  FOR  PRESSURE  STABLIZED  ARCH 


RADIUS  RATIO*  T9. OOOO 
REFEEENCE  RATIO*  313.0000 


PRESSURE 


.159744-01 

.159744-01 

.159744-01 

.319469-01 

.319489-01 

.319489-01 

.479233-01 

.479233-01 

.479233-01 

.638978-01 

.638978-01 

.638978-01 

.798722-01 

.798722-01 

.798722-01 

.958466-01 

.958466-01 

.958466-01 


NONDIMENSIONAl 

FLEXIBILITY 

OUT 

.4641 1 3+02 
.468046+02 
.451330+02 
.355951+02 
.346118+02 
.353984+02 
.277288+02 
.274338+02 
.270405+02 
.236973+02 
.234023+02 
.229106+02 
.209441+02 
.204524+02 
.2025 j8+02 
.190758+02 
.185842+02 
.183875+02 


FLEXIBILITY 

IN 

.664704+02 

.668637+02 

.659787+02 

.395282+02 

.379550+02 

.378566+02 

.311703+02 

.300887+02 

.297937+02 

.255655+02 

.251722+02 

.248772+02 

.213374+02 

.210424+02 

.208457+02 

.178959+02 

.176009+02 

.172076+02 


WRINKLING  LOAO 


.711896-02 

.711896-02 

.711896-02 

.132209-01 

.132209-01 

.137294-01 

.198314-01 

.198314-01 

.198314-01 

.249164-01 

.244079-01 

.244079-01 

.289843-01 

.294928-01 

.289843-01 

.371203-01 

.371203-01 

.371203-01 


NGNDIMENSIONAL  DATA  FOR  PRESSURE  STABLIZED  ARCH 


RADIUS  RATIO*  25.0000 


REFEEENCE  RATIO*  313.000 

NOND I ME  NS  I ON A L 

PRESSURE 

flexibility 

out 

.119808-01 

.702069+02 

.119808-01 

.672570+02 

.119808-01 

.672570+02 

.239617-01 

.564408+02 

.239617-01 

.529010+02 

.239617-01 

.506394+02 

.359425-01 

.416915+02 

.359425-01 

.393316+02 

.369425-01 

.386433+02 

.479233-01 

.332352+02 

.479233-01 

.321536+02 

.479233-01 

.321536+02 

.599042-01 

.281221+02 

.599042-01 

.268438+02 

.599042-01 

.274338+02 

.718850-01 

.268438+02 

.718850-01 

.255655+02 

.718850-01 

.255655+02 

flexibility 

IN 

WRINKLING  LOAD 

.140610+03 

.000000 

. 140610+03 

.000000 

.140610+03 

.000000 

.766632+02 

.000000 

.799414+02 

.000000 

. 799414+02 

.000000 

.643071+02 

.000000 

.610623+02 

.000000 

.618489+02 

.000000 

.499511+02 

.000000 

.479845+02 

.000000 

.476895+02 

.000000 

.427731+02 

.000000 

. -'.06099+02 

.000000 

.410032+02 

.000000 

.361851+02 

.000000 

.351034+02 

.000000 

.349068+02 

.000000 

NONOIMENSIONAL  DATA  FOR  PRESSURE  STABllZEO  ARCH 


RAOlUS  RATIO* 
REFEEENCE  RATIO* 


13.0000 

313.0000 


NONOIMENSIONAL 

PRESSURE  FLEXIBILITY  FLEXIBILITY  WRINKLING  LOAD 

OUT  IN 


.239617-01 
.239617-01 
.239617-01 
.479233-01 
.479233-01 
.479233-01 
.718850-01 
.718850-01 
.718850-01 
.958466-01 
.958466-01 
.958466-01 
.1 1 9808*00 
.11 9808+00 
.11 9808+00 
. 1 43770+00 
.143770+00 
.143770+00 


.188792+02 
.188792+02 
.188792+02 
.148477+02 
.144544+02 
.139627+02 
.114062+02 
. 118978+02 
.118978+02 
.102262+02 
.100296+02 
.983289+01 
.904626+01 
.875128+01 
.845629+01 
.786632+01 
.766966+01 
.747300+01 


.196658+02 

.188792+02 

.192725+02 

.136677+02 

.134711+02 

.134711+02 

.107179+02 

.104229+02 

.104229+02 

.914459+01 

.875128+01 

.894793+01 

.757133+01 

.737467+01 

.737467+01 

.668637+01 

.656304+01 

.658804+01 


.176279-01 
. 1 '<6279-01 
.000000 
.271199-01 
.271199-01 
.000000 
.362728-01 
.355948-01 
.000000 
.457648-01 
.457648-01 
.000000 
.576297-01 
.576297-01 
.000000 
.711896-01 
.694946-01 
.694946-01 


NONDIMENSIONAL  data  for  pressure  stablized  arch 

RADIUS  RATIOS  11.2857 
REFEEENCE  RATIO=  313.0000 


PRESSURE 


.279553-01 
.279553-01 
.279553-01 
.559105-01 
.559105-01 
.559105-01 
.838658-01 
.838658-01 
.838658-01 
.111821+00 
. 1 1 1821+00 
.111821+00 
.139776+00 
.139776+00 
.139776+00 
.167732+00 
.167732+00 
.167732+00 


NONDIMENSIONAL 

FLEXIBILITY 

OUT 

.117011+02 
. 117011+02 
.114062+02 
.914459+01 
.875128+01 
.845629+01 
.737467+01 
.727634+01 
.707969+01 
.639130+01 
.629305+01 
.619472+01 
.560475+01 
.550642+01 
.5506  12+01 
.491645+01 
.491645+01 
.491645+01 


FLEXIBILITY 

IN 

.165193+02 
.165193+02 
.165193+02 
. 116028+02 
. 117011+02 
. 116028+02 
.914459+01 
.875128+01 
.865295+01 
.776799+01 
.747300+01 
.737467+01 
.668637+01 
.639138+01 
.639138+01 
.589974+01 
.589974+01 
.599807+01 


WRINKLING  LOAD 


.238267-01 
.239267-01 
.238267-01 
.383552-01 
.383552-01 
.383552-01 
.493969-01 
.493969-01 
.493969-01 
.668311-01 
.66831 1-01 
.668311-01 
.770010-01 
.770010-01 
.770010-01 
.000000 
.000000 
.000000 


•DATA, l 16. 


DATA  T7 

RL70-5  11/03-15:34:52 

1 . 

.190000+02 

.159744-01 

2. 

. 190000+02 

.159744-01 

3. 

.190000+02 

.159744-01 

4. 

. 190000+02 

.319489-01 

5. 

. 190000* ^2 

.319489-01 

6. 

.190000+02 

.319489-01 

7. 

.190000+02 

.479233-01 

8. 

.190000+02 

.479233-01 

9. 

.190000+02 

.479233-01 

10. 

.190000+02 

.638978-01 

1 1 . 

.190000+02 

.638978-01 

12. 

.190000+02 

.638978-01 

13. 

.190000+02 

.798722-01 

14. 

.190000+02 

.798722-01 

15. 

.190000+02 

.798722-01 

16. 

.190000+02 

.958466-01 

17. 

. 190000+02 

.958466-01 

18. 

.190000+02 

.958466-01 

19. 

.250000+02 

. 1 19808-01 

20. 

.250000+02 

. 1 19808-01 

21  . 

.250000+02 

. 119808-01 

22. 

.250000+02 

.239617-01 

23. 

.250000+02 

.239617-01 

24. 

. 250000+02 

-239C57-01 

25. 

.250000+02 

.359425-01 

26. 

.250000+02 

.359425-01 

27. 

.250000+02 

.359425-01 

28. 

.250000+02 

.479233-01 

29. 

.250000+02 

.479233-01 

30. 

.250000+02 

.479233-01 

31. 

.250000+02 

.5990-2-01 

32. 

.250000+02 

.599042-01 

33. 

.250000+02 

.599042-01 

34. 

.250000+02 

.718850-01 

35. 

.250000+02 

.718850-01 

36. 

.250000+02 

.718850-01 

37. 

. 130000+02 

.239617-01 

38. 

.130000+02 

.239617-01 

39. 

.130000+02 

.239617-01 

40. 

.130000+02 

.479233-01 

41  . 

.130000+02 

.479233-01 

42. 

.130000+02 

.479233-01 

43. 

.130000+02 

.718850-01 

44. 

.130000+02 

.718850-01 

45. 

.130000+02 

.718850-01 

46. 

.130000+02 

.958466-01 

47. 

.130000+02 

.958466-01 

48. 

.130000+02 

.958466-01 

4'.. 

.130000+02 

. 119808+00 

50. 

.130000+02 

. 119808+00 

51  . 

.130000+02 

. 119808+00 

52. 

.130000+02 

.143770+00 

53. 

.130000+02 

.143770+00 

54. 

.130000+02 

.143770+00 

55. 

.112857+02 

.279553-01 

.4641 13*02 
.468046+02 
.451330+02 
.365951+02 
.3461 18+02 
.353984+02 
.277288+02 
.274338+02 
.270405+02 
.236973+02 
.234023+02 
.229106+02 
.209441+02 
.204524+02 
.202558+02 
.190758+02 
.185842+02 
.183875+02 
.702069+02 
.672570+02 
.672570+02 
.564408+02 
.529010+02 
.506394+02 
.416915+02 
.393316+02 
.386433+02 
.332352+02 
.321536+' ’ 
.321536+02 
.281221+02 
.268438+02 
.274338+02 
.268438+02 
.255655+02 
.255655+02 
.188792+02 
.188792+02 
.188792+02 
.148477+02 
.144544+02 
.139627+02 
.114062+02 
. 118978+02 
.118978+02 
.102262+02 
.100256+02 
.983289+01 
.904626+01 
.875128+01 
.845629+01 
.786632+01 
.766966+01 
.747300+01 
. 117011+02 


1 


.664704+02 
.666637+02 
.659787+02 
.395262+02 
.379550+02 
.376566+02 
.311703+02 
.300887+02 
.297937+02 
•255655+02 
.251722+02 
.248772+02 
.213374+02 
.210424+02 
.208457+02 
.178959+02 
.176009+02 
.172076+02 
.140610+03 
.140610+03 
.140610+03 
.786632+02 
.799414+02 
.799414+02 
.643071+02 
.610623+02 
.618489+02 
.49951 1+02 
.479845+02 
.476895+02 
.427731+02 
.406099+02 
.410032+02 
.361851+02 
.351034+02 
.349068+02 
.196658+02 
.188792+02 
.192725+02 
.136677+02 
.134711+02 
.134711+02 
.107179+02 
.104229+02 
.104229+02 
.914459+01 
.875128+01 
.894793+01 
.757133+01 
.737467+01 
.737467+01 
.668637+01 
.658804+01 
.658804+01 
.165193+02 


.711896-02 
.711696-02 
.711896-02 
.132209-01 
.132209-01 
.137294-01 
.198314-01 
.198314-01 
.198314-01 
.249164-01 
.244079-01 
.244079-01 
.289843-01 
.294928-01 
.289843-01 
.371203-01 
.371203-01 
.371203-01 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.000000 
.176279-0 1 
.176279-01 
.000000 
.271199-01 
.271199-01 
.000000 
.362728-01 
.355948-01 
.000000 
.457648-01 
.457648-01 
.000000 
.576297-01 
.576297-01 
.000000 
.711896-01 
.694946-01 
.694946-01 
.238267-01 


56. 

. 112857+02 

.279553-01 

.117011+02 

.165193+02 

57. 

. 112857+02 

.279553-01 

.114062+02 

.165193+02 

58. 

.112857+02 

.559105-01 

.914459+01 

.1  16028  + 02 

59. 

. 112857+02 

.559105-01 

.875128+01 

. 1 17011+02 

60. 

.112857+02 

.559105-01 

.845629+01 

.1 16028+02 

61 . 

. 112857+02 

.836653-01 

.737467+01 

.914459+01 

62. 

.112857+02 

.838658-01 

.727634+01 

.875128+01 

63. 

. 112857+02 

.838658-01 

.707968+01 

.865295+01 

64. 

.112857+02 

.111821+00 

.639138+01 

.776799+01 

65. 

. 112857+02 

.11 1821+00 

.629305+01 

.747300+01 

66. 

. 112857+02 

.111821+00 

.619472+.’. 

.737467+01 

67. 

.112857+02 

.139776+00 

.560475+01 

.668637+01 

68. 

.112857+02 

.139776+00 

.550642+01 

.639138+01 

69. 

. 112857+02 

.139776+00 

.550642+01 

.6391 38+01 

70. 

.112857+02 

.167732+00 

.491645+01 

.589974+01 

71  . 

.112857+02 

.167732+00 

.491645+01 

.589974+01 

72. 

ENO  DATA. 

.112857+02 

.167732+00 

.491645+01 

.599807+01 

VERS 

FURPUR  27R2  RL72R1  11/03/77  15:34:56 
ENO  ERS. 


VERS  15. 

^ ENO  ERS. 

<0 

VASG.T  10..F14  . FILE  WITH  GAUSSIAN  INTEGRATION  CONSTANTS 


.238267-01 

.238267-01 

.383552-01 

.383552-01 

.383552-01 

.493969-01 

.493969-01 

.493969-01 

.66831 1-01 

.66831 1-01 

.668311-01 

.770010-01 

.770010-01 

.770010-01 

.000000 

.000000 

.000000 


•FOR, IS 

FOR  00E3-1 1/03/77-15:42:08 


MAIN  PROGRAM 

STORAGE  USED:  CODE(1)  000041;  DATA(O)  000321;  BLANK  C0MM0N(2)  C00000 
EXTERNAL  REFERENCES  (BLOCK,  NAME) 


• J03 

LGAUSS 

0004 

NINTHS 

0005 

NWDUS 

0006 

N 1 02$ 

0007 

NI01S 

0010 

NSTOPS 

STORAGE  ASSIGNMENT  (BLOCK,  TYPE,  RELATIVE  LOCATION,  NAME) 

0001  000012  1 L 0001  000030  114G  0000  000312  2F  0000  000313  3F  0000  I 000311  I 

0000  I 000310  N 0000  R 000144  W 0000  R 000000  X 


00101 

1* 

DIMENSION  X(100).W(100) 

000000 

00103 

2* 

N=90 

000001 

00104 

3* 

CALL  LGAUSS(X,W,N.$1 ) 

000003 

00105 

4» 

1 

CONTINUE 

000012 

00106 

5* 

WRI TE( 10,2)  N 

000012 

001 1 1 

6* 

2 

FORMAT (14) 

000021 

00112 

7* 

WRITE( 10,3)  (X(I).W(I),I«1 ,N) 

000021 

00121 

8* 

3 

FORM AT (2E1 0.8) 

000034 

00122 

9» 

STOP 

000034 

00123 

10* 

END 

000040 

END  OF  COMPILATION:  NO  DIAGNOSTICS. 


•XQT 

MAP2BR1  RL71-3  11/03/77  15:43:13 


ADDRESS  LIMITS  001000  010071  3642  IBANK  WORDS  DECIMAL 

040000  044220  2193  DBANK  WORDS  DECIMAL 

STARTING  ADORESS  010031 


SEGMENT  SMAINS  001000  010071  040000  044220 

NSWTCS/F0R69  S ( 1 ) 001000  001024 

NRBLKS/F0R-E2  $(1)  001025  001047 

NRWNDS/F0R-E3  S(1)  001050  001133  S(2)  040000  040011 


NHEFS/F0R-E2 

NBDCVS/ FOR-EO 

NFTVS/F0R-E2 

NCNVTS/F0R68 

NCLOSS/FOR-EU 

NWBLKS/FOR68 

NBSBLS/F0R68 

NUPDAS/F0R68 

NBFOOS 

NOT IN$/ FOR-E3 
N0UT$/F0R-E3 
NFMT  5/F0R-E3 
NI0ERS/F0R-E3 
NFCHKS/F0R-E3 

NTABS/FOR— TE3 

FORCOMS/FORFTN 

ERUS/SYS72-8 

NERCOMS/FOR-TE3 

F0RVC0MS/F0R-TE3 

NERRS/FOR-E3 

NSTOPS/FOR-TE3 

NIER5/FOR-E3 

N0BUFS/F0R68 

NINTRS/FOR-E3 

L GAUSS 


Sd) 
*(  1 ) 
$(1) 
SO) 
Sio 
SO) 
SO) 
SO) 

SO) 
SO) 
SO) 
SO) 
SO) 
S(  3) 


SO) 

SO) 

SO) 

SO) 

SO) 

SO) 


BLANKSCOMMON(COMMONBLOCK) 
NAMES  SO) 


001134  001341 
001342  001472 
001473  001515 
001516  001737 
001740  002175 
002176  002307 
002310  002350 
002351  002404 

002405  002701 
002702  004416 
004417  005301 
005302  005521 
005522  006513 
006514  006514 


006515  006574 

006575  007136 
007137  007201 
007202  007357 
007360  007420 
007421  007473 
007474  010030 


010031  010071 


SYSS«RLI8S.  LEVEL  72-8 
END  MAP 


S(2)  040012  040031 

$(2)  040032  040107 

$(2)  0401 10  040204 

S(2)  040205  040232 


S(2)  040233  042460 
S(2)  042461  042464 
S(2)  042465  042526 
S ( 2 ) 042527  042603 
$(2)  042604  042753 
S ( 2 ) 042754  043124 
S ( 4 ) 043125  043176 
S( 2 ) 043177  043235 
S(2)  043236  043243 

S(2)  043244  043257 
S (2)  043260  043267 
S(2)  043270  043447 
$(2)  043450  043457 
S(2)  043460  043577 

S(2)  043600  043615 
S(0)  043616  043677 

$(2)  blankscommon 


S(0)  043700  044220 

$(2)  BLANKSCOMMON 


PDA  T A , L 10. 


OATA  T7 

RL70-5  11/03-15:46:06 

1 . 

90 

2. 

-.99964698+00 

.90525075-03 

3. 

-.99814039+00 

.21074554-02 

4. 

-.99543183+00 

.33086729-02 

5. 

-.99152394+00 

. 45059G90-02 

6. 

-.98642138+00 

.56979301-02 

7. 

-.98013026+00 

.68829324-02 

8. 

-.97265817+00 

.80596809-02 

9. 

-.96401411+00 

.92266372-02 

10. 

-.95420849+00 

.10382656-01 

11  . 

-.9432531 1+00 

. 1 1525941-C1 

12. 

-.93116120+00 

. 12655360-01 

13. 

-.91794731+00 

.13769669-01 

14. 

-.90362736+00 

. 14867329-01 

15. 

-.88821861+00 

.15947057-01 

16. 

-.87173963+00 

.17007603-01 

17. 

-.85421026+00 

.18047633-01 

18. 

-.83565165+00 

.19065906-01 

19. 

-.81608614+00 

.20061191-01 

20. 

-.79553731+00 

.21032342-01 

21  . 

-.77402992+00 

.21973135-01 

22. 

-.75158988+00 

.22897445-01 

23. 

-.72824423+00 

.23789182-01 

24. 

-.704021 1 1+00 

.24652231-01 

25. 

-.67894970+00 

.25485585-01 

26. 

-.65306021+00 

.26288226-01 

27. 

-.62638382+00 

. 27C592 1 3-0 1 

28. 

-.59895270+00 

.27797571-01 

2 •• 

-.57079988+00 

.28502445-01 

30. 

-.54195929+00 

.29172971-01 

31  . 

-.51246569+00 

.29808359-01 

32. 

-.48235460+00 

.30407821-01 

33. 

-.45166232+00 

.30970-540-01 

34. 

-.42042581+00 

.31496140-01 

35. 

-.38868273+00 

.31993702-01 

36. 

-.35647131+00 

.32432720-01 

37. 

-.32383037+00 

.32842656-01 

38. 

-.29079925+00 

.33213020-01 

39. 

-.25741773+00 

.33543365-01 

40. 

-.22372605+00 

.33833297-01 

41  . 

-.18976479+00 

.34082457-01 

42. 

-.15557488+00 

.34290551-01 

43. 

-.12119751+00 

.34457333-01 

44. 

-.866741 12-01 

.34582590-01 

45. 

— . 520462  77-01 

.34666181-01 

46. 

-.17355731-01 

.34707999-01 

47. 

.17355731-01 

.34707999-01 

48. 

.52046277-01 

.34666181-01 

49. 

.86674112-01 

.34582590-01 

50. 

.121 19751+00 

.34457333-01 

51. 

.15557^38+00 

.34290551-01 

02. 

.18976479+00 

.34082457-01 

53. 

.22372605+00 

.33833297-01 

54. 

.25741773+00 

.33543365-01 

55. 

.29079925+00 

.33213020-01 

56. 

.32363037+00 

.32842656-01 

57. 

.35647131+00 

.32432720-01 

58. 

.38868273+00 

.31983702-01 

59. 

.42042581+00 

.31495140-01 

60. 

.45165232+00 

.30970640-01 

61  . 

.48235460+00 

.30407021-01 

62. 

.51246569+00 

.29808359-01 

63. 

.54195929+00 

.29172971-01 

64. 

.57079988+00 

.28502445-01 

65. 

.59895270+00 

.27797571-01 

6 •. 

.62638382+00 

.27059213-01 

67. 

.65306021+00 

. 26286226-01 

68. 

.67894970+00 

.25485585-01 

69. 

.704021 1 1+00 

.24652231-01 

70. 

.72824423+00 

.23789182-01 

71  . 

.75158988+00 

.22897445-01 

72. 

.77402992+00 

.21978135-01 

73. 

.79553731+00 

.21032342-01 

74. 

.81608614+00 

.20061191-01 

75. 

.83565165+00 

. 19065906-01 

76. 

.85421026+00 

.18047633-01 

77. 

.87173963+00 

. 17007603-01 

78 . 

.88821861+00 

.15947057-01 

79. 

.90362736+00 

.14867329-01 

80. 

.91794731+00 

.13769669-01 

81 . 

.93116120+00 

.12655380-01 

82. 

.94325311+00 

.11525941-01 

83. 

.95420849+00 

.10382556-01 

84. 

.96401411+00 

.92266372-02 

85. 

.97265817+00 

.80596809-02 

86. 

.98013026+00 

.68829324-02 

87. 

.98642138+00 

.56979301-02 

88. 

.99152394+00 

.45059690-02 

89. 

.99543183+00 

.33086729-02 

90. 

.99814039+00 

.21074554-02 

91 . 

.99964698+00 

.90525975-03 

END  DATA. 

•ERS 

FURPUR  27R2  RL72R1  11/03/77  15:46:15 
END  ERS. 


•ASG.T  TP..6C9.1964  . TAPE  WITH  ARCH  ANALYSIS  PROGRAMS 


CMOVE  TP  7 

FURPUR  27R2 ' RL72R1  11/03/77  15:46:37 


•COPIN  TP., 
3 SYM 


• FREE 


TP 


•FOR.S  TPFS.MAINARCH.M 

FOR  00E3-1 1/03/77-15:57:21  (0.) 


MAIN  PROGRAM 

STORAGE  USED:  CODE(I)  003227;  OATA(O)  004616;  BLANK  COMMON(2)  001S41 


EXTERNAL  REFERENCES  (BLOCK.  NAME) 

0003  PLOTS 

0004  NEWPEN 

0005  OGJR 

0006  CALCA 

0007  CALPAR 

0010  SCALE 

0011  AXIS 

0012  SYMBOL 

0013  NUMBER 

0014  PLOT 
00.5  FLINE 

0016  LINE 

0017  EXIT 

0020  NINTRS 

0021  NRDCS 

4k  0022  N 102$ 

0023  NRDUS 

0024  NL03S 

0025  DSORT 

0026  NPRT$ 

0027  DSIN 

0030  DCOS 

0031  NREWS 

0032  NI01 $ 

0033  NWEFS 

0034  NSTOPS 


STORAGE  ASSIGNMENT  (BLOCK,  TYPE,  RELATIVE  LOCATION.  NAME) 


0001 

000027 

10L 

0000 

004274 

1001  F 

0000 

004300 

1002F 

0000 

004335 

1005F 

0000 

004346 

1006F 

0000 

004220 

1009F 

0000 

004204 

1012F 

0000 

004373 

1013F 

0000 

004430 

1 0 1 4F 

0000 

004464 

1017F 

0000 

004471 

1020F 

0000 

004175 

1 021  F 

0000 

004153 

1041  F 

0000 

004147 

1042F 

0001 

002062 

1 10SG 

0001 

000125 

20L 

0000 

004473 

2001  F 

0001 

001076 

250L 

0001 

001054 

252L 

0001 

000415 

255G 

0001 

001036 

300L 

0000 

004500 

3003F 

0000 

004502 

3004F 

0000 

004172 

3005' 

0001 

000535 

322G 

0001 

000537 

325G 

0001 

001266 

328L 

0001 

000576 

346G 

0001 

000605 

353G 

0001 

000625 

365G 

r.001 

C01676 

377L 

0001 

001733 

395L 

0001 

001734 

396L 

0001 

002066 

401  L 

0001 

000651 

404G 

0001 

000165 

45L 

0001 

002067 

500L 

0001 

001040 

507G 

0001 

002200 

5 1 0 L 

0001 

002365 

525L 

0001 

002375 

530L 

0001 

002415 

540  L 

0001 

002660 

5B0L 

0001 

002740 

600L 

0001 

001204 

605G 

0001 

001223 

62 1 G 

0001 

001224 

624G 

0001 

001244 

635G 

0001 

001273 

656G 

0001 

001313 

670G 

0001 

001324 

674G 

0000 

004316 

1003F 

0000 

004324 

1004F 

0000 

004216 

1 01  OF 

0000 

004352 

1 0 1 1 F 

0000 

004450 

1 01 5F 

0000 

004457 

1016F 

0000 

004200 

1022F 

0000 

004160 

1040F 

0001 

002104 

1 122G 

0001 

002112 

1 132G 

0001 

000411 

25 1 G 

0001 

001200 

25 1 L 

0000 

004475 

300 1 F 

0000 

004476 

3002F 

0000 

004173 

3006F 

0001 

001364 

306L 

0001 

001335 

329L 

0001 

000542 

331G 

000' 

000626 

370G 

0001 

001653 

376L 

0001 

000144 

40L 

0001 

000650 

40 1 G 

0001 

000755 

460G 

0001 

000177 

SOL 

0001 

002334 

5 1 BL 

0001 

002342 

520L 

0001 

002417 

550L 

0001 

002464 

570L 

0001 

001205 

61  OG 

0001 

002743 

6 1 OL 

0001 

001245 

640G 

0001 

003150 

6S0L 

0001 

001353 

710G 

0001 

001576 

763G 

r 


1 


0001 

001577 

766G 

0001 

001610 

776G 

0001 

003166 

950L 

0001 

003161 

960L 

0001 

003205 

990L 

0000 

0 

000033 

A 

0000 

0 

004050 

AC 

0000 

0 

004054 

AE 

0000 

D 

004035 

ALPHA 

0000 

R 

003307 

AN 

0000 

R 

003306 

AR 

0003 

0 

000007 

ARR 

0000 

R 

003310 

AW  I 

0000 

R 

00331 1 

AMO 

0000 

R 

00331 2 

AWRL 

0000 

D 

000473 

B 

0000 

D 

004053 

BC 

0000 

0 

004C55 

BE 

0002 

0 

000053 

BETA 

0002 

D 

000175 

BLAM2 

0003 

D 

00031 1 

BLAM3 

0000 

0 

004063 

C 

0000 

R 

001760 

CN 

0002 

0 

00011 1 

CSA 1 

0002 

D 

000115 

CSA3 

0000 

R 

001613 

CM 

0000 

R 

002136 

CWL 

0003 

D 

000161 

0 

0002 

D 

000225 

DA 

0007 

0 

000005 

E 

0000 

R 

003140 

ELN 

0000 

R 

002634 

LN 

0000 

R 

002310 

EMI 

0000 

R 

002772 

EWL 

0000 

R 

002456 

EwO 

0000 

0 

004043 

F 

0003 

0 

000075 

G 

0002 

D 

000 1 25 

H 

0002 

0 

000121 

HCA2 

0002 

D 

000123 

HSA2 

0000 

I 

004075 

1 

0003 

1 

001537 

160UN 

0002 

I 

000004 

IC 

0000 

I 

004102 

ICL 

0000 

I 

004034 

IDUM 

0000 

I 

004146 

1 EOF 

0003 

I 

001540 

IdUNC 

0000 

1 

004033 

I LOAD 

0002 

I 

001535 

INTEG 

0000 

I 

004144 

IP 

0000 

I 

004033 

I PLOT 

0000 

I 

004141 

ISAV 

0002 

I 

001536 

1TVP 

0000 

I 

004077 

J 

0000 

I 

001576 

JC 

0000 

1 

004143 

JP 

0000 

I 

004076 

K 

0000 

1 

004106 

KSET 

0000 

I 

004103 

LIN 

0000 

I 

004145 

NEXP 

0000 

1 

004060 

NIS 

0000 

I 

004047 

UMAX 

0000 

I 

004142 

NVAL 

0000 

D 

004107 

PA 

0000 

0 

004111 

PB 

0000 

D 

004133 

QCB 

0000 

0 

004137 

QCH 

0000 

D 

004131 

QG1 

0000 

D 

004133 

0G2 

0000 

D 

004135 

0G3 

0000 

0 

C04137 

QG4 

0000 

D 

000000 

01  DEN 

0009 

D 

004121 

QSB 

0000 

0 

004125 

OSH 

0002 

0 

000000 

RHO 

0000 

D 

004067 

RHOSQ 

0000 

D 

004071 

RHCSQI 

0000 

0 

003314 

RlIM 

0000 

R 

003313 

RRHO 

0000 

D 

004045 

SVG 

0003 

0 

000113 

SNAJ 

0003 

D 

000117 

SNA3 

0002 

0 

000715 

SOL 

0000 

0 

001132 

T 

0000 

D 

003342 

TARR 

0000 

0 

004104 

TEMPA 

0000 

D 

004115 

THEHAR 

0000 

0 

004041 

THEHAT 

0000 

0 

004037 

THETA 

0000 

D 

004117 

thetar 

0000 

0 

004113 

THETR 

0000 

D 

001573 

V 

0000 

D 

004100 

WF  A 

0000 

R 

002304 

WLSCAL 

0000 

R 

002274 

MSCAL 

0000 

D 

003330 

X 

0000 

0 

004061 

xc 

0000 

D 

004065 

XE 

0002 

0 

000145 

XI 

0002 

D 

000002 

XN 

0000 

D 

003333 

xnlim 

0000 

R 

003300 

xnscal 

0002 

0 

000153 

XR 

0002 

C 

000067 

V LAM 

0000 

0 

004073 

YLAMSQ 

0000  0 003334  Z 


N> 

CJI 


00100 

1* 

C 

000000 

00101 

2* 

IMPLICIT  DOUBLE  PRECISION(A-H.O-Z) 

000000 

00101 

3* 

c 

000000 

00101 

4* 

c 

ELASTIC  ANALYSIS  OF  PRESSURE  STABILIZED  ARCH 

000000 

0010 

5* 

c 

000000 

00101 

6* 

c 

REOUESTED  BY  EARL  C STEEVES,  GEPL  X2406 

000000 

00101 

7* 

c 

000000 

00101 

8* 

c 

PURPOSE  - SOLUTION  OF  SET  OF  SECOND  ORDER  DIFFERENTIAL  EOUATIONS 

000000 

00101 

9* 

c 

WITH  CONSTANT  COEFFICIENTS 

000000 

00101 

10* 

c 

000000 

00101 

11* 

c 

C CHRETIEN 

000000 

00101 

12* 

c 

OCTOBER  1973 

000000 

00101 

13* 

c 

CHANGE  FEB  1974 

000000 

00101 

14* 

c 

000000 

00101 

15* 

c 

000000 

00101 

16* 

c 

000000 

00101 

17* 

c 

H=  ARRAY  CONTAINING  CONSTANT  COEFFICIENTS  IN  EOUATIONS 

000000 

00101 

18* 

c 

E.RHO.C.XN  - INPUT  POSSIBLY  VARYING  OVER  AN  INTERVAL 

000000 

00101 

19* 

c 

ARRAYS  EL1M. RLIM, CLIM, XNLIM. ARE  USED  TO  READ  INI T , F INAL, DELTA 

000000 

00101 

20* 

c 

OF  EACH 

000000 

00101 

21* 

c 

X INT =ARRAY  OF  3 VARIABLES  DEFINED  IN  TERMS  OF  INTEGRAL 

000000 

00101 

22* 

c 

NUMERICALLY  EVALUATED 

000000 

00103 

23* 

COMMON  RHO, XN , IC . E 

000001 

00104 

24* 

COMMON  ARR (3,6) , BET A( 6 ) , YLAM( 3 ) ,G(6) , CSA1 ,SNA1 .CSA3.SNA3 

000001 

00105 

25* 

COMMON  HCA2.HSA2.Hl 8) 

000001 

00106 

26* 

COMMON  XI (3) ,XR(3) .0(6) , BLAM2 ( 6 ) , BLAM3 ( 6 ) 

000001 

00107 

27* 

COMMON  DA( 12. 13) ,SOL(200) , INTEG. I TYP, I BOUN , I JUNC 

000001 

00110 

28* 

DIMENSION  Q I DEN ( 13) 

000001 

00111 

29* 

DIMENSION  A(12.12).B(12.12),T(12,12).V(2),JC(12) 

000001 

00112 

30* 

REAL  CM (102). CN (102), CWL( 102) 

000001 

00113 

31* 

REAL  WSCAL(4) , XNSCAL ( 4 ) , WLSCAL ( 4 ) 

000001 

00114 

32* 

REAL  EWI( 102) ,EW0l 102) ,EN( 102) ,EWL(102) ,ELN(102) 

000001 

00115 

33* 

REAL  AR.AN, AWI ,AWO. AWRL.RRHO 

000001 

00116 

34* 

DIMENSION  RL1M(3) , XNL IM( 3 ) , X ( 2 ) , Z( 3 ) 

000001 

00117 

35  • 

DIMENSION  TARR (12,13) 

000001 

00120 

36* 

READ  3001  I PLOT , I LOAD . I NT EG , I BOUN 

000001 

00120 

37* 

C 

IBOUN=0  CLAMPED  ENDS 

000001 

00120 

38* 

C 

I BOUN* 1 SIMPLY  SUPPORTED  ENDS 

000001 

00120 

39* 

C 

INTEG=0  IMPLIES  NO  INTEGRATION  OF  GREEN  S FUNCTION 

000001 

00120 

40* 

c 

iNTEGs 1 INPLIES  INTEGRATION  OF  GREEN  S FUNCTION 

000001 

00120 

41* 

c 

ITYP=0  CONCENTRATED  LOAO  AND  INTEGRATION  OF  GREEN  FUNCTION  000001 

0012 

42* 

c 

ITYP*  1 UNIFORM  LOAD 

000001 

00120 

43* 

c 

I L0AD=0  IMPLIES  NORMAL  LOAD 

000001 

00120 

44* 

c 

I LOAD: 1 IMPLIES  VERTICAL  LOAD 

000001 

00126 

45* 

IF  ( IPLOT  .EQ.  1 ) GO  TO  10 

000011 

00130 

46* 

I DUM= 1 00 

000014 

00131 

47* 

CALL  PLOTS  ( IDUM, IDUM.20) 

000016 

00132 

48* 

CALL  NEWPEN  ( 1 ) 

000023 

00132 

49* 

c 

CALL  PLOT  (0.0, -36. 0,-3) 

000023 

00132 

50* 

c 

CALL  PLOT  (0.0,2. 0,-3) 

000023 

00133 

51  * 

10  READ!  5,3002.END=990)  RLIM.XNLIM 

000027 

00137 

52* 

READ  1001 .ALPHA, THETA,  THEHAT 

000042 

00144 

53* 

READ  2001,  F.SMG.NMAX 

000051 

00151 

54* 

READ  3002  AC , BC , AE, BE 

000060 

00157 

55* 

NISsO 

000070 

00160 

56* 

READ  1010 , QIDEN 

000071 

00163 

57* 

X(1)=0.0 

000100 

00164 

58* 

X(2)=2. 0*3. 14159 

000102 

00165 

59* 

RH0=RLIM( 1 ) 

000104 

00166 

60* 

XC=AC*XNLIM(1 )+BC 

000106 

00167 

61* 

C=XC 

000112 

00170 

62* 

XE*AE*XNLIM(1 )*BE 

000113 

00171 

63* 

E = XE 

000117 

00172 

64* 

XN  = XNLIM(  1 ) 

000120 

00173 

65* 

IC=0 

000122 

00174 

66* 

GO  TO  50 

000123  ‘ 

00175 

67* 

20  IF  (RHO  .GE.  RLI M( 2 ) ) GO  TO  10 

000125 

00177 

68* 

RHOsRHO+RLlM( 3) 

000130 

00200 

69* 

c=xc 

000133 

00201 

70* 

E = XE 

000135 

00202 

71* 

XN  = XNLIM(  1 ) 

000137 

00203 

72* 

IC=0 

000141 

00204 

73* 

GO  TO  50 

000142 

00205 

74* 

40  IF  (XN  .LT.  XNLIM(2 ) ) GO  TO  45 

000144 

00207 

75* 

IF  (SKG  .EQ.  0.0  .OR.  IPLOT  . EQ.  1)  GO  TO  20 

000147 

00211 

76* 

GO  TO  500 

000163 

00212 

77* 

45  XN=XN*XNL1M(3I 

000165 

00213 

78* 

C=AC*XN+BC 

000167 

00214 

79* 

E=AE*XN+BE 

000172 

00215 

80* 

50  RHOSQ=RHO*RHO 

000177 

00216 

81* 

RHOSQI* 1 .O/RHOSQ 

000201 

00217 

82* 

Z ( 1 ) =RHOSQ*RHOSQ* ( 2 . 0/ ( SQRT ( 1.0  RHOSQI ) ) -2 . 0-RHOSQI ) 

000204 

00220 

83* 

Z( 2) * 1 . 0+RHOSQI*Z( 1 ) 

000222 

00221 

84* 

Z ( 3) * ( 1 . 0+ RHOSQI * ( 1 .0-Z(1 ) )+RHOSQI*RHOSQI*Z(  1 ) ) 

000225 

00222 

85* 

H( 1 )=2.0*E/RHO+E*Z(2)/(RHOSQ*RHO) 

000236 

00223 

86* 

H( 2 ) *C*Z ( 3)/RH0+XN*Z ( 2 ) /RHO 

000251 

00224 

87* 

H(3)=E*Z( 2 )/RHOSQ 

000261 

00225 

88* 

H( 4 ) *C*Z ( 3 )+XN»Z( 1 )/RHOSQ 

000264 

00226 

89* 

H( 5 ) =H( 1 )+H( 2 ) 

000271 

00227 

90* 

H( 6 ) =RH0*H( 3 ) 

000274 

00230 

91* 

H( 7 ) =RHO*H ( 4 ) 

000277 

00231 

92* 

H(8)=H(3)-fH(4) 

000302 

00232 

93* 

YLAMSQ=-(H(1 )*H(4)*(H( 4 )— RH0*H( 2 ) ) )/(H( 2 ) *H(3) * ( RH0*H( 1 )-H(3) ) ) 

000305 

00233 

94* 

YLAM( 1 ) =SQRT ( Y LAMSQ ) 

000323 

00234 

95* 

Y '.AM  ( 2 ) =0 . 0 

000327 

00235 

96* 

Y LAM ( 3 ) * 1 . 0 

000331 

00235 

97* 

C 

000331 

00235 

98* 

C CALCULATION  OF  BETA 

000331 

00235 

99* 

C 

000331 

00235 

100* 

C 

000331 

00235 

101* 

C 

000331 

0023G 

102* 

A< 1 ,1 ) = Y LAMSQ*H( 1 )“H( 2 ) 

000333 

00237 

103* 

A(2,1)=YLAMSQ*H(3)-H(4) 

000337 

00240 

104* 

A( 1 , 2 ) =A(  2 ♦ 1 ) 

000343 

00241 

105* 

A ( 2 , 2 ) * Y LAMSQ*H( 6)“H( 7 ) 

000344 

00242 

106* 

183  B(1,1)=-YLAM(1)*H(5) 

000350 

00243 

107* 

8(2,1 ) *— Y LAM ( 1 ) «H(B ) 

000354 

00244 

108* 

V ( 1 ) = 1 .0 

000360 

00245 

109* 

CALL  DGJR( A. 12,12,2,2, $960 , UC , V) 

000361 

00246 

110* 

IF  (JC(1 ) .NE.  2)  GO  TO  950 

000373 

00250 

111* 

DO  190  1=1 ,2 

00041 1 

00253 

112* 

BETA( I)=0.0 

00041 1 

00254 

113* 

DO  185  K=1  ,2 

000415 

00257 

114* 

BETA(I)=BETA(I)+A(I,K)*B(K,1 ) 

000415 

00260 

1 15* 

185  CONTINUE 

0.0421 

00262 

116* 

BETA(I)=-BETA(I) 

000421 

00263 

117* 

190  CONTINUE 

000430 

00265 

1 18* 

A ( 1 , 1 J ="H( 1 ) -H( 2 ) 

000430 

00266 

119* 

A(2,1)=-H(3)-H(4) 

000432 

00267 

120* 

A(3, 1 )=0.0 

000434 

00270 

121* 

A(4, 1 )=0.0 

000436 

00271 

122* 

A(1.2)=A(2.1) 

000437 

00272 

123* 

A(2,2)=-H(6)-H(7) 

000441 

00273 

124* 

A(3, 21=0-0 

000445 

00274 

125* 

A(4,2)=0.0 

000446 

00275 

126* 

A ( 1 ,3)=-2.0*H(1 ) 

000447 

00276 

127* 

A ( 2 , 3 ) =“2 . 0*H( 3 ) 

000453 

00277 

128* 

A ( 3 , 3 ) ="H( 1 ) -H( 2 ) 

000457 

00300 

129* 

A(4.3)=-H(3)-H(4) 

000460 

00301 

130* 

A ( 1 , 4 ) = A ( 2 , 3 ) 

000461 

00302 

131* 

A(2,4)=-2.0*H(6) 

000463 

00303 

132* 

A(3,4)=A(2.1 ) 

000467 

00304 

133* 

A(4,4)=A(2.2) 

000471 

00305 

134* 

8(1,1 )=H(5) 

000473 

00306 

135* 

B( 2, 1 )=H(8) 

000475 

00307 

136* 

B( 3 . 1 )=0.0 

000477 

00310 

137* 

B(4,1)=0.0 

000500 

0031 1 

138* 

B(1.2)=-H(5) 

000S01 

0031" 

139* 

8(2,2) =-H( 8 ) 

000503 

00313 

140* 

B(3,2)*-H(5) 

000505 

00314 

141* 

8(4,2 ) =“H ( 8 ) 

000507 

00315 

142* 

V( 1 ) = 1 .0 

00051 1 

00316 

143* 

CALL  DGJR( A, 12,12,4,4, $960 , jC , V) 

000513 

00317 

144* 

IF  (JC(1 ) .NE.  4)  GO  TO  950 

000525 

00321 

145* 

DO  200  1=1 ,4 

000537 

00324 

146* 

DO  200  J*1 ,2 

000537 

00327 

1 47* 

T(  I . J)*0.0 

000537 

00330 

148* 

DO  200  K=1 .4 

000542 

00333 

149* 

T(I.J)*T(I,J)*A(I,K)*B(K,J) 

000542 

00334 

1 50* 

200 

CONTINUE 

000561 

00340 

151* 

8ETA(3J*-T(1 . 1 ) 

0n0561 

00341 

152* 

BETA(4)*-T(2. 1 ) 

000563 

00342 

153* 

SETA(5)=-T(1 ,2) 

000565 

00343 

154* 

BETA(6)*-T(2.2) 

000567 

03343 

155* 

C 

000567 

00343 

156* 

c 

UNIFORM  LOAD 

000567 

00343 

157* 

c 

000567 

00344 

158* 

IJUNC*0 

000571 

00345 

159* 

DO  210  1*1 ,6 

000576 

00350 

160* 

BLAM2(I)=BETA(I)*YLAM(2) 

000S76 

00351 

161* 

BLAM3 (I)*BETA(I ) • YLAM( 3 ) 

000600 

00352 

162* 

DO  210  U*1 ,6 

000605 

00355 

163* 

A(I.J)*0.0 

000605 

00356 

164* 

210 

CONTINUE 

000614 

00361 

165* 

IF  (F  .EO.  0.0)  GO  TO  300 

000614 

00363 

166* 

CALL  CALCA  (ALPHA) 

000616 

00364 

167* 

DO  212  1*1 .3 

000626 

00367 

168* 

DO  21 1 J*1 ,6 

000626 

00372 

169* 

TARR( 1+3, J)*ARR( I , J) 

000626 

00373 

170* 

211 

CONTINUE 

000637 

00375 

171* 

212 

CONTINUE 

000637 

00377 

172* 

CALL  CALCA  (-ALPHA) 

000637 

00400 

173* 

DO  217  1*1 .3 

000651 

00403 

174* 

DO  216  J=1 ,6 

000651 

00406 

175* 

TARR(I,J)»ARR(I,J) 

000651 

0040 

176* 

216 

CONTINUE 

000662 

0041 1 

177* 

217 

CONTINUE 

000662 

00413 

178* 

PRINT  1011 

000662 

00415 

179* 

I F ( IEOUN.EQ. 1 ) PRINT  1041 

000666 

00420 

180* 

IF< IBOUN.E0.0)  PRINT  1042 

000675 

00423 

181* 

1042 

FORMAT ( 59X ‘ C LAMPED  ENDS') 

000703 

00424 

182* 

1041 

FORMAT ( 56X ' SIMPLY  SUPPORTED  ENDS') 

000703 

00425 

183* 

PRINT  1010  01  DEN 

000703 

00430 

184* 

PRINT  1040  AC.BC.AE.BE 

000712 

00436 

185* 

1040 

FORMAT ( /5X, *C*' , F6-4. '*N*' . F6.4/5X , '£■' , F6. 3 , ' *N* ' .F6.4) 

000722 

00437 

186* 

PRINT  1013,  ALPHA . RHO , E , C , XN 

000722 

00446 

187* 

PRINT  1014 

000733 

00446 

188* 

C 

000733 

00446 

189* 

c 

SOLVE  FOR  G 

000733 

00446 

190* 

c 

000733 

00450 

191* 

G( 1 ) *-F/W( 1 ) 

000737 

00451 

192* 

1 ARR ( 1 . 7 ) =0 . 0 

000743 

00452 

193* 

T ARR ( 2 . 7 ) >0 . 0 

000745 

00453 

194* 

T ARR( 3 , 7 ) *-G( 1 ) 

000746 

00454 

195* 

TARR(4.7)«0.0 

000747 

00455 

196* 

T ARR ( 5 , “ ) *0 . 0 

000750 

00456 

197* 

T ARR (6,7) *-G( 1 ) 

000751 

00456 

198* 

C UNIFORM  LOAD  IS  REDUCED  TO  A PROBLEM  OF  ORDER  3 BV  STMMETRY 

000751 

00457 

199* 

00  221  1*1,3 

000755 

00462 

200* 

TARR( I . 2)*TARR( 1,4) 

000755 

00463 

201* 

T ARR( 1 , 3 ) *TARR( 1.5) 

000756 

00464 

202* 

221 

TARR( I,4)*TARfl(I.7) 

000760 

r 


00466 

203* 

IF(IBOUN.E0.1 ) T ARM (2,4)*  G(1) 

000763 

00470 

204* 

V ( 1 ) *4 . 0 

000770 

00471 

205* 

CALL  DGJR(TARR, 1 3 . 1 2 , 3 . 4 , $960 . JC.V) 

000772 

00472 

206* 

0(1) *TARR (1,4) 

001004 

00473 

207* 

D( 2) *0. 0 

001006 

00474 

208* 

D( 3) *0 • 0 

001010 

00475 

209* 

0(4 ) =TARR( 2 . 4 ) 

00101 1 

00476 

210* 

0(5) =TARR( 3, 4 ) 

001013 

00477 

211* 

D(  6) *0 ■ 0 

001015 

00500 

212* 

ITYP=1 

001016 

0050 

213* 

CALL  CALPAR(-ALPHA.NMAX.THETA,*FA,0,0) 

001020 

00502 

214* 

IF  (SMG  .EO.  0.0)  GO  TO  40 

001032 

00502 

215* 

C 

001032 

00502 

216* 

C GREEN'S  FUNCTION 

001032 

00502 

217* 

C CONCENTRATED  LOAD 

001032 

00502 

218* 

c 

001032 

00504 

219* 

300  CONTINUE 

001040 

00505 

220* 

3005  FORMAT ( 14 ) 

001040 

00506 

221* 

DO  303  ICL*1 i 200 

001040 

0051 1 

222* 

303  SOL( ICL ) *0 . 0 

001040 

00513 

223* 

LIN*0 

001042 

00514 

224* 

I F ( INTEG.EQ.O)  GO  TO  250 

001043 

00516 

225* 

READ( 1 0 . 3005 ) NIS 

001045 

00521 

226* 

252  LIN.LIN+1 

001054 

00522 

227* 

I JUNC*0 

001056 

0C523 

228* 

RE AD( 1 0 . 3006 ) THETA. MFA 

001057 

00527 

229* 

1 HET  A*THETA* ALPHA 

001066 

00530 

230* 

HFA*0.0174S*ALPHA*HFA 

001071 

00531 

231* 

3006  FORMAT  (2E18.8) 

001076 

00532 

232* 

250  CONTINUE 

001076 

00533 

233* 

1 F ( INTEG. EO. 1 ) GOTO  251 

001076 

00535 

234* 

PRINT  1002 

001100 

00537 

235* 

1 F ( ILOAD.EQ.O) PRINT  1021 

001104 

00542 

236* 

1 F ( I LOAD . EO* 1 ) PRINT  1022 

001112 

00545 

237* 

1021  FORMAT ( 61 X . 8H( NORMAL ) ) 

001 I2t 

00546 

238* 

1022  FORMAT ( 60X , 1 OH ( VERT ICAL ) ) 

001121 

00547 

239* 

IF( IBOUN.E0.1 ) PRINT  1041 

001121 

00552 

240* 

IF(IBOUN.EQ.O)  PRINT  1042 

001130 

00555 

241* 

PRINT  1010  oiden 

001136 

00560 

242* 

PRINT  1040  AC.BC.AE.BE 

001145 

00566 

243* 

print  1013,  ALPHA. RHO.E.C.XN 

001155 

00575 

244* 

PRINT  1003,  THETA 

001166 

00600 

245* 

PRINT  1014 

001173 

00602 

246* 

251  CONTINUE 

001200 

00603 

247* 

G( 1 )*0.0 

001200 

00604 

248* 

DO  310  1*1 ,3 

001205 

00607 

249* 

DO  310  J*7.12 

001205 

00612 

250* 

A( I , J ) *0 ■ 0 

001205 

00613 

251* 

A(I*-3,J-6)=0.0 

001206 

00614 

252* 

310  CONTINUE 

001214 

00617 

253* 

CALL  CALCA(O.O) 

001214 

00620 

254* 

DO  320  1*1,3 

001224 

00623 

255* 

DO  315  J*1 ,6 

001224 

00626 

256* 

A(I,J)*ARR(I,J) 

001224 

00627 

257* 

315  CONTINUE 

00123S 

00631 

258* 

320  CONTINUE 

001235 

00633 

259* 

CALL  CALCA(O.O) 

001235 

130 


00634 

260* 

00  324  1*4,6 

001245 

00637 

261* 

DO  323  d=7 , 12 

001245 

00642 

262* 

A(I.d)=ARR(I-3.d-6) 

001245 

00643 

263* 

323 

CONTINUE 

001256 

00645 

264* 

324 

CONTINUE 

001256 

00647 

265* 

IdUNC=1 

001256 

00650 

266* 

T"MPA=THETA 

001260 

00651 

267* 

KSET  *1 

001262 

00652 

268* 

I F ( KSET . EQ. 1 ) GO  TO  329 

001263 

00654 

269* 

328 

KSET =KSET+1 

001266 

00655 

270* 

DO  360  J=1  .6 

001273 

00660 

271* 

A(10,d+6)*-A(10,d) 

001273 

00661 

272* 

A ( 1 1 , d+6 ) =-A (11 ,d) 

001274 

00662 

273* 

A( 12,d+6)*A( 12, d) 

001276 

00663 

274* 

360 

CONTINUE 

001301 

00665 

275* 

THETA=ALPHA+TEMPA 

001301 

00666 

276* 

CALL  CALCA(THETA) 

001304 

00667 

277* 

DO  330  1*1 ,3 

001313 

00672 

278* 

K=  I +6 

001313 

00673 

279* 

DO  325  d=1  ,6 

001316 

00676 

280* 

A(K,d)*ARR(I,d) 

001324 

00677 

281* 

325 

CONTINUE 

001333 

00701 

282* 

330 

CONTINUE 

001333 

00703 

283* 

GO  TO  306 

001333 

00704 

284* 

329 

CONTINUE 

001335 

00705 

285* 

T HE TA= ALPHA-THETA 

001335 

00706 

286* 

CALL  CALCA(THETA) 

001337 

00707 

287* 

DO  305  d=1  .6 

001353 

00712 

288* 

■A  ( I , d+6  ) *— A(  I , d ) 

001353 

00713 

289* 

A ( 7 , d+6 ) = ARR ( 1 ,d) 

001354 

00714 

290* 

A ( 8 , d+6 ) = ARR ( 2 , d) 

001356 

00715 

291* 

A ( 9 , d+6) *- ARR ( 3 , d ) 

001360 

00716 

292* 

305 

CONTINUE 

001364 

00720 

293* 

304 

CONTINUE 

001364 

00721 

294* 

306 

CONTINUE 

001364 

00722 

295* 

PA=H(1 )*(BETA(1) *YLAM( 1 ) — 1 • 0 ) +H (3)*BETA(2)*Y  LAM ( 1 ) 

001364 

00723 

296* 

A ( 1 0 , 1 ) *PA*C5A1 

001375 

00724 

297* 

A( 10,2)=PA*SNA1 

001377 

00725 

298* 

PA  = BETA(1 )*YLAM( 1 )-1 . 0+RHO*BETA ( 2 ) *YLAM(  1 ) 

001402 

00726 

299* 

A ( 1 1 ,1 )*PA*CSA1 

001407 

00727 

300* 

A ( 1 1 ,2)*PA*SNA1 

001411 

00730 

301* 

P A*H( 2 ) * ( YLAM( 1)+8ETA(1)) +HI 4) *BET  A(2) 

001414 

00731 

302* 

A ( 12.1 ) * PA*SNA1 

001423 

00732 

303* 

A( 12.2)*PA*CSA1 

001425 

00733 

304* 

PA=(BETA(3)-1 .0)*H( 1 )+H( 3 ) * BETA ( 4 ) 

001430 

00734 

305* 

PB*(BETA(5)-BETA(3))*H(1 )+H( 3) * ( BETA( 6)-BETA (4) ) 

001440 

00735 

306* 

THETR= THETA 

001452 

00736 

307* 

THET  A = 0 . 01 745* THETA 

001454 

00737 

308* 

A( 10.3)=PA*SNA3 

001456 

00740 

309* 

A( 1 0 , 4 ) *PA*CSA3 

001460 

0074' 

310* 

A(  10.5) *PB«CSA3+THETA*A(  10.3) 

001463 

00742 

311* 

A< 10,6) *-PB*SNA3+THETA*A( 10,4) 

001467 

00743 

312* 

PA*BETA(3)  — 1 . 0+RHO*  8ETA  ( 4 ) 

001475 

00744 

313* 

PB*(BETA(5)-BETA(3) )+RHO* ( BETA( 6 )-BET A( 4 ) ) 

001501 

00745 

314* 

A ( 1 1 ,3) *PA*SNA3 

001S05 

00746 

315* 

A ( 1 1 ,4)*PA*CSA3 

001507 

00747 

316* 

A ( 1 1 ,5)*PB*CSA3+THETA*A(1 1 ,3) 

001512 

00750 

317* 

A( 1 1 ,6)*-PB*SNA3+THETA*A( 11,4) 

001517 

00751 

318* 

PA=H(2)*(1 .0-BETA(3) )-H( 4 ) * BETA ( 4 ) 

001525 

00752 

319* 

PB=H(2)*(1.0+BETA(5))+H(4)*BETA(6) 

001535 

00753 

320* 

A(  1 2 , 3 ) =PA*CSA3 

001544 

00754 

321* 

A( 1 2 , 4 ) -—PA* SNA 3 

001547 

0075'’ 

322* 

A(  1 2 , 5>  = Aj  12,3)  *THET  A-fPB*  SNA3 

001553 

0075b 

323* 

A(12.6)=A(12,4)*ThETA+PB*CSA3 

001560 

00757 

324. 

IFIKSET.EQ.1 ) GO  TO  328 

001S66 

00761 

325* 

theta=tempa 

001571 

0-7762 

326* 

DO  370  1*1,12 

001577 

00765 

327* 

DO  365  d=1 ,12 

001577 

00770 

328* 

DA(I.J)*A(I,d) 

001577 

00771 

329* 

365 

CONTINUE 

001610 

00773 

330* 

370 

CONTINUE 

001619 

00775 

331* 

DO  375  1=1,11 

001610 

01000 

332* 

DA( I , 1 3)=0 . 0 

001610 

01001 

333* 

375 

CONTINUE 

001612 

01003 

334* 

THEHAR=0.01745«THEHAT 

001612 

01004 

335* 

THETAR=0 . 01 745*THET  A 

001615 

01005 

336* 

QSBsS 1N( THET  AR ) 

001620 

01006 

337* 

0C8=C0S ( THET  AR ) 

001624 

01007 

338* 

QSH=SIN(THEHAR) 

001630 

01010 

339* 

QCH=C0S(THEHAR) 

001634 

0101 1 

340* 

IF(ILOAD.EQ.I)  GO  to  376 

001640 

01013 

341* 

0G 1:0.0 

001643 

01014 

342* 

QG2*0.0 

001645 

01015 

343* 

QG3-SMG 

001646 

01016 

344* 

0G4=0 . 0 

001650 

01017 

345* 

GO  TO  377 

001651 

01020 

346* 

376 

0G1 =-SMG*OSM*OSH*OCB 

0016S3 

0102 

347* 

0G2=-SMG*0S8 

001660 

01022 

348* 

QG3*SMG*QCH*QCH*QCB 

001664 

01023 

349* 

0G4*SMG*0CH*0SB 

001672 

01024 

350* 

377 

DA( 10,1 3) *002 

001676 

01025 

351* 

DA ( 1 1 , 1 3 ) =_QG4/H( 3) 

001677 

01026 

352* 

DA ( 1 2 , 1 3 ) =0G1 “0G3 

001703 

01027 

353* 

V(  1 )=4.0 

001706 

01030 

354* 

CALL  DGJR( DA ,13,12,12,1 3,5960 ,JC,V) 

001710 

01031 

355* 

G( 1 ) =0 .0 

001722 

01032 

356* 

IF  ( IP LOT  -1)  394,395,394 

001724 

01035 

357* 

394 

ISAV=1 

001727 

01036 

358* 

GO  TO  396 

001731 

01037 

359* 

395 

ISAV=0 

C01733 

01040 

360* 

396 

NVAL=0 

001734 

01041 

361* 

I TYP=0 

001734 

01042 

362* 

CALL  CALPAR(-ALPHA,NMAX,THETA.WFA,ISAV,NVAL) 

001735 

01043 

363* 

Ir((LIN.LT.NIS). AND . ( I NT EG . EO- 1 ) ) GO  TO  252 

001747 

01045 

364* 

IF(INTEG.EO.O)  GO  TO  401 

001764 

01047 

365* 

REWIND  10 

001766 

01050 

366* 

PRINT  1012 

001771 

01052 

367* 

IF( 180UN.E0. 1 ) PRINT  1041 

001775 

01055 

368* 

IF(IBOUN.EO.O)  PRINT  1042 

002004 

01060 

369* 

PRINT  1010.QIDEN 

002012 

01063 

370* 

PRINT  1040  AC.BC.AE.BE 

002021 

01071 

371* 

PRINT  1013, ALPHA. RHO.E.C.XN 

002031 

01100 

372* 

PRINT  1014 

002042 

01102 

373* 

JP*8«(NMAX+1  ) 

002046 

ono3 

374* 

PRINT  1009. ( (SOL( IP) ) . IP*1 . dP) 

002052 

01111 

375* 

1012  FORMAT ( ' 1 • ,45X , ' DEFORMATION  OFPRESSURE  STABILIZED  ARCHES') 

002066 

01112 

376* 

1010  FORMAT ( 1 X , 1 3A6 ) 

002066 

01113 

377* 

401  CONTINUE 

002066 

01114 

378* 

1009  FORMAT  1 /OPF 10.3.1P7E17.8) 

002066 

01115 

379* 

GO  TO  40 

002066 

01116 

380* 

500  END  FILE  IS 

002067 

01117 

381* 

REWIND  IS 

002071 

01120 

382* 

RRHO=RHO 

002074 

01121 

383* 

00  503  1*1.4 

002104 

0112  • 

384* 

WSCAL< I)*0.0 

002104 

01125 

385* 

XNSCA  L ( I ) *0 . 0 

002104 

01126 

386* 

WLSCAL ( I ) =0 . 0 

002105 

01127 

387* 

503  CONTINUE 

002112 

01131 

388* 

DO  508  d*1 . IC 

002112 

01134 

389* 

READ  (15.3003)  CN( d ) . CW( d ) . CWL( d ) 

002112 

01141 

390* 

IF  ( CN( J ) .LT.  XNSCAL(I))  XNSCA L ( 1 ) *CN( J ) 

002121 

01143 

391* 

IF  ( CN ( d ) .GT.  XNSCAL ( 2 ) ) XNSC A L ( 2 ) =CN( d ) 

002127 

01145 

3S2* 

IF  (CW(d)  .LT.  WSCAL(I))  WSCAL ( 1 ) =CW ( d ) 

002135 

01147 

393* 

IF  (CW(d)  .GT.  WSCAL ( 2 ) ) WSCAL ( 2) *CW( d ) 

002143 

01151 

394* 

IF  (CWL(d)  .LT.  WLSCAL(I))  WLSCAL( 1 ) *CWL( d ) 

002151 

01153 

395* 

IF  (CWL(d)  .GT.  WLSCAL(2))  WLSCAL ( 2 ) *CWL( d ) 

002157 

01155 

396* 

508  CONTINUE 

002167 

01157 

397* 

REWIND  15 

002167 

01160 

398* 

NEXP=0 

002172 

01161 

399* 

IE0F*0 

002173 

01162 

400* 

REWIND  16 

002174 

01163 

401* 

510  READ  ( 16.3004.END*530)  AR . AN . AWI , AWO, AWRL 

002200 

01172 

402* 

IF  ( ABS(AR-RRHO)  .GT.  0.005)  GO  TO  540 

002212 

01174 

403* 

NEXP*NEXP*1 

002220 

01175 

404* 

IF  (AN  .LT.  XNSCAL(I))  XNSCAL(1)*AN 

002223 

01177 

405* 

IF  (AN  .GT.  XNSCALI 2 ) ) XNSCAL(2)*AN 

002231 

01201 

406* 

IF  (AWI  .LT.  WSCAL(I)  .AND.  AWI  .NE.  0.0) 

WSCAL( 1 )*AWI 

002237 

01203 

407* 

IF  (AWO  .LT.  WSCAL! 1)  .AND.  AWO  .NE.  0.0) 

WSCALt 1 )*AWO 

002254 

01205 

408* 

IF  (AWI  .GT.  WSCALt 2 ) .AND.  Awl  .NE.  0.0) 

WSCAL ( 2 ) « AWI 

002271 

01207 

409* 

IF  (AWO  .GT.  WSCAL( 2 ) .AND.  AWO  .NE.  0.0) 

WSCAL ( 2 ) • AWO 

002306 

0121 1 

410* 

IF  (AWRL  .NE.  0.0)  GO  TO  520 

002323 

01213 

411* 

IF  (NEXP  .NE.  1)  GO  TO  516 

002325 

01215 

412* 

ELN( 1 ) *0 . 0 

002330 

01216 

413* 

EWL( 1 ) *0 . 0 

002331 

01217 

414* 

GO  TO  525 

002332 

01220 

415* 

518  EWL( NEXP ) *EWL ( NEXP- 1 ) 

002334 

01221 

416* 

ELN(NEXP)«ElN(NEXP-1  ) 

002336 

01222 

417* 

GO  TO  525 

002340 

01223 

418* 

520  IF  (AWRL  .LT.  WLSCAL(I))  WLSCAL ( 1 ) *AWRL 

002342 

01225 

419* 

IF  (AWRL  .GT.  WLSCAL(2) ) WLSCAL ( 2 ) *AWRL 

002351 

01227 

420* 

EWL (NEXP ) * AWRL 

002357 

01230 

421* 

E LN ( NEXP ) * AN 

002362 

01231 

422* 

525  EN(NEXP)*AN 

002365 

O' 232 

423* 

EWO( NEXP)* AWO 

002367 

01233 

424* 

EWI (NEXP)*AWl 

002371 

01234 

425* 

GO  TO  510 

002373 

01235 

426* 

530  IF  (NEXP  .NE.  0 .OR.  IEOF  .NE.  1)  GO  TO  550 

002375 

01237 

427* 

I EOF* 1 

002406 

01240 

428* 

REWINO  16 

002410 

01241 

429* 

GO  TO  510 

002413 

01242 

430* 

540  IF  (NEXP  .EO.  0)  GO  TO  510 

002415 

01244  431* 

01245  432* 

01 240  433* 

01247  434* 

01250  435» 

01251  436* 

01252  437* 

01253  438* 

01254  439* 

01255  440* 

01257  441* 

01260  442* 

01261  443* 

01261  444* 

01262  445* 

01262  446* 

01263  447* 

01264  448* 

01265  449* 

01266  450* 

01267  451* 

01270  452* 

01271  453* 

01272  454* 

01274  455* 

01275  456* 

01277  457* 

01300  458* 

01301  459* 

01302  460* 

01303  461* 

01304  462* 

01305  463* 

0130G  464* 

01310  465* 

01311  466* 

01312  467* 

01313  468* 

01314  469* 

01315  470* 

01316  471* 

01317  472* 

01320  473* 

01321  474* 

01321  475* 

01322  476* 

01322  477* 

01323  478* 

01324  479* 

01325  480* 

0132:  481* 

01327  482* 

01330  483* 

01331  484* 

01332  485* 

01334  486* 

01335  487* 


550  CALL  SCALE  ( ANSCAL . 8 . 0 . 2 , 1 ) 

CALL  SCALE  (a^SCAL. 8.0,2, 1 ) 

CALL  SCALE  ( WSCAL , 8 ■ 0 , 2 , 1 ) 

CN( IC+1 )*XNSCAL(3) 

CM  IC+2 ) * XNSCAL( 4 ) 

CM  IC+1  ) *WSC A L(  3 ) 

CW( IC+2 ) *WSC AL ( 4 ) 

CWL( IC+1 )*WLSCAl(3> 

CWL( 1C+2)=WLSCAL(4) 

IF  (NEXP  . EQ . 0)  GO  TO  570 
EMNEXP+1  ) *CN(  IC+1  > 

EM NEXP+2 ) *CM  IC+2) 

570  CALL  AXIS  (0.0.0.0, 'N0ND1MENSI0NAL  PRESSURE ', -23 . 5 . 0 . 0 . 0 , CM IC+1 ) , 
1 CM  IC+2)  ) 

CALL  AXIS  (0.0. 0.0, 'NONDIMENSIONAL  FLEXI B1 LITY • . 26 ,6. 0 , 90 . 0 , 

1 CW(  IC+1  ) .CM  IC+2)  ) 

CALL  SYMBOL  ( 3 . 0 . 9 - 5 . 0 . 21 , ' RHO  **,0.0,5) 

CALL  NUMBER  ( 4 . 5 , 9- 5 . 0 . 21 , RRHO, 0 . 0 , 4 ) 

CALL  NEWPEN  (2) 

CALL  PLOT  (1 .0,9. 0.3) 

CALL  PLOT  (1 .5. 9. 0.2) 

CALL  SYMBOL  ( 1 . 6 . 9. 0 , 0 . 1 4 . * CALC  VAL',0.0,8) 

CALL  FLINE  ( CN , CW , I C , 1 , 0 , 0 ) 

IF  (NEXP  .EQ.  0)  GO  TO  600 
CALL  SCALE  (EM ,8.0. NEXP, 1 ) 

IF  ( EWI INEXP+ I ) .EQ.  0.0  .ANO.  EWI(NEXP+2)  .EQ.  0.0)  GO  TO  580 
EWI(NEXP+1 ) =CW( IC+1  ) 

EW I ( NEXP+2 ) *CW( IC+2 ) 

CALL  NENPEN  (3) 

CALL  LINE  ( EN , EM , NEXP , 1 , - 1 , 0 ) 

CALL  SYMBOL  ( 3 . 605 . 9 . 07 . 0 . 1 4 . 0 . r . 0 , -1 ) 

CALL  SYMBOL  ( 3 . 8 , 9 . 0 , 0 . 1 4 , * EXP  IN', 0.0, 6) 

580  CALL  SCALE  ( EWO , 8 . 0 . NEXP , 1 ) 

IF  ( EWO( NEXP+1 ) .EQ.  0.0  .ANO.  EW0(NEXP+2)  .EQ.  0.0)  GO  TO  600 
EW0(NEXP+1 ) = CM IC+1 ) 

EWO( NEXP+2 ) = CW( IC+2) 

CALL  NEWPEN  ( 1 ) 

CALL  SYMBOL  (5.4.9.07,0.14,1,0.0,-1) 

CALL  SYMBOL  ( 5 . 595 . 9 . 0 , 0 . 1 4 , * EXP  OUT',0.0,7) 

CALL  LINE  ( EN , EWO ,NEXP,1 ,-1 ,2) 

GO  TO  610 

600  CALL  NEWPEN  (1  ) 

610  CALL  PLOT  ( 10. 0,0. 0.-3) 

CALL  AXIS  (0.0, 0.0. ‘NONDIMENSIONAL  PRESSURE ’ ,-23 , 5 . 0 , 0 . 0 ,CM IC+1 ) , 
1 CM  IC+2)) 

CALL  AXIS  (0.0, 0.0. ‘NONDIMENSIONAL  WRINKLING  LOAD ' , 29 , 6 . 0 , 90 . 0 , 

1 CWL( IC+1 ) ,CWL( IC+2) ) 

CALL  SYMBOL  ( 3 . 0 , 9 . 5 , 0 . 21 , ' RHO  *',0.0, 5) 

CALL  NUMBER  ( 4 . 5 . 9. 5 . 0 . 21 . RRHO, 0 . 0 , 4 ) 

CALL  NEWPEN  (2) 

CALL  PLOT  (1 .0.9. 0.3) 

CALL  PLOT  (1 .5. 9. 0.2) 

CALL  SYMBOL  ( 1 . 6 . 9 • 0 . 0 . 1 4 . ‘ CALC  VAL‘,0.0.8) 

CALL  FLINE  ( CN.CWL. IC , 1 , 0 , 0 ) 

IF  (NEXP  .EQ.  0)  GO  TO  650 
CALL  SCALE  ( EWL , 8 . 0 , NEXP , 1 ) 

IF  ( EWL ( NEXP+1 ) .EQ.  0.0  .ANO.  EWL(NEXP+2)  .EQ.  0.0)  GO  TO  650 


002417 

002424 

002432 

002440 

002443 

002445 

002447 

002451 

002453 

002455 

002457 

002462 

002464 

002464 

002502 

002502 

002521 

002531 

002541 

002544 

002551 

002556 

002566 

002576 

002600 

002606 

0-2620 

002622 

002624 

002627 

002637 

002647 

002660 

002665 

002677 

002701 

002703 

002706 

002716 

002726 

002736 

002740 

002743 

002747 

002747 

002766 

002766 

003005 

003015 

003025 

003030 

003035 

003042 

003052 

003062 

003064 

003072 


0 1 337 

480* 

EWL(NEXP+1  ) sCWL ( IC+1  ) 

003104 

01340 

489* 

EWL(NEXP+2)=CWL( IC+2) 

003106 

01344 

490  • 

ELN( NEXP+1 ) = CN( IC+1 ) 

003110 

01342 

491  • 

ELN( NCXP+2 ) *CN( IC+2) 

003112 

01343 

492* 

CALL  MEWPEN  (3) 

003114 

01344 

493* 

CALL  LINE  (ELN.EWL.NEXP. 1 ,-1 ,5) 

003117 

0 1 345 

494* 

CALL  SYMBOL  (3.605.9.07,0.14,14.0.0,-1) 

003127 

01346 

495* 

Call  symbol  (3. 8, 9.0,0. 14. ' EXP  VAL', 0.0,7) 

003137 

01347 

496* 

650  CALL  NEWPEN  ( 1 ) 

003150 

01350 

497* 

CALL  PLOT  (10. 0,0. 0,-3) 

003152 

01351 

498* 

GO  TO  20 

003157 

01352 

499* 

960  PRINT  1016 

003161 

01354 

500* 

GO  TO  40 

003164 

01355 

501* 

950  PRINT  1017 

003166 

01357 

502* 

GO  TO  40 

003171 

01360 

01360 

•DIAGNOSTIC* 

503* 

CONTROL  CAN  NEVER  REACH  THE  NEXT  STATEMENT 

970  PRINT  1004 

003172 

Cl  362 

504* 

GO  TO  40 

003176 

01363 

01363 

•DIAGNOSTIC* 

505* 

CONTROL  CAN  NEVER  REACH  THE  NEXT  STATEMENT 

980  PRINT  1005 

003177 

01365 

506* 

GO  TO  40 

003203 

01366 

507* 

990  PRINT  1006 

003205 

01370 

508* 

IF  ( I PLOT  .ED.  1)  CALL  EXIT 

003210 

01372 

509* 

CALL  PLOT  (4.0.0.0,999) 

003215 

01373 

510* 

CALL  EXIT 

003222 

01374 

511* 

1001  FORMAT  (3F10.0. 10X.3F10.0) 

003226 

01375 

512* 

1002  FORMAT  (• 1' , 45X ,* DE FORMAT  ION  OF  PRESSURE  STABILIZED  ARCHES'// 

003226 

01375 

513* 

1 57X, 'CONCENTRATED  LOAD') 

003226 

01376 

514* 

1003  FORMAT  (13X,'L0A9  LOCATION  =',F10.3) 

003226 

01377 

515* 

1004  FORMAT  ( ' 1 ROOTS  DO  NOT  CONVERGE  IN  100  ITERATIONS') 

003226 

0140  • 

516* 

1005  FORMAT  ('1  INTEGRAL  DOES  NOT  CONVERGE  IN  100  ITERATIONS') 

003226 

01401 

517* 

1006  FORMAT  ('1  END  OF  JOB') 

003226 

01402 

518* 

1011  FORMAT  (' 1 ' ,45X, 'DEFORMATION  OF  PRESSURE  STABILIZED  ARCHES'// 

003226 

01402 

519* 

1 49X, 'NORMAL  UNIFORM  LOAD, DIRECT  SOLUTION') 

003226 

01403 

520* 

1013  FORMAT  ( //1 1 X , ’ HALF  SPAN  ANGLE  *', FI 0 . 3/1 4X ' RADIUS  RATIO  *',F10.3/ 

003226 

01403 

521* 

1 5X , ' ELAST 1C  MODULUS  RATIO  «',F10.3/ 

003226 

01403 

522* 

2 7X , ' SHEAR  MODULUS  RATIO  s',Fl0.3/ 

003226 

01403 

523* 

3 8X, 'PRESSURE  PARAMETER  «',F10.3) 

003226 

01404 

524* 

1014  FORMAT (///3X, ' PSI ' . 1 3X , 'U ' , 1 6X , ' W , 15X , ■ PHI ' , 1 5X , ' E( 0) ' 

1 1 2X , ' E(  PI  ) ' , 1 1 X,  ' N ' , 1 7X , ' M 1 ) 

1015  FORMAT  ('  ERROR  IN  LAMBDA-SQ  ' .6E18.8) 

003226 

01404 

525* 

003226 

01405 

526* 

003226 

01406 

527* 

1016  FORMAT  ('  OVERFLOW  IN  GJR  ') 

003226 

01407 

528* 

1017  FORMAT  ('  SINGULARITY  IN  GUR  ') 

003226 

01410 

529* 

1020  FORMAT  (7F10.3.I5) 

003226 

0141 1 

530» 

2001  FORMAT  (2F10.0.I5) 

003226 

01412 

531* 

3001  FORMAT (411) 

003226 

01413 

532* 

3002  FORMAT  (7F10.0) 

003226 

01414 

533* 

3003  FORMAT  (3E20.8) 

003226 

01415 

534* 

3004  FORMAT  (5E15.6) 

003226 

01416 

535* 

END 

003226 

END  OF  COMPILATION 


2 DIAGNOSTICS 


•FOR.S  TPFS. ARCHRESULTS, R 
FOR  OOE3-1 1/03/77-16:01 :36  (0.) 

SUBROUTINE  CALPAR  ENTRY  POINT  001046 

STORAGE  USED:  CODE(1)  00107S;  DATA ( 0)  000222;  BLANK  COMMON(2)  001541 


E TERNAL  REFERENCES  (BLOCK,  NAME) 

0003  DCOS 

0004  DS IN 

0005  DCOSH 

0006  DSINH 

0007  NPRT$ 

0010  NI02S 

0011  NWOUS 

0012  NERR3S 


STORAGE  ASSIGNMENT  (BLOCK,  TYPE,  RELATIVE  LOCATION,  NAME) 


0001 

000750 

1 0L 

0001 

001025 

1 OOL 

OOCO 

000101 

1 001  F 

0000 

000104 

1002F 

0001 

000764 

1 1 L 

0001 

000150 

1 16G 

0001 

000174 

125G 

0001 

000201 

3L 

0001 

000203 

4L 

0001 

000222 

5L 

0001 

000702 

6L 

0002 

D 

000007 

ARR 

0002 

D 

000053 

BETA 

0002 

D 

000175 

BLAM2 

0002 

D 

00021 1 

BLAM3 

0000 

D 

000075 

CINT 

0000 

D 

000024 

CS 

0002 

D 

00011 1 

CSA1 

0002 

0 

000115 

CS  A3 

0000 

D 

000032 

CSP 

0002 

D 

000161 

D 

0002 

0 

000225 

DA 

0000 

D 

000006 

DELTH 

0002 

D 

000005 

E 

0000 

D 

000052 

EPI 

0000 

D 

000050 

EZERO 

0000 

D 

000002 

FL 

0002 

D 

000075 

G 

0002 

0 

000125 

H 

0002 

D 

CO0 1 2 1 

HCA2 

0002 

D 

000123 

HSA2 

0002 

001537 

I BDUN 

0002 

I 

000004 

IC 

0000 

I 

000015 

ID 

0002 

001540 

I JUNC 

0000 

000123 

1NJP$ 

0002 

I 

001535 

INTEG 

0002 

I 

001536 

ITYP 

0000 

I 

000010 

N 

0000 

I 

000074 

N I 

0000 

D 

000044 

PA 

0000 

D 

000046 

PB 

0000 

D 

000060 

PC 

0000 

D 

000062 

PD 

0000 

D 

000064 

PE 

0000 

D 

000066 

PF 

0000 

D 

000070 

PG 

0000 

D 

000072 

PH 

0000 

D 

000040 

PHI 

0000 

D 

000016 

PSI 

0000 

D 

000013 

PSIDEG 

0000 

0 

000030 

PI 

0000 

D 

000022 

P2 

0000 

D 

000020 

P3 

0000 

D 

000054 

OM 

0000 

D 

000056 

ON 

0000 

D 

000004 

QNSTEP 

0002 

D 

000000 

RHO 

0002 

D 

000002 

SMN 

0000 

D 

000026 

SN 

0002 

D 

000113 

SNA1 

0002 

D 

000117 

SNA3 

0000 

D 

000034 

SNP 

0002 

D 

000715 

SOL 

0000 

D 

000036 

U 

0000 

3 

000042 

W 

0000 

D 

000077 

WRL 

0002 

D 

000145 

XI 

0000 

D 

00001 1 

XN 

0002 

D 

000153 

XR 

0000 

D 

000000 

XXTH 

0002 

D 

000067 

Y LAM 

00101 

1* 

SUBROUTINE  CALPAR  ( DEGREE , NSTEP . THETA , WF A . I SAV . IND ) 

000000 

00103 

2* 

IMPLICIT  DOUBLE  PRECISION(A-H.O-Z) 

000000 

00104 

3* 

COMMON  RHO , SMN , IC , E 

onoooo 

00105 

4* 

COMMON  ARR(3,6) ,BETA(S) , YL AM (3) , G ( 6 ) , CSA 1 ,SNA1 , CSA3 , SNA3 

000000 

00106 

5* 

COMMON  HC A2 , HSA2 ,H( 8 ) 

000000 

00107 

6* 

COMMON  XI(3),XR(3),D(6), BLAM2( 6 ) , BLAM3 ( 6 ) 

000000 

00110 

7* 

COMMON  DA( 12, 13) ,SOL(200) , INTEG, ITYP, IBOUN, I JUNC 

000000 

00111 

B» 

XXTH=. 01745 *THETA 

000000 

00112 

9* 

FL* . 0326»RHO/ ( RHO-1 . )»COS(XXTH) 

000002 

00112 

10* 

C FL  CAN  BE  ENTERED  AS  A FUNCTION  OF  PSI 

000002 

00113 

11* 

ONSTEP*NSTEP 

000016 

00114 

12* 

DELTH=-2 . OOEGREE/QNSTEP 

000024 

00115 

13* 

DO  100  N*0, NSTEP 

000150 

00120 

14* 

XN*N 

000150 

00121 

15* 

00122 

16* 

00124 

17* 

00127 

18* 

00131 

19* 

00132 

20* 

00133 

21* 

00134 

22* 

00136 

23* 

00137 

24* 

00141 

25* 

00142 

26* 

00144 

27* 

00145 

28* 

00146 

29* 

00147 

30* 

00150 

31* 

00151 

32* 

00152 

33* 

00153 

34* 

00154 

35* 

00155 

36* 

00156 

37* 

00157 

38* 

00160 

39* 

00161 

40* 

00162 

41* 

00163 

42* 

00164 

43* 

00164 

44* 

00164 

45* 

00165 

46* 

00165 

47* 

00165 

48* 

00166 

49* 

00167 

50* 

00170 

51* 

00171 

52* 

00172 

53* 

00173 

54* 

00174 

55* 

00175 

56* 

00176 

57* 

00177 

58* 

00177 

59* 

00200 

60* 

00200 

61* 

00201 

62* 

00202 

63* 

00203 

64* 

0020': 

65* 

00204 

66* 

00204 

67* 

00205 

68* 

00206 

69* 

00207 

70* 

0021 1 

71* 

PSIDEG*  DEGREE+XN*DELTH 
IF( ITYP.EQ. 1 ) GO  TO  5 
DO  4 I D = 1 ,6 
IF(PSIDEG.GT.THETA)GO  TO  3 
D( ID)=DA( ID , 1 3) 

GO  TO  4 

3 E { I D ) = DA ( ID+6, 13) 

4 CONTINUE 

PSI=(-DEGREE+PSIDEG) *0.01 745 

IF( PSIDEG.GT .THETA)  PSI = ( -DEGREE-PSIDEG) *0 . 01 745 

5 CONTINUE 

IFdTYP  .EQ.1)  PSI  = 0.01745*PSIDEG 
P3  = PSI * Y LAM( 3 ) 

P2  = PSI * Y LAM( 2 ) 

CS=C0S(P3) 

SN=SIN ( P3 ) 

P 1 = PS I * Y L AM ( 1 ) 

CSP=CCSH ( P 1 ) 

SNP  = S I NH (PI ) 

U=D( 1 )*BETA( 1 )*SNP+D(2)*BETA( 1 )*CSP 
PHI =D( 1 )*BETA(2)*SNP+D(2)*BETA(2)*CSP 
W=D( 1 ) *CSP+D( 2 ) *SNP 
PA=BETA( 1 ) *YLAM( 1 )-1 .0 
PB=BETA(2)*YLAM( 1 ) 

EZERO=(PA-PB)*CSP«D( 1 )+(PA-PB)*SNP*D(2) 

EPI=(PA+PB)*CSP*D( 1 )♦( PA+PB)*SNP*D(2) 
OM=-RHO*H(3)*(PA+RHO*PB)*(D( 1 ) *CSP+D( 2 ) *SNP ) 

Qf.'=  ( PA*RHO*H(  1 ) + RHO*H(  3 ) * PB ) * ( D ( 1 ) *CSP+D(2)  *SNP) 
U=U-D(3)*BETA(3)*CS+D(4)*BETA(3)*SN 

1 +D(5)*(BETA(5)*SN-BETA(3)*PSI*CS) 

2 +D(6)*(BETA(5)‘CS+BETA(3).P-.:»SN) 
PHI=PHI-D(3)*BETA(4)*CS+D(4)*BETA(4)*SN 

1 +D(5)*(BETA(6)*SN-BETA(4)*PSI*CS) 

2 +D(6)*(BETA(6)*CS+BETA(4)*PSI*SN) 
W=W+D(3)*SN+D(4)*CS+D(5)*PSI*SN+D(6) * PSI *CS+G( 1 ) 

PA  = BET  A ( 3 ) - 1 . 0-BETA ( 4 ) 

PB=BETA(5)-BETA(3)-BETA(6)+BETA(4) 

PC=BETA( 3 ) — 1 . 0+BETA ( 4 ) 

PD=BETA(5)-BETA(3)+BETA(6)-BETA(4) 

PE= ( BET  A ( 3 )-1 . 0 ) *RHO*H( 1 ) +RHO*H( 3 ) *BET A ( 4 ) 
PF=(BETA(5)-BETA(3) )*RH0*H(1 )+RH0*H(3)*(BETA(6)-BETA(4) ) 
PG=BETA(3)-1 . 0+RHO*BETA(4 ) > 

PH=BETA(5)-BETA(3)+RH0*(BETA(6)-BETA(4) ) 
EZER0=EZER0+D(3)*PA*SN+D(4)*PA»CS-G(1 ) 

1 +D(5)*(PB*CS+PA*PSI*SN)+D(6)*(-PB*SN+PA*PSI*CS) 

EPI  = EPI+D(3)»  PC*SN+D( 4 ) *PC*CS-G ( 1 ) 

1 +D(5)*(PD*CS+PC*PSI*SN)+D(6)*(-PD*SN+PC*PSI*CS) 

EZER0=EZER0/(RH0+1 .0) 

EPI=EPI/(RH0-1 .0) 

0N=0N+D(3)*SN*PE+D(4)*PE*CS+D(5)*(PF*CS+PE*SN*PSI )+D(6)*(-PF 
1*SN+PE*CS*PSI)  — RHO*H( 1 ) *G( 1 ) 

QM=QM-RH0*H(3)*(D(3)*PG*SN+D(4)*PG*CS+D(5)*(PH*CS+PG*SN*PSI) 
1+D(6)*(-PH*SN+PG*CS*PSI)-G(1 )) 

ON=QN/RHO 

OM=QM/RHO 

IF((PSIDEG.LE. THETA). OR. (ITYP.E0.1))  GO  TO  6 
U = -U 


000156 

000161 

000174 

000174 

000175 

000177 

000201 

000204 

000204 

000210 

000222 

000222 

000227 

000232 

000235 

000241 

000245 

000250 

000254 

000260 

000270 

000300 

000306 

000310 

000312 

000325 

000336 

000346 

000356 

000356 

000356 

000407 

000407 

000407 

000440 

000461 

000463 

000465 

000467 

000471 

000473 

000475 

000477 

000501 

000501 

000533 

000533 

000565 

000567 

000571 

000571 

000621 

000621 

000653 

000655 

000660 

000675 


* 4 


00212 

72* 

PHI*— PHI 

000677 

00213 

73* 

6 

CONTINUE 

000702 

00214 

74* 

IF(INTEG.EQ.O)GO  TO  10 

000702 

00216 

75* 

NI*8*N 

000703 

00217 

76* 

S0L(NI+1 )=PSIDEG 

000707 

00220 

77* 

CINT  =WFA*FL 

0007 1 1 

00221 

78* 

S0L(NI+2)=S0l(NI+2l+U*CINT 

0C0713 

00222 

79* 

SOL( NI +3 ) *S0l ( NI»3)*W*CINT 

000716 

00223 

80* 

SOL ( NI+4)*S0L(NI+4)+PHl*CINT 

000722 

00224 

81* 

SOL ( NI+5 )*S0L(NI*5)*EZER0*CINT 

000726 

00225 

82* 

S0L(NI+6)*SCL(NI*-6i*EPI*CINT 

000732 

00226 

83* 

S0L( NI  + 7 ) sSOL ( NI  + 7 ) *ON*C I NT 

000736 

00227 

84* 

S0L( NI+8 )*S0L(NI+8 ) +QM*CINT 

000742 

00230 

85* 

GO  TO  11 

000746 

00231 

86* 

10 

CONTINUE 

000750 

00232 

87* 

PRINT  1001,  PSIDEG.U.W.PHI , EZERO . EPI , QN , QM 

000750 

00244 

88* 

1 1 

CONTINUE 

000764 

00245 

89* 

IF  (ISAV  . EQ.  0)  GO  TO  100 

000764 

00247 

90* 

IF  (ABS(PSIDEG)  .GT.  0.5  .OR.  IND  .EO.  1)  GO  TO  100 

000765 

00251 

91  • 

WRL=ABS(SMN/EZERO)/E 

001002 

00252 

92* 

WRITE  (15,1002)  SMN.W.WRL 

001007 

00257 

93* 

IC=IC+1 

001017 

00260 

94* 

IND*1 

001022 

00261 

95* 

100 

CONTINUE 

001026 

00263 

96* 

RETURN 

001026 

00264 

97* 

1001 

FORMAT  (/F10.3.1P7E17.8) 

001074 

00265 

98* 

1002 

FORMAT  (3E20.8) 

001074 

00266 

99* 

END 

001074 

END  OF  COMPILATION 


NO  DIAGNOSTICS 


MFOR.S  TPFS.BOUNDARYARCH.B 
FOR  00E3-1 1/03/77-16:02:55  (0.) 


SUBROUTINE  CALC A ENTRY  POINT  000250 


STORAGE  USED:  C0DE(1)  000255:  DATA ( 0 ) 000034;  BLANK  C0MMQN(2)  001541 


EXTERNAL  REFERENCES  (BLOCK.  NAME) 


0003 

DCOS 

0004 

DSIN 

0005 

DSINH 

0006 

DCOSH 

0007 

NERR3S 

STORAGE  ASSIGNMENT  (BLOCK.  TYPE.  RELATIVE  LOCATION.  NAME) 


0001 

000056 

1 OL 

0001 

000073 

1 1 L 

0001 

000166 

20L 

0001 

000226 

21  L 

0000 

D 

000000 

ANGLE 

0002 

D 

000007 

ARR 

0000 

D 

000004 

A1 

0000 

D 

000002 

A3 

0002 

0 

000053 

beta 

0002 

D 

000175 

BLAM2 

0002 

D 

00021 1 

B LAM3 

0000 

D 

000010 

CONI 

0000 

D 

000012 

C0N2 

0000 

0 

000006 

C01 

0002 

0 

00011 1 

CSA1 

0002 

D 

000115 

CSA3 

0002 

0 

000161 

D 

0002 

D 

000225 

DA 

0002 

D 

000005 

E 

0002 

0 

G0007S 

G 

0002 

0 

000125 

H 

0002 

D 

000121 

HCA2 

0002 

D 

000123 

HSA2 

0002 

I 

001537 

I BO'JN 

0002 

000004 

IC 

0002 

I 

001540 

IJUNC 

0000 

000020 

INOPS 

0002 

001535 

INTEG 

0002 

001536 

ITYP 

0002 

D 

000000 

RHO 

0002 

D 

000113 

SNA1 

0002 

D 

000117 

SNA3 

0002 

D 

000715 

SOL 

0002 

0 

00b  * 45 

XI 

0002 

D 

c-roooa 

XN 

0002 

D 

000153 

XR 

0002 

D 

000067 

YLAM 

00101 

1* 

SUBROUTINE  CALCA  (DEGREE) 

OOOoC'i 

00103 

2* 

IMPLICIT  DOUBLE  PRECISION(A-H.O-Z) 

000000 

00104 

3* 

COMMON  RHO.XN.IC.E 

ocoooo 

00105 

4* 

COMMON  ARR ( 3 , 6 ) , BET  A ( 6 ) , Y LAM( 3 I ,G(6) . JSA1 . SNA1 . CSA3 , SNA  3 

000000 

00106 

5* 

COMMON  HCA2 , HSA2 , H( B ) 

GCCC90 

00107 

6* 

COMMON  XI (3) .XR(3) .0(6) ,BLAM2(6) , B L/ 43(61 

000000 

00110 

7* 

COMMON  0A( 12. 13) ,SOL(200) , INTEG. ITV° , I BOUN . I JUNC 

*,10000 

001 1 1 

8* 

ANGLE=0 .01 745 ‘DEGREE 

00b  TOO 

00112 

9* 

A3=YLAM(3)*ANGLE 

000002 

00113 

10* 

CSA3=C0S( A3) 

00000-1 

00114 

11* 

SNA3=SIN(A3) 

00001 c 

00115 

12* 

A 1 = Y LAM( 1 ) * ANGLE 

0000  !' 

00116 

13* 

SNAI =S1NH( A1  ) 

0"001 t 

00117 

14* 

CSA1 =COSH( A1 ) 

0000:3 

00120 

15* 

A«R( 1 . 1 ) =BETA( 1 (-SNAI 

000027 

00121 

16* 

ARR( 1 .2)*BETA( 1 )*CSA1 

000032 

00122 

17* 

IF( ( IBOUN.EO. 1 ) .AND. ( IJUNC. EO.O) ) GO  TO  10 

000034 

00124 

18* 

ARR ( 2 , 1 ) * BETA ( 2 ) *SNA  1 

000046 

00125 

19* 

ARR(2.2)*BETA(2)*CSA1 

000051 

00126 

20* 

GO  TO  1 1 

000054 

00127 

21* 

10 

C01 * BET  A ( 1 ) • Y LAM( 1 )-1 . 0«-RH0*  BET  A ( 2 ) *YLAM( 1 ) 

000056 

00130 

22* 

ARR ( 2 . 1 ) =C01 *CSA1 

000065 

00131 

23* 

ARR{2.2)r“0l*SNA1 

000067 

00132 

24* 

11 

CONTINUE 

000073 

00133 

25* 

ARR(3.1 )=CSA1 

000073 

00134 

26* 

ARR(3.2)=SNA1 

000074 

00135 

27* 

ARR  ( 1 ,3)=-BETA(3)*CS.-3 

000076 

00136 

28* 

ARR | 1 , 4 ) = BET  A ( 3 ) • SNA3 

000102 

00137 

29« 

ARR ( 1 , 5)  = 8ET  A( 5 ) *SNA3-BET  A ( 3 )• ANGLE *CS A3 

000105 

00140 

30* 

ARR ( 1 . 6 ) *BET  A ( 5 ) »CSA3*BET  A( 3) • ANGLE *SN A3 

0001 16 

00141 

31* 

c(  ( IB06N.E0. 1 ) . AND.  ( I JJNC.EO.O)  ) GO  TO 

20 

000124 

00143 

32* 

ARR(2,3)=-BETA(4)*CSA3 

000136 

00144 

33* 

ARR(2.4)=BETA(4)*SNA3 

000142 

00145 

34* 

ARR(2,5)=BETA(6)*SNA3-BETA(4) • ANGLE *CS A3 

000145 

00146 

35* 

ARR(2,6)=8ETA(6)*CSA3+BETA(4) • ANGLE *SN A3 

000156 

00147 

36* 

GO  TO  21 

000164 

00150 

37* 

20  CONI =BETA(3)”1 . O+RHO*  BET  A ( 4 ) 

000166 

00151 

38* 

C0N2  = BET A ( 5 ) -BETA(  3 )*-RHO*  ( BETA(  6 )-BET A ( 4 ) 

1) 

000173 

00152 

39* 

ARR( 2 , 3 ) =CON 1 » SNA3 

0CC202 

00153 

40* 

ARR(2,4)=C0N1»CSA3 

000204 

00154 

41* 

ARR ( 2 , 5 ) =C0N2«CSA3+C0N1 *ANGlE*SNA3 

00C207 

00155 

42* 

ARR(2,6)=-C0N2*SNA3+C0N1*ANGLE*CSA3 

000216 

00156 

43* 

21  CONTINUE 

000226 

00)57 

44* 

ARR( 3, 3) =SNA3 

000226 

00160 

45* 

ARR(3,4)*CSA3 

000227 

00161 

46* 

ARR ( 3 , 5 ) = ANGLE *SNA3 

000231 

00162 

47* 

ARR(3,6)=ANGLE»CSA3 

000233 

00163 

48* 

RETURN 

000235 

00164 

49* 

END 

000254 

ENO  OF 

COMPILATION:  NO  DIAGNOSTICS. 

DEFORMATION  OF  PRESSURE  STABILIZED  ARCHES 


CONCENTRATED  LOAD 
(NORMAL) 

SIMPLY  SUPPORTED  ENDS 


FINITE  ELEMENT  CHECK  CASE 

C*  . 0000*N+  .0500 

E*  ,000*N+1 .0000 

HALF  SPAN  ANGLE  * 

90.000 

RADIUS  RATIO  * 

20.000 

ELASTIC 

MODULUS  RATIO  > 

1 .000 

SHEAR 

MODULUS  RATIO  = 

.050 

PRESSURE  PARAMETER  = 

.057 

LOAD  LOCATION  > 

.000 

PSI 

U 

w 

PHI 

-90.000 

-2.03287907-20 

0.00000000 

1 58475030-04 

-72.000 

-3.25139188-04 

-1 .20542675-03 

8.97332786-05 

-54.000 

-8 . 60 732065- U4 

-1 .24301006-03 

-5.05428694-05 

-36.000 

-1.20240372-03 

-7.94067007-06 

-1 .81650596-04 

-IB. 000 

-9.79840769-04 

2.18096021-03 

-2.18952714-04 

.000 

4.74338450-20 

4.50843540-03 

-2.26089141-20 

18.000 

9.79840769-04 

2.18096021-03 

2.18952714-04 

36.000 

1.20240372-03 

-7.94067007-06 

1 .81650596-04 

54.000 

8.60732065-04 

-1 .24301006-03 

5.05428694-05 

72.000 

3.25139188-04 

-1 .20542675-03 

-8.97332786-05 

90.000 

0.00000000 

0.00000000 

-1 .58475030-04 

E(0) 

E(PI) 

N 

M 

-1 .65779550-05 

-1  .65779550-05 

-3.31559099-05 

4.23516474-23 

8.00368673-07 

-3.95025582-05 

-3.76939858-05 

2.01388530-05 

4.60105860-06 

-4.43697882-05 

-3.85436923-05 

2.44701009-05 

-1 .15403683-06 

-3.53225971-05 

-3.56218854-05 

1 .70735891-05 

-2.19261753-05 

-6.91271813-06 

-2.92144648-05 

-7.50203102-06 

-7.90181502-05 

6.26127394-05 

-1 .99483988-05 

-7.07711297-05 

-2.19261753-05 

-6.91271813-06 

-2.92144648-05 

-7.50203102-06 

-1 .15403683-06 

-3.53225971-05 

-3.56218854-05 

1.70735891-05 

4.60105860-06 

-4.43697882-05 

-3.85436923-05 

2.44701009-05 

8.00368673-07 

-3.95025582-05 

-3.76939858-05 

2.01388530-0S 

-1.65779550-05 

-1  .65779550-05 

-3.31559099-05 

8.47032947-23 

DEFORMATION  of  pressure  stabilized  arches 


CONCENTRATED  LOAD 
(NORMAL) 

SIMPLY  SUPPORTED  ENDS 


FINITE  ELEMENT  CHECK  CASE 

C*  . 0000*N+  .0500 

E*  . 000*N+1 .0000 

HALF  SPAN  ANGLE  » 

90.000 

RADIUS  RATIO  « 

20.000 

ELASTIC 

MODULUS  RATIO  > 

1.000 

SHEAR 

MODULUS  RATIO  ■ 

.050 

PRESSURE  PARAMETER  * 

.067 

LOAD  LOCATION  « 

.000 

PS  I 

U 

N 

PHI 

-90.000 

-1.01643954-20 

0.00000000 

1 .36283977-04 

-72.000 

-2.96422217-04 

-1.03989002-03 

7.63284971-05 

-54.000 

-7.72215360-04 

-1 .05442386-03 

-4.59015665-05 

-36.000 

-1.07282052-03 

3.92254658-05 

-1 .60391672-04 

-18.000 

-8.73666662-04 

1.97038106-03 

-1 .93859" 12-04 

.000 

-2.71050543-20 

4.05357628-03 

-2.41201731-20 

18.000 

8.73666662-04 

1 .97038106-03 

1 .93859342-04 

36.000 

1.07282052-03 

3.92254658-OS 

1 .60391672-04 

54.000 

7.72215360-04 

-1.05442386-03 

4.59015665-05 

72.000 

2.96422217-04 

-1 .03989002-03 

-7.63284971-05 

90.000 

0.00000000 

0.00000000 

-1.36283977-04 

E(0) 

E(  PI ) 

N 

M 

-1.65779255-05 

-1 .65779255-05 

-3.31558510-05 

-1.48230766-22 

-1 .55356452-06 

-3.69647999-05 

-3.76325295-05 

1.76945378-05 

1.69830150-06 

-4.11982382-05 

-3.84^68521-05 

2.14348479-05 

-3.00068467-06 

-3.32162718-05 

-3.54610941-05 

1.50983394-05 

-2.06971344-05 

-8.01096887-06 

-2.90254559-05 

-6.33911337-06 

-7. 16352752-05 

5.50548433-05 

-1.97496669-05 

-6.33054190-05 

-2.06971344-05 

-8.01096887-06 

-2.90254559-05 

-6.33911337-06 

-3.00068467-06 

-3.32162718-05 

-3.54610941-05 

1.50983394-05 

1.69830150-06 

-4.11982382-05 

-3.84268521-05 

2.14348479-05 

-1.55356452-06 

-3.69647999-05 

-3.76325295-05 

1.76945378-05 

-1.65779255-05 

-1  .65779255-05 

-3.31558510-05 

1.05879118-22 

DEFORMATION  OF  PRESSURE  STABILIZED  ARCHES 


CONCENTRATED  LOAD 
(NORMAL) 

SIMPLY  SUPPORTED  ENDS 

FINITE  ELEMENT  CHECK  CASE 

C*  . OOOO'N-f  .0500 
E*  . 000*N+1 .0000 


HALF  SPAN  ANGLE  « 90.000 

RADIUS  RATIO  * 30.000 

ELASTIC  MODULUS  RATIO  » 1.000 

SHEAR  MODULUS  RATIO  * .050 

PRESSURE  PARAMETER  » .057 

LOAD  LOCATION  « .000 


W 


PSI 

U 

W 

PHI 

E(  0 ) 

E(  PI ) 

N 

M 

-90.000 

1 .62630326-19 

4.65868121-21 

2.19743349-04 

-1  .65779550-05 

-1 .65779550-05 

-3.31559099-05 

1.41172158-22 

-72.000 

-5.26901584-04 

-2.02705377-03 

1 .23351702-04 

-1.57546747-06 

-3.67041428-05 

-3.76939696-05 

1.75594560-05 

-54.000 

-1 .41243755-03 

-2.11982231-03 

-5.65979756-05 

9.73197916-07 

-4.02033255-05 

-3.85436615-05 

2.05825396-05 

-36.000 

-1 .98626875-03 

-9.85614664-05 

-2.30614024-04 

-1.81803947-06 

-3.43460883-05 

-3.56218429-05 

1 .62595041-05 

-. J.000 

-1 .62621295-03 

3.55609708-03 

-3.08653672-04 

-1.73204319-05 

-1 .18019832-05 

-2.92144149-05 

-2.75845749-06 

.000 

2.57498016-19 

7.38875194-03 

-8.99441353-19 

-8.85435716-05 

7.12593476-05 

-1 .99483463-05 

-7.98792524-05 

18.000 

1.62621295-03 

3.55609708-03 

3.08653672-04 

-1 .73204319-05 

-1 .18019832-05 

-2.92144149-05 

-2.75845749-06 

36.000 

1.98626875-03 

-9.85614664-05 

2.30614024-04 

-1.81803947-06 

-3.43460883-05 

-3.56218429-05 

1 .62595041-05 

54.000 

1 .41243755-03 

-2.11982231-03 

5.65979756-05 

9.73197916-07 

-4.02033255-05 

-3.85436615-05 

2.05825396-05 

72.000 

5.26901584-04 

-2.02705377-03 

-1 .23351702-04 

-1.57546747-06 

-3.67041428-05 

-3.76939696-05 

1.75594560-05 

90.000 

0.00000000 

0.00000000 

-2.19743349-04 

-1 .65779550-05 

-1 .65779550-05 

-3.31559099-05 

2.82344316-23 

DEFORMATION  OF  PRESSURE  STABILIZED  ARCHES 


CONCENTRATED  LOAD 
(NORMAL) 

SIMPLY  SUPPORTED  ENDS 

FINITE  ELEMENT  CHECK  CASE 

C*  . 0000  *N+  .0500 
E*  . 000*N+1 . 0000 


HALF  SPAN  ANGLE  * 90.000 

RADIUS  RATIO  = 30.000 

ELASTIC  MOOULUS  RATIO  = 1.000 

SHEAR  MODULUS  RATIO  > .050 

PRESSURE  PARAMETER  * .067 


LOAD  LOCATION  > 

.000 

PSI 

U 

W 

PHI 

-90.000 

-2.43945489-19 

2.54109884-21 

1 .87444223-04 

-72.000 

-4.74621036-04 

-1.72613997-03 

1 .04025766-04 

-54.000 

-1 .25093116-03 

-1.77526265-03 

-5.06792375-05 

-36.000 

-1 .74888422-03 

-8.60272085-06 

-2.00449914-04 

-18.000 

-1 .43089193-03 

3.17093732-03 

-2.69625203-04 

.000 

-1.34847645-18 

6.56372347-03 

2-65411643-18 

18.000 

1.43089193-03 

3.17093732-03 

2.69625203-04 

36.000 

1 .74888422-03 

-8.60272085-06 

2.00449914-04 

54.000 

1 .25093116-03 

-1 .77526265-03 

5.06792375-05 

72.000 

4.74621036-04 

-1 .72613997-03 

-1 .04025766-04 

90.000 

-2.03287907-20 

0.00000000 

-1 .87444223-04 

E(  0 ) 

E(PI) 

N 

M 

-1.65779255-05 

-1  .65779255-05 

-3.31558510-05 

-3.10578747-22 

-3.83076676-06 

-3.43098749-05 

-3.76325154-05 

1.52353185-05 

-1 .74853405-06 

-3.72704883-05 

-3.84268252-05 

1.77560408-05 

-3.80151175-06 

-3.21318488-05 

-3.54610570-05 

1.41612316-05 

-1.65029362-05 

-1  .24549916-05 

-2.90254123-05 

-2.02340980-06 

-7.97383034-05 

6.23576056-05 

-1 .97496211-05 

-7.10282080-05 

-1.65029362-05 

-1 .24549916-05 

-2.90254123-05 

-2.02340980-06 

-3.80151175-06 

-3.21318488-05 

-3.54610570-05 

1 .41612316-05 

-1.74853405-06 

-3.72704883-05 

-3.84268252-05 

1.77560408-05 

-3.83076676-06 

-3.43098749-05 

-3.76325154-05 

1.52353185-05 

-1.65779255-05 

-1 .65779255-05 

-3.31558510-05 

2.82344316-23 

SYMBOLS 


Aj,  Bj,  Cj 

b 

a 

cii 

c 

c,d 

E 


fi 

Gi 


9i 


* 


9 

Hi 

Ji 

Lj 

!9 

Mx,  My,  Mz 


mx,  my,  mz 


constants  of  integration 
force  vector 
cross-section  radius 
material  stiffness  parameters 

vector  whose  elements  are  the  integration  constants  Cj 
nondimensional  material  stiffnesses 

nondimensional  form  of  e and  the  energy  density  function 
axial  strain 

applied  surface  loading 
applied  force 

generalized  forces  associated  with  Fj 

nondimensional  form  of  fj 

applied  line  load  intensity 

generalized  forces  associated  with  Gj 

nondimensional  form  of  g; 

magnitude  of  concentrated  load 

coefficients  of  the  governing  differential  equations 

constants  used  in  writing  the  boundary  condition  matrix 

constants  defined  in  Appendix  E 

gage  length 

nondimensional  form  of  mx,  , my,  mz 
internal  moments 
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stress  resultants 


Qy.  Qz 

V qz 

R, , R2 


u.  V.  w 


u,  V,  w 


stress  resultants  due  to  pressurization 

first  order  stress  resultants  due  to  applied  load 

second  order  stress  resultants  due  to  applied  load 

nondimensional  pressure  parameter 

pressure 

nondimensional  form  of  qy,  qz 
internal  shear  forces 
principal  radii  of  curvature 
arch  radius 
coefficient  matrix 
elements  of  matrix 
arch  length 

nondimensional  form  of  t 
internal  axial  force 
displacements  due  to  applied  load 
cross-section  displacements 
nondimensional  form  of  U,  V,  W 
general  displacements 
displacements  due  to  pressurization 
constants  resulting  from  integration  w.r.t.  0, 
half  span  angle  of  arch 
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«2 

ft 

rY-  rz 
7y.  7Z 
7 

eij 

*°ii 

cij 

Mi 

Ox.  02 

h 

tx 

Kx,  Ky,  Kz 

KX,  Ky,  KZ 

X 


<t>x.  «y,  «Z 


CO 


parameter  defining  arch  length 

constants  used  in  the  relocations  among  A,,  B,,  and  C-, 

nondimensional  form  of  yz 

transverse  shear  strains 

nondimensional  flexibility 

strain  components 

strain  components  due  to  pressurization 

first  order  strain  components  due  to  applied  load 

second  order  strain  components  due  to  applied  load 

coordinates  locating  position  on  arch 

location  of  applied  line  load 

location  of  concentrated  load 

nondimensional  forms  of  xx,  Ky,  k2 

curvatures 

parameter  specifying  the  solution  of  the  characteristic 
equation 

radius  ratio 

cross-section  rotations 

angle  defining  region  over  which  load  acts  and  characteristic 
number 


Subscript  w denotes  value  parameter  at  wrinkling 
( )'  denotes  differentiation  w.r.t.  02  except  for  ey  and  Njj 

6 denotes  variation 
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